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FOREWORD 

The prospect of undertaking a reusable launch vehicle development led the NASA 
Office of Manned Space Flight (OM8F) to rec^a st the Office of Advanced Research and 
Technology (OART) to organize and direct a program to develop the technology that 
would aid in selecting the best system alternatives and that would support the ultimate 
development of an earth -to -orbit shuttle. Such a Space Transportation System Tech- 
nology Program has been initiated. OART, OMSF, and NASA Flight and Research 
Centers with the considerable inputs of Department of Defense personnel have genera- 
ted the program through the efforts of several Technology Working Groups and a Tech- 
nology Steering Group. Funding and management of the recommended efforts is being 
accomplished through the normal OART and OMSF line management channels. The 
work is being done in government laboratories and under contract with industry and 
universities. Foreign nations have been invited to participate in this work as well. 
Substantial funding, from both OART and OMSF, was applied during the second half of 
fiscal year 1970. 

The Space Transportation System Technology Symposium held at the NASA Lewis 
Research Center, Cleveland, Ohio, July 15-17, 1970, was the first public report on 
that program. The Symposium goals were to consider the technology problems, their 
status, and the prospective program outlook for the benefit of the industry, govern- 
ment, university, and foreign participants considered to be contributors to the pro- 
gram. In addition, it offered an opportunity to identify the responsible-individuals al- 
ready engaged in the program. The Symposium sessions were intended to confront 
each presenter with his technical peers as listeners, and this, I believe, was substan- 
tially accomplished. 

Because of the high interest in the material presented, and also because the people 
who could edit the output are already deeply involved in other important tasks, we have 
elected to publish the material essentially as it was presented, utilizing mainly the il- 
lustrations used by the presenters along with brief words of explanation. Those who 
heard the presentations, and those who are technically astute in specialty areas, can 
probably put this story together again. Wc hope that more will be gained by compil- 
ing the information in this form now than by spending the time and effort to publish 
a more finished compendium later. 


A. O. Tischler 
Chairman, 

Space Transportation System 
Technology Steering Group 
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Dynamics and Aeroelastlcity 
Introductory Statement 
by Harry L. Runyan 

The Space Shuttle, being an hybrid - an airplane and a launch vehicle - 
represents the greatest challenge that the dynamioist and the aeroelasti- 
cian have faced. Some specific problem areas related to the Space Shuttle 
are listed on figure 1. 

Dynamics and aeroelasticity envelop many disciplines, including aero- 
dynamics, vibration, random processes, structures, fluid flow, mechanics, 
etc., but, of more importance, they involve the interaction and coupling 
of many of these various disciplines. Fundamentally, we are concerned 
with structural integrity and safe flight, i.e., trying to ensure that 
the vehicle will remain structurally intact as well as function properly 
in the presence of the many faceted dynamic environment. 

A new area which may have an impact on our task is the effect of the 
high temperature environment. In the past, we have been able to success- 
fully decouple the temperature effects from our problem formulation. For 
the Space Shuttle, this problem must be closely reexamined. 

As has bsen pointed out in the opening remarks of the Conference, the 
Dynamics and -Aeroelasticity Technology Group comprises one of several tech- 
nology groups which are attempting to provide the necessary research to 
support a successful and safe vehicle. The group has members from most of the 
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NASA Centers as well as from the Air Force groups. We meet periodically 
vo review ongoing work, search for new problem areas; and we are constantly 
updating and revising our program. The group la organised Into three 
panels as shown on figure 2: a panel on Dynamic Loads and Response, one 

on Aeroelastlclty, and one on Flight Dynamics and Environment. The 
Conference papers accordingly have been grouped In these same three areas, 
with each Panel Chairman acting as moderator for his particular resslon. 
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TECHNOLOGY AREA 
DYNAMICS AND AEROELAST1C 1TY 
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The direct solution procedure is a vital part of other eigenvalue routines (Assumed Functions. Matrix 
Reduction, and Component Modes). It is a standard subroutine in most large-scale structural analysis 
programs (NASTRAN, STARDYNE, ASTRA, REXBAT, SNAP, FORMAT, etc.). 
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FIXED-WING QRBITER (FIRST MODE) 








2n functions when one wants n eigenvalues. 


ASSUMED FUNCTION 
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APPLICATION OF DYNAMIC MODELS 
by 

H. Wayne Leonard and Homer G. Morgan 
NASA-Langley Research Center 

SUMMARY 

Structural dynamic, models have been widely used to obtain vibration 
data required in launch vehicle development. A generalized model is now 
ing used to obtain lateral response data on space shuttle configurations 
to account for new factors such as parallel staging, lifting surfaces, and 
interface stiffness. A second generation generalized model that would 
have representative longitudinal responses is also believed to be necessary. 
Finally, a detailed nearly-repl ica dynamic model of the final space shuttle 
configuration is proposed as the most cost-effective technique for obtain- 
ing timely vibration data to establish confidence in the design. Such a 
model would be an integral part of the development program, and would be 
closely coordinated with analysis and limited full-scale vibration tests. 
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INTRODUCTION 

Many kinds of dynamic models have been used to solve special problems 
in the development of airplanes and spacecraft. This paper is concerned 
with one particular type of dynamic model - the subscale model for predict- 
ing vibration properties and structural dynamic response of a complete 
vehicle. Slide 1 shows models of this type that have been used in recent 
launch vehicle development programs. These models of Saturn I, Titan III, 
and Saturn V have demonstrated their versatility, usefulness, and cost 
effectiveness for obtaining the modal vibration data that are required to 
verify analytical methods used in response and stability studies. Each 
of these models employed different degrees of sophistication in their 
design and construction, arf each was related to the development program 
in a different way, as discussed in references 1 through h. 

This background of structural dynamic models technology can be 
utilized in the space shuttle program to develop both new technology and 
the vehicle itself. However, such models can only be utilized effectively 
when closely coordinated with analysis. In this paper, a series of 
structural dynamic models is described vhich will complement the develop- 
ment of analytical methods for the space shuttle program. The series of 
models would culminate in a near-replica model that would provide most 
of the experimental data that have usually been obtained on full-scale 
test vehicles. 


^8l8h 


30 








FIRST GENERATION MODEL 


The complexity of dynamic models in use at a given time must reflect 
the definition of the space shuttle system at that time. The accompanying 
analysis will generally be no more complex then the model to which it is 
applied. At the present time, very few definitive details are known about 
the final space shuttle configuration. Thus, any model at this time must 
be a generalized structure for studying classes of problems rather than 
for generating quantitative data for specific application. 

A first generation generalized model has been built, and is shown in 
Slide 2. Now being prepared for testing at Langley Research Center, it 
will be used in conjunction with analysis to study the coupled dynamic 
response of winged paral lel-staged bodies. The interface connection spring 
rate will be varied parametrically over a wide range. The model is approxi 
mately 1/15-scale, based on 6-mohth old mass and stiffness data from a 
Manned Spacecraft Center 12,500-pound payload vehicle. However, it is 
generalized in that the orbiter can be mounted forward or aft on the 
booster, as illustrated here, and both orbiter and booster can have either 
straight wings (as 'shown) or delta wings. The model contains no liquids, 
but propellant mass and irertia are simulated by ballast masses to obtain 
liftoff, maximum dynamic pressure, and booster burnout weight conditions. 
Since the parallel bodies are scaled as beams, this model will exhibit repre 
sentative coupled lateral vibrations of importance in flight loads and 
control system analyses. However, its longitudinal response will not be 
representative of an actual vehicle. 
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SECOND GENERATION M3DEL 

As the space shuttle configuration becomes better defined, an advanced 
second generation structural dynamics model, such as shown schematically 
in Slide 3, will be used to study coupled behavior of a more representative 
vehicle. Such a model will be of 1/10- to 1/5-scale in size, and can be 
buiU very early in Phase C of space shuttle development. Stiffness and 
mass distributions will be refined, and 1 iqu id -carrying tankage will be 
included to obtain the first meaningful experimental data on coupled 
longitudinal responses. Thus, early data on the structural part of the 
Pogo problem for the space shuttle configuration could be obtained. To 
be useful for this purpose, the model should have representative liquid- 
shell interactions, propellant feedline characteristics, and thrust struc- 
ture details. Other configuration features such as the interface connec- 
tion, lifting surface geometry, and relative position of the orbiter and 
booster would be incorporated in this model. 



ADVANCED DYNAMICS MODEL FOR SPACE SHUTTLE 


MASS-STIFFNESS DISTRIBUTION 
LONGITUDINAL DYNAMICS 
POGO 

LIQUID-SHELL INTERACTIONS 
PROPELLANT-FEEDLINE DYNAMICS 
THRUST STRUCTURE DETAILS 

INTERFACE CONNECTION 

LIFTING SURFACE GEOMETRY 

BOOSTER-ORBITER POSITIONING 




THIRD GENERATION M3DEL 


First and second generation models will be used essentially to develop 
structural dynamics technology for space shuttle configurations. However, 
the primary purpose of the final or third generation model would be to 
directly support development and testing of the selected vehicle. Design 
and fabrication of this model should parallel final design, should lead 
fabrication of the prototype, and could be most efficiently performed by 
the prototype contractor. It would be a subscale replica of the prototype 
structure within practical limits. Final scale selection would depend on 
tradeoffs involving available facilities, ease of construction, testing 
convenience, and cost, but would probably be in the range of 1/3- to 1/5- 
scale. Some advancements in modeling technology may be required to pro- 
perly account for the thermal protective system. 

The role of such a model in the development of the space shuttle is 
illustrated in Slide 4. Model tests, full-scale tests, and flight opera- 
tions are shown to be unified by a constantly refined analytical model of 
the vehicle. The time sequence of tests is roughly indicated from left 
to right. Early test and analysis of the model, sitting airplane-fashion 
on its wheels, followed by test and analysis of the mated vehicle would 
lead to a verified analytical model prior to the availability of full- 
scale hardware. Some full-scale airplane-type ground vibration tests 
on flight hardware would then be conducted as checks and to develop more 
confidence in the analysis. 

The complexity and importance in dynamic response of the structural 
interface between the orb iter and the booster will require some limited 
testing on full-scale hardware to check both thr model and analysis. . 
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MODEL APPROACH 








These tests could be static deflection measurements or vibration tests on 
the mated vehicles. Such tests could be done in the empty weight condition 
with the vehicle horizontal, if the configuration permits, or cantilevered 

from the launch tower. 

The primary data source for supporting flight operations would be 
analysis. However, the model, being more inexpensive tc maintain on 
standby than a full-scale dynamic test vehicle, would become the primary 
tool for evaluation of structural modifications and payload changes, and 
for investigation of operational problems. 

RECOMMENDED DYNAMIC M3DEL PRO CHAM 

Slide 5 presents recommendations for the utilization of a detailed 
structural dynamic model to achieve the most cost-effective and reliable 
space shuttle. First, the model must be an integral part of the develop- 
ment program. Then, it must be closely coordinated with analysis and 
full-scale tests. If these factors are achieved, then only limited full- 
scale testing is required. Analysis, thoroughly verified by model tests 
and further checked by some full-scale testing, would provide the final 
data for all flight decisions. However, since configurations change and 
maintaining a model on standby in -datively inexpensive, the model should 
be kept current so that it c*n be used throughout the program for evaluat- 
ing configuration or payload changes, and for operational problem solving. 
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RECOMMENDATION S 

REPLICA MODEL TO BE INTEGRAL PART OF DEVELOPMENT PROGRAM 
COORDINATED MODEL— ANALYSIS— FULL-SCALE PROGRAM 
LIMITED FULL-SCALE TESTS REQUIRED 

VERIFIED ANALYSIS PROVIDES FINAL DATA FOR FLIGHT DECISIONS 

MODEL CONFIGURATION KEPT CURRENT AND USED FOR EVALUATION OF 
CONFIGURATION OR PAYLOAD CHANGES AND FOR OPERATIONAL 
PROBLEM SOLVING 
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ALTERNATIVE TO DYNAMIC MODELS 


The alternative to an integrated dynamic model program in space shuttle 
development lends to the situation described in Slide 6. A more extensive 
full-scale test program would be required, This would probably require a 
special dynamic test vehicle and a new test facility. An order-of-magnitude 
increase in manpower is required. Test operations and logistics also 
become more complex so that the overall cost-effectiveness of this approach 
is questionable. 

A significant factor this approach is that data generated wauld be 
available late in the development cycle. Any changes to design resulting 
from the tests could only be implemented with considerable delay to the 
total program. 

The scope and complexity of an extensive full-scale test program 
necessitates detailed long-range planning with a resultant loss of flexi- 
bility. The rate of data generation is also slower than in a model program 
due tc che time required to effect configuration or test condition changes. 


ALTERNATIVE TO MODEL APPROACH 


EXTENSIVE FULL-SCALE TEST PROGRAM REQUIRED 

SPECIAL TEST VEHICLE 
LARGE TEST FACILITY 
INCREASED MANPOWER 
COMPLEX TEST OPERATIONS 
INCREASED COST 


LOSS OF TIMELINESS 

DATA AVAILABLE LATE IN PROGRAM 
LESS IMPACT ON DESIGN 


LESS FLEXIBILITY IN TEST PROGRAM 

COMPLEXITY OF OPERATION 
NUMBER OF PEOPLE INVOLVED 
DIFFICULTY OF MODIFICATIONS 
LONG LEAD TIME PLANNING 
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CONCLUDING REMARKS 


Structural dynamic models have been successfully applied in previous 
launch vehicle development programs to obtain vibration data. Models can 
also be used early in space shuttle development to quickly and inexpensive- 
ly obtain the data required for verifying analytical methods for treating 
new configurations. 

A detailed structural dynamic model of the final space shuttle vehicle 
can also be used effectively and inexpensively to provide most of the 
experimental data for verifying analytical design methods. Only limited 
testing of full-scale flight hardware vro u Id be needed to supplement model 
data. The alternative to a dynamic model approach, without sacrifice of 
confidence in the design, is extensive, expensive, and less timely full- 
scale testing. 
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BASEL VIBRATION AND RANDOM LOADS 

by William C. Walton, Jr. N 7 Q - g ^ q 

and * * 8 

Eugene c. Neumann 

HftSA-Langley Research Center 
Hampton, Virginia 


RESEARCH OBJECTIVES FOR TPS PANELS 


We are working to establish procedures to predict basic vibration 
characteristics and random response of TPS panels. The purpose of this 
talk is to tell you what we are doing and why. Please consider my 
r emar ks as an invitation to constructive criticism on our efforts. 

I emphasize that the work on vibration characteristics is general 
and not restricted to providing Inputs to our own random response 
studies. For example, we are-committed to provide structural inputs to 

support the panel flutter work which will later be discussed by Dr. 

Dixon of the Structures Research Division. 


RESEARCH OBJECTIVES FOR TPS PANELS 

CAPABILITY TO PREDICT BASIC-VIBRATION- 
CHARACTERISTICS 

CAPABILITY TO PREDICT RESPONSE TO 
FLUCTUATING PRESSURE LOADS 
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SINGLE PANEL SYSTEM 


This slide shows two things. It shows our delineation of the panel 
vibration problem which is necessarily somewhat arbitrary 
think logical. It shows also what we presently see as the possibly 
significant effects to be considered in characterizing the TPS skin o 
vibrations . We have decided to concentrate for the time beingonwhat 
we cal 1 a single panel system - defined as a portion of elastic skin 
bounded by thermal expansion joints and resting on elastic standoff 
fixtures which in turn rest on a rigid underlying structure. Whjt are 
our reasons for defining the problem within th ese particular lindts. 

First, like anyone else, we are motivated by the desire o ™ . 

simple and reasonably sized piece of structure as we can, consistent 
with retaining essential engineering information. Second, since the 
purpose of thlrmal expansion Joints is to statically isolate a portion of 
skin from surrounding skin, it is plausible that such joints will de— 
cS^brrtS Twell.' Thirdly, the TPS skinandataadoffmtures 
carry only local pressure loads and, therefore, should be relatively light 
stStSe She structure vbtcC carries integrated l«.de 

generally be relatively stiff making the assumption of * ri B* d the 

reasonable. Finally, even if it should turn out that vibrations of the 
s?S” pa^l Astern do couple significantly to surrounding skin or under- 
lying structure, an understanding of the single- panel system will atill 
be necessary before the behavior of the coupled system can be intelli- 
gently assessed. 

Now I will go over all the effects which we believe should be 
considered. We know we will be dealing with very thin skins - down 
around .01 inch. This means that random local sha^ im^rfections may^ 
be larger than the skin thickness, substantially affecting local stiff 
ness ai stresses. This has already caused us problems in “ttemptato 
measure vibration modes of very thin flat plates ao a subsequent slide 
Jill show. Next, we are going to have to deal with skin conf 1 ^®^" 8 
and support fixtureo which have complicated geometries. Typical candi- 
date configurations feature corrugated skins. We have Been P r °P° B d 
standoff fixtures ranging from simple posts to *™ 8B h t 

wobs with cutouts. The different types of standoff fixtures a * 

variously located points on the skin, sometimes on the edge 1 ^tJ£Js 
in the interior. Proposed methods for securing ek*n to standoff fixtures 
have Included fixed attachments, such as B pot welds, and sliding it** 6 * - 
£".!Tlcb nUo»-moic. free the™ol cxpcnclon. 

tiiin walled stiffeners running perpendicular to the corrugawlono* Major 
effects which it will 6train our capacities to understand and represent, 
rtem fJ™ fact that these panels will get hot, aom^pnr^giL In 
excess of 2000°F. We shall have to deal, therefore, with the material 
changes^ euch^M lowered modulus of elasticity which occur at elevated 
temueiatures and rerheps with permanent material changes caused by 

XS It 1. true that .very effort .111 rte-Ksn 
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panels which can expand so as to relieve thermal stress hut we cannot 
assume that such designs will he completely * h * 8e 

nanels will also he affected hy the presence of a static pressure 
differential. It has been proposed to reduce pressure differentials 
5»riSr atmospheric flight by introducing a gas flow in the space between 

the TPS skin and the primary structure “ an f a^PS^kin^? 3 ^ 6 

aaualizins system may he necessary in order to allow- a TPS skin oi 
sufficiently light weight which can still endure the steady aerodynamic 
loSs in the atmospheric portions of the Shuttle trajectory. However, 
onmniete eaualization will certainly not he realized and the prestress 
the remaining pressure differential could profoundly affect vibra- 
tion response. 


APPROACH FOR BASIC VIBRATION CHARACTERISTICS 


This slide states the two main ideas in our approach as Presently 

“SSS =u. 

structural description of the panel To ^ ’ wh l ch can he attached 

We will measure success in terms of the conriaenc of 

to these modes which we determine. I emphasize inaix me 
vibration while the system is in a e^^^ pres ^ p re; iict 

accomplish our objective^ we require supporting capab * * b a V e 

static stress Cue to MoM psrtlculsrir 

within a panel. We want historielTand find 

the response analyst can add up •> other things vr hope to provide 

nil noir.tB of significant, stress. Among other tmngn w '* r~ r 

in this way a more rational basis for fatigue life predic on. 

With regard to structural damping, we recognize that In addition 

jr sg rr 

damping plays. 
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APPROACH FOR BASIC VIBRATION 
CHARACTERISTICS 

STRUCTURAL DESCRIPTION BY VIBRATION MODES 
DAMPING DESCRIPTION BY MODAL DAMPING 


SLIDE 3 



PROJECTS FOR VIBRATION TEST 
CAPABILITY 


DISPLACEMENT AND TEMPERATURE PROBE 
HIGH TEMPERATURE STRAIN GAUGES 
HOLOGRAPH 

SHAKE TEST EQUIPMENT FOR HOT PANELS 


SLIDE 4 
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SCHEMATIC REPRESENTATION OF TEST SETUP 


Let me emphasize that this representation is schematic and for the 
purposes of illustration. Some items of equipment, particularly heat 
shields are not shown. This is a half section viev with the cut along 
the shaded portions. Shown at the top is a quartz tube radiant heater 
capable of sustaining temperatures at the panel specimen here of up to 
1200°F. This frame here represents a hydraulically actuated clamping 
system vhich will accommodate panels up to two feet square. The unheated 
panel is placed in this support system with the clamps not actuated. 

Heat iB applied and the panel brought to steady temperature. Clamps may 
then be applied along any combination of edges. It is felt that with 
this system ve can simulate conditions for TPS non-ablative panels which 
will get hot but which according to present thinking will be free to 
expand. This is a l—l/2 lb force electrodynamic shaker mounted through 
soft springs upon a pneumatic lift for positioning. There is a vacuum 
cup attachment connecting the force rod to the panel and there is a 
force gauge in the force rod. The vacuum attachment and pneumatic lift 
allow quick disengagement of the shaker for damping decay measurements. 

The shaking point can be located anywhere in the middle 75 percent of the 
panel. This is the dual sensor probe which is servo-controlled to maintain 
a constant distance of the tip from the mean position of the vibrating 
panel. This box contains a remotely controlled motor which moves the probe 
along this track. The track upon which the probe rusts in turn moves in 
the perpendicular direction along this track beneath. The second track 
rests on a turntable allowing the entire system to be rotated. These 
movements allow scanning of over 90 percent of the panel uurface. 

Actually there is a second dual sensor probe on a second track not shown. 
The two moving probes are operated simultaneously for last scanning. It 
is believed that it will be possible to moke a complete scan in about five 
minutes. In anticipation of eventual higher temperature requirements the 
support system, probe tracking systems and shaker system are designed for 
modal temperatures of 2000°F. 
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NASTRAN ANALYSIS 


For the analytical aide of the Job, we expect to use the NASTRAN 
computer program - supporting it with other programs If necessary, nils 
slide illustrates an application of NASTRAN which we are attempting in 
order to ascertain the capabilities of the program to predict vibration 
irociee and give accurate and detailed descriptions of the modal stress 

distributions. The panel is corrugated and has two channel stiffeners. 

Our initial goal is to compute the free-body modes', assuming no prestress. 
We take advantage of the four quadrants of symmetry so that we only have 
to deal with a quarter of the panel. We make a very fine subdivision of 
the panel into elements, from the element library of the NASTRAN program. 
Note that even the floor of the channel stiffener is subdivided. 

With a six— degree— of— freedom mass at each node, the quarter panel 
system has about ^600 degrees of freedom. Utilizing the capabilities of 
the NASTRAN program to provide advance estimates of computing time, we 
have determined that a direct attack on the J600 x 3600 eigenvalue problem 
would require from about 10 to SO hours on our CDC t>600 computer to get the 
first mode. 'The time would very probably be closer to 10 hours than to 
20. This amount of required time we consider to be prohibitive, but for 
two reasons we do not consider it discouraging. First, the time reflects 
the fact that the present veroiem of NASTHAN uoea double precision 
arithmetic. Double precision is probably not necessary on the CDC. 

Sinple precision on this machine is almost the equivalent of double 
precision on most machines. Steps are being taken to provide a single 
precision version of NASTRAN and this alone should reduce the time by a 
factor of four or more. Secondly, we feel there is a good chance we can 
reduce the time substantially by Guyan reduction techniques and we are 
positive we can reduce it drastically by substructuring. We have a high 
priority in-house program underway exploring both these avenues. 
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RESPONSE TO FLUCTUATING PRESSURE LOADS 


I’m going to switch now from basic vibrations to response to fluctur- 
ating pressure loads* Charles Coe of Ames is going to speak on the 
aerodynamic aspects of fluctuating pressure loads* My purpose is to state 
our needs, as we presently see them, with regard to load description for 
predicting vibration stress response* 1 think it is fair to say that 
engineering techniques for predicting panel response to flight fluctuating 
pressure loads have been profoundly influenced by the assumption that 
there is always extremely limited information about the actual load 
environment and nearly always limited information about the structure* We 
favor an attempt to break out of this pattern* As I have indicated, at LRC 
we are laying the groundwork for a structural description of the TPS 
surface which which should allow prediction of the details of vibration 
response. It would be most gratifying to see a wind tunnel test program 
aimed at predicting with comparable detail the fluctuating pressure fields 
which will occur throughout the flight trajectory. Assuming adequate 
scaling lows, we think we could make effective use of wind tunnel pressure 
measurements in which a dense array of pressure gauges is packed within 
a scaled panel dimension giving details of the distribution of pressure 
in space as well as time. Utilizing such measured pressures and the 
computed panel modes, we would compute directly the generalized forces 
and from them the panel time history of response. From such response 
time histories corresponding to a number of points in the trajectory, we 
should be able to count stress cycles and from them make rations! estimates 
of fatigue life. 



RESPONSE TO FLUCTUATING PRESSURE 
LOADS 

LIMITED UNDERSTANDING OF TURBULENT FLOWS 
MORE SEVERE TURBULENCE 
THINNER SKINS 
REPEATED FLIGHTS 


SLIDE 7 
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VIBRATIOII OF A SQUARE CLAMPED PLATE 


As indicated in the discussion of a previous slide, we have some 
data which illustrate the effe' • of shape imperfections on vibrations of 
sauare panels. FJ^tted here are ratios of experimentally and analytically 
determined natural frequencies for 2^inch by 2^-inch flat aluminum alloy 
plates with several values of plate thickness. The plates we^e clamped 
on all four edges. The horizontal axis gives the plate thickness as a 
multiple of the thickness >f the thinnest plate. Here beneath is noted 
the corresponding approximate ratios of thickness to edge length. The 
analytical results were obtained using the NASl'RAN computer program. 
Comparisons were made for the first six modes of each plate. The “operate 
errors for the thickest panel are attributed to observed boundary motions 
due to the fact that the plate was quite heavy, neglecting the thickest 
plate, the agreement passes from very good to unacceptable as the plates 
are made thinner and thinner. This same trend was observed for free- 
body vibration tests of the same panels. For each plate measurements 
were made of the deviations from flatness. The maximum deviations 
increased from 20 percent of the thickness for the 2nd thickest plate to 
approximately 2.5 times the thickness for the thinnest plate. It is 
noted that for the two thinnest plates, which had the largest percentage 
deviations from flatness, the first mode exhibits the greatest error. 

This is consistent with the observed long-wave-length character of the 
imperfections. 
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VIBRATION OF A SQUARE CLAMPED PLATE 
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PROPELLANT DYNAMICS PROBLEMS IN SPACE SHUTTLE VEHICLES 

H. Norman Abramson, Franklin T. Dodge, and Daniel D. Kana 
Southwest Research Institute 
San Antonio, Texas 


ABSTRACT 


The dynamic behavior of liquid propellants in launch vehicles and spacecraft 
has been of great design importance for past and current generations of 
vehicles, and will be of equal or even greater importance for future space 
shuttle vehicles. Large masses of propellants, possibly contained in 
nonaxisymmetric tanks and subjected to many varying forms and levels of 
acceleration, may lead to slosh forces of such magnitude as to require 
extensive baffling or may lead directly to other important dynamic loads 
by interaction or coupling with the elastic structure. 7hi3 paper will out* 
line some of these problems and indicate the directions in which preliminary 
efforts are being made to delineate the problems in more detail or to effect 
their solution. 


INTRODUCTION 

Seventy to eighty percent of the lift-off mass of proposed space 
shuttle vehicles consists of cryogenic liquids; the "orbiter" vehicle itself 
will cai ry a larger amount of liquid into orbit than has any previous vehicle. 

It is therefore evident that the dynamic loads exerted by the contained liquids, 
which were of much importance in the design of present-generation vehicles, 
will be even more important in the design of the space shuttle. In fact, the 
nonsymmetric configuration of most proposed shuttle vehicle structures, the 
requirement of airplane-like maneuverability, and the stringent requirements 
imposed on structures and materials by a planned 100 re-use cycle all 
emphasize the need to understand and design for the influence of liquid 



propellant dynamics on shuttle performance. Southwest Research Institute, 
under contract to NASA-Langlty Research Center*, is presently attempting 
to delineate and evaluate some of these liquid dynamics loads and liquid- 
structure interactions that may be important in the design of the shuttle 
vehicle. These problems are discussed briefly in this paper. 

We will first review some of the liquid dynamics problems that 
appear to be of consequence with respect to loads, mention some of the 
means of controlling those loads, and then discuss some preliminary plans 
for studying liquid-structure interaction problems peculiar to shuttle 
vehicles. We will specifically omit discussing pogo phenomena and prob- 
lems comprising what has come to be known as "propellant management in 
orbit. " Pogo problems applicable to space shuttles have been outlined by 
Runyan [1] who pointed out that pogo oscillations cannot be tolerated because 
of the magnitude of the dynamic loads exerted on the vehicles as a result 
of the large offset masses. Certain pump and feed line phenomena, such 
as water hammer and bubble collection in the impellers and inducers, are 
known to couple with pogo instabilities. Propellant management in orbit 
encompasses such problems as (1) "recapturing" the liquid propellant 
after an extended period of zero-g coasting and repositioning at least some 
minimum amount of propellant over the feed line intakes prior to engine 
restarting, and (2) ensuring that only vapor is released during venting 
operations meant to relieve excessive buildup of pressure in the tanks 

♦Contract NASI -9890. 
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through boil-off. Needless to say, all these problems must be recognized 
and accounted for in the design and operation of the shuttle; in some cases, 
substantial developments in technology will be required. 

LOADS ASSOCIATED WITH LIQUID PROPELLANT MOTIONS 

In the present part of this paper, we will describe a number of 
general problems associated with liquid motions, assuming that the struc- 
tures and tankage containing the liquid are rigid. The arrangement of the 
discussion follows the general sequence of events during a mission. 
Propellant Slosh During Ascent 

Although liquid propellant sloshing has been studied extensively 
in the past [2-5], the nonsymmetric geometry of the space shuttle intro- 
duces several new facets to the problem. 

First of all, a side-by-side arrangement of the booster and orbiter 
causes the roll mode of the vehicle and the propellant slosh to be coupled; 
that is, the slosh can induce roll motions of the vehicle, and vice versa. 
Thus, large, resonant slosh torques should be expected at several excita- 
tion frequencies that depend upon the propellant tank diameters and the 
thrust acceleration. Large, rigid-body-like torques, because of the offset 
of the liquid masses from the roll axis, also should be expected. This 
effect was only of minor importance with respect to the outboard tanks of 
tlig Saturn I vehicle and was handled quite readily; however, in the space 
shuttle, such torques will be very large and difficult to control. 


In some designs, the tanks themselves are not symmetric, and an 
orbiter-tanker probably will have compartmented payload tanks, thus 
creating further asymmetry. Nonaxi symmetric tanks will have slosh cou- 
pling to roll excitation about their own axes in addition to that caused by 
offsetting the tank from the vehicle roll axis. Consequently, shuttle vehicle 
propellant sloshing will be susceptable to roll, yaw, and pitch excitations, 
which is a further complication in booster design over past experience [6]. 
Equivalent mechanical models (mathematical) of the slosh characteristics, 
as commonly employed in guidance and control analyses, will be complicated 
by the fact that several slosh modes (having different resonant frequencies) 
can be excited, depending upon the orientation of the excitation to the tank. 
Figure 1 shows a typical -slosh model for a tank chosen to fit the structural 
envelope of the shuttle; this mechanical model contains two "slosh masses ' 
(i.e., two slosh modes are important), a rigidly attached mass m c which 
simulates the rigid body inertia of the liquid to lateral excitations, and two 
mass moment-of-inertias I D and J e which simulate the rigid body reactions 
to pitch, roll, and yaw. 

Another problem caused by asymmetrical conditions, and one which 
can occur even with an axisymmetrical tank and an in-line booster and 
orbiter, arises during maneuvering whenever the thrust axis is not aligned 
with the tank axis. For this condition, the steady- state free surface of 
the liquid is not perpendicular to t ie tank axis; this, in effect, creates 
asymmetry and therefore will bring additional slosh modes into being. 
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Many of these problems are being investigated at SwRI under sponsor- 
ship of NASA -Langley Research Center and NASA-Marshall Space Flight 
Center*. In addition to analytical studies, slosh force -frequency response 
characteristics will be measured in the laboratory for several typical non- 
axisymmetric tank geometries in order to provide information comparable 
to what has been made available previously for symmet-ic tank slosh. 

Liquid Motions Induced at Separation 

Rather large, transverse, impulsive forces will be required to 
separate the booster and the orbiter in the side-by-side configuration. 
Assuming that the thrust level for the orbiter is small during and just 
after separation, transverse impulses will cause gross liquid motions 
and large disturbing forces and torques at a critical period in the control 
of the orbiter. It is probable that some of the propellants will impact on 
tank bulkheads. In any event, significant residual liquid motions will 
remain when thrusting begins and it important that these not be amplified. 

Liquid Behavior in Orbit and During Docking 

Besides those propellant management problems mentioned in the 
Introduction, other fluid motions during orbital operations can cause large 
loads on the vehicle structure and exert disturbing effects on the guidance 
and control system. There may be a large mass of liquid onboard at this 
time, especially when the orbiter is used as a tanker so that it contains 
liquids other than its own propellants. 

^Contract NAS8- 25920. 
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It is likely that the orhiter will "stop" several times, each time 
docking with some other vehicle. When the payload or propellant tanks are 
partially empty, the docking deceleration has the potential to move the 
liquid bodily from one end of a tank to the other, and the impact loads that 
result might cause "de-latching" (i.e., an aborted docking) or perhaps 
even damage the tank bulkheads. On the other hand, the impact loads might 
be used to advantage to eliminate the rocket thrust requirement for com- 
pletion of a "hard" docking. But the fluid mechanics of docking have so far 
been almost entirely neglected, especially as regards analysis, and conse- 
quently our knowledge of liquid motions during docking is meager . * It is 
believed that the motions that can occur are similar to those shown in 
Figure 2. The nonsymmetrical type of motion shown in Figure 2c is the 
most probable, for free-surface stability reasons, unless extreme care is 
taken to insure that the docking is aligned perfectly with the tank axes, 
these nonsymmetric fluid motions obviously impose the largest delatching 
torques. It is known, also, that a certain minimum deceleration is required 
to disrupt the steady free surface before the liquid can move bodily. Thus, 
if the docking deceleration could be maintained below this value then only 
small disturbances would occur; unfortunately, the destabilizing accelera- 
tion is very small. For example, the acceleration that will destabilize 
a zero-g interface with a zero degree contact angle in a cylindrical tank is 

♦A* study of this problem is being sponsored by NASA-MSFC (Contract 
NAS8-257 12). 
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FIGURE 2. TYPICAL LIQUID MOTI 




0.84«r/pr2 ^ - surface tension, p = density, r D = tank radius) which, for a 
lox tank of say 250 cm diameter, is less than 10 ^ cm/sec^. 

During the periods when an or biter -tanker is thrusting and in the 
process of changing orbits and speeds, both the propellants and the pay- 
load can slosh, similarly as experienced during the vehicle boost phan 
and thereby impose significant disturbing forces and torques on the orbiter. 

Flyback of Orbiter or Booster 

During the normal "flyback" of an orbiter or booster, residual pro- 
pellants in the amount of five to ten percent of capacity may remain in the 
tanks even after venting. Although this a relatively small mass, it can 
cause large disturbances because of the way in which the liquids are oriented 
during both the high angle-of-attack return and maneuvering just prior to 
landing. In both cases, the effective gravity is directed more or less 
perpendicularly to the sides of the tanks rather than along the tank axis, 
and hence the liquids, during flyback, are contained in what are effectively 
long, shallow tanks, as shown in Figure 3. Tests have shown that liquids 
in partially filled, long shallow tanks are very susceptible to small excita- 
tions and response in violent, large amplitude nonlinear motions (7]. The 

response can take the form of normal sloshing but, more commonly, 
traveling waves occur which reflect back and forth along the tank, impacting 
on the bulkheads each time; these possible reponses are shown in Figure 4. 
Such liquid responses can impose severe disturbances on vehicle guidance 
and control. 
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In some designs, water impact landing is contemplated. Large load- 
ings on the tank structure due to the gross motions of the rapidly decelerated 
propellants will occur during such landings, and may be of equal importance 
in the event of any "hard" landing. 

Aborted Flybacks 

If a mission is aborted, it is likely that at least a brief period of 
maneuvering and a? rplane-like flying will be necessary before the propellant 
tanks can be dumped or vented. Because of the large mass of liquid still 
onboard during this period, various propellant motions could be extremely 
deleterious to vehicle stability. 

If the main engines are shut down while the vehicle is flying thrvugh 
the sensible atmosphere, the liquids may impact on the upper bulkheads 
of the tanks because of the sudden deceleration resulting from air drag. 

This "dome impact" c.an rupture the bulkheads, and a fireball may result 
if hyper golic propellants are employed [8,9]. 

During the flying phase of abort, the large shifts in liquid eg and the 
fluctuating forces and torques caused by propellant motions such as shown 
in Figures 3 and 4 would be difficult to compensate for by the pilots. Obviously 
then, the possible propellant motions occurring during an aborted flyback 
need to be carefully investigated, and preventive measures provided for In 
the design of the space shuttle. 

CONTROL OF LIQUID PROPELLANT MOTIONS 
Because of the 100 re-use cycle requirement, the flyback capability, 
and the need to maximize the payload, the space shuttle is a weight-critical 


vehicle. Therefore, the devices used to control or dampen the propellant 
motions must be kept as light as possible. Flexible baffles, which effectively 
control sloshing and other forms of motion, are attractive in this respect. 

Although they have been shown to be feasible in tests with model tanks 
and conventional liq- .ds [10-12], flexible baffles need to be investigated 
with regard to (1) fatigue life, '2) compatibility with cryogenic liquids, 

(3) fastening techniques, (4) manufacturing techniques, and (5) in-place 
cleaning. c wRI, jnder contract to NASA-Langley Research Center#, is 
conducting such an engineering design analysis for space shuttle applications, 
which will attempt to answer these questions. This study is oriented toward 
eventually conducting a large-scale, proof-of-concept test program. 

Other specific measures to control liquid motions will undoubtedly 
rely upon adaptations of baffle and deflector configurations successfully 
employed in present vehicles [5]. Tank compartmentation may be par- 
ticularly effective in shifting liquid frequencies and in reducing the magnitudes 
of sloshing masses. 

INTERACTION OF LIQUID PROPELLANTS 
WITH ELASTIC STRUCTURE 

Experience to date has indicated that the overall dynamic response 
of essentially axisymmetric launch vehicle structures can be separated 
into basically (1) longitudinal accordion-like motions (pogo) and (2) lateral 
bending motions. Botli kinds of response are highly influenced by the presence 


♦Contract NAS1-1G074. 


of liquid propellants in the system; however, because of symmetry, these 
responses have, in the past, been essentially decoupled from one another 
(except for possible localized motions). On the other hand, the radically 
different and nonaxisymmetric designs for shuttle vehicles will no longer 
permit this simplification. Therefore, studies have already been initiated 
with the objective of making a relatively simplified parametric analysis of 
coupling between lateral and longitudinal structural motions for a typical 
shuttle design and, further, to determine the influence of the liquid pro- 
pellant dynamics on the coupled response. 

One study is being performed in-house at NASA-Langley Research 
Center and is described in detail elsewhere. Its objective, however, is to 
investigate the effects of joint stiffness between the booster and orbiter 
vehicles on the coupled vehicle vibration response, emphasizing rigid 
body and lateral bending modes and considering only solid mass representa- 
tions of the contained propellants. 

Complementary efforts are in progress at SwRI* to determine the 
influence of the liquid propellants on the coupled system of elastic shells. 

For this latter study, a sketch of the design of a simplified shuttle vehicle 
model is shown in Figure 5 and the model booster is shown in Figure 6. 

The two-tank booster is joined to a smaller two-tank orbiter by means of a 
common spar, and intermediate skirts and various rigid masses are utilized 
to simulate various components. This arrangement is an extension of some 

♦Contract NASI -9890. 
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recent research on purely axisyrnmetric longitudinal dynamics of a simplified 
launch vehicle dynamic model sponsored by NASA-Marshall Space Flight 
Center (13). Since the principal objective of this study is to determine the 
influence of liquid propellants on coupled longitudinal-lateral bending modes, 
natural frequencies will be measured on the model suspended in a simulated 
free-free condition (Fig. 6) for a variety of parameters such as propellant 
depths, ullage pressures, and orbiter attachment spring rates and position. 
Initial tests will be conducted with the tanks containing a simulated frozen 
liquid (using materials such as a mixture of sand and flou»-, or plastic 
beads); subsequent tests will be performed with liquids having the same 
density. Forced responses will also be determined by excitation at various 
simulated gimbal angles, required by the antisymmetric design. The gimbal 
angle changes with the va rying lateral position of the center of gravity as a 
result of propellant depletion during the boost phase. Also, both longitudinal 
and lateral components are required. It is anticipated that two preferential 
axes of motion will exist; one in the plane of, and the other normal to the 
plane of the drawing in Figure 5. A comparison of results using both 
simulated frozen and free liquids will permit an overall determination 
of the significance of the liquid coupling as well as provide some perspective 
regarding the increased complexity which will be required for pogo analyses 
of shuttle vehicle designs. To aid in these analyses, a complementary ana- 
lytical model is also being developed consisting of a series of linear and 
torsional springs, flexible beams, and rigid masses. Fach mass is to be 
capable of two degrees of translational and one degree of rotational motion, 




except for those representing the lateral sloshing and longitudinal liquid 
contributions. Mechanical models for lateral sloshing are those previously 
given in [2], while those for longitudinal liquid motion were recently given 
in 1 141. Initial plana are to study motion in the plane of the model only. 

CONCLUDING REMARK 

We have pointed out in the foregoing discussion some of the areas 
in which developmental efforts are required in order to insure that problems 
of propellant dynamics in shuttle vehicles can be satisfactorily analyzed 
and solved. Many of these are already being studied, either to delineate 
them in more detail or to effect their solution. Some of the more significant 
problems include: (1) propellant slosh characteristics in nonsymmetric tank 

configurations; (2) propellant impact on tank bulkheads, either at high-g as 
a consequence of abort, hard landing, etc., or at low-g as a consequence of 
docking, maneuvering, etc.; (3) flexible baffle technology; and (4) influence 
of propellant dynamics on elastic shell and vehicle vibrations. The coupling 
of the latter with other elements of the system, thereby leading to overall 
vehicle instability (pogo), is ofprimary importance in shuttle vehicle design. 
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TRANSIENT LOADING CONSIDERATIONS.' LOR SHUTTLE DESIGN 


E. F. Baird 


Grumman Aircraft 


Introduction 

This paper will review some of the transient-response problems of the space 
shuttle vehicle, and also introduce a combined test and analytic technique that could 
be a valuable supplement to the standard methods for calculating transient response. 
The technique relies on using experimentally determined impulse transfer functions to 
predict structural responses to longer -duration inputs. Before going into this tech- 
nique, the transient loadings the shuttle will encounter from liftoff to landing will be 
summarized. 


Fig. 1 lists the more important transient loadings due to launch, separation, 
docking, and landing. Since other papers will discuss so/nr. of these conditions in 
detail, such as launch and separation, comments on these items will be limited . 
Pyrotechnics have been omitted since these will probably not be important to over- 
all structural response. 
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Fig. 2 lists recommendations for limiting transient response during ignition and 
release at liftoff. Reference 1 will cover this subject in some detail, so only some of 
the more important considerations will be mentioned. These include flexibility of the 
launch pad structure, including the operation of the release mechanism, control sys- 
tem interactions, and lateral excitation due to wind effects and thrust build-up differ- 
ences among the various engines. 

Since these are principally ground conditions, it has been suggested that research 
in this area be directed primarily toward the design of ground-load alleviation devices. 
These studies, including dynamic-model and wind-tunnel tests, should focus on ? *ch 
items as shielding systems to protect against wind loads, and holddown and release 
systems which would limit the response of the vehicle elastic modes during engine 
start and at liftoff. 
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Booster engine shutdown transients should be mentioned briefly. On the Apollo 
Saturn V, the largest lateral accelerations on the LM during the launch phase occur 
after engine shutdown as a result of the longitudinal-lateral coupling (X-Z) between 
booster and spacecraft. 

In the shuttle, the relatively large eg offset between booster and or biter is ex- 
pected to result in fundamental modes which exhibit considerable coupling. A rela- 
tively refined clastic analysis and model program is, therefore, recommended to 
evaluate this transient condition. The large sizes of the matrices required can maki 
this analysis expensive, so improved matrix reduction techniques will become highly 
desirable . In addition, analytical and experimental studies should be performed to 
determine an optimum interstage arrangement and stiffness to minimize the adveise 
transient loads due toX-Z coupling. 
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The next problem to be reviewed, separating the orbiter and booster, 
the most critical in the shuttle mission. 


is one of 


One concept for accomplishing separation is shown in Fig. 4. The prime objec- 
tives of the scheme arc to provide reliable operation for both normal and abort situa- 
tions; to avoid, where possible, the harmful plume impingement of the rocket exhaust 
of one vehicle on the other; and to give positive assurance of no recontact of the ve- 
hicles under widely varying aerodynamic conditions. 

The dynamics of separation involve a number of complex problems: transients 
due to mechanism operation, engine shutdown and startup, and aerodynamic interac- 
tions of the two vehicles in proximity to each other. A combination of wind-tunnel 
testing and analytic backup is recommended to gain a fundamental understanding of 
the aerodynamic problems. Reference 2 will cover this topic in more detail. 
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Another problem In the shuttle mission is docking, as shown in Fig. 5. While 
docking-impact loads will probably not design major portions of the orbiter structure, 
they will be important in developing docking mechanism hardware and designing local 
structure, and will influence payload handling and orbiter control system design. 
Propellant behavior, as affected by Peking loads, is another important consideration 
that will require study. 

Two basic schemes for docking the orbiter to a space station are shown, togeth- 
er with some of the advantages of each. 

There are many other schemes for docking and payload deployment that range 
from using the payload structure as the docldng interface to designing for free-flying 
payload modules. Tradeoff studies involving the interaction of control-system dynam- 
ics, docking impact, and post-latching loads must be performed to establish a practi- 
cal system. Manned simulations of the approach phase of docking should be consider- 
ed, as well as model testing again backed by an analytic program. 
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The final step in the shuttle miscion is landing. In the early stages, shuttle 
gear design will probably follow transport-aircraft practices. Studies of the weight 
penalties associated with landing gear design criteria, as illustrated in Fig. C, should 
be performed to learn where significant savings can be made. The sensitivity of strut 
weight to sink speed is shown here. Since the landing gear itself is ot sizable weight 
and since landing imposes design-level loads on local structure and is of importance 
to overall vehicle response, advanced landing systems should be investigated. These 
should range from replacing wheels with brushes or skids to studies of air-cushion- 
type landing schemes. 
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In addition to analyses, conventional aerospace practice makes extensive use of 
testing. This would include mode surveys, vihraticn testing, and dynamic simulations. 
Modal information obtained from test is still difficult to use to predict internal load 
distributions, stress levels, or component accelerations. Dynamic simulations may 
verify the integrity of the structure, but may not yield basic information which can be 
applied to conditions that were not tested or which arise after the testing is completed. 
To overcome these limitations, Grumman has been developing, under MSC sponsorship, 
a combined analytical and testing technique that is called "impulse testing." This is 
expected to become a valuable supplement to the standard tests and analyses planned 

for the shuttle. 

The technique Is based on experimentally determining the responses, such as 
strains, forces, or accelerations, to short-duration pulses that approximate an im- 
pulse. These responses or transfer functions may then be used to predict the res- 
ponses to long -duration loadings. Fig. 7 and 8 give a quick review of the principles 
involved in the method. 

Suppose we are concerned with determining an internal response at R due to a 
transient loading at I. The loading of interest could arise from landing or separation 
dynamics or docking impact. The structure, which may be a portion of the SSV hard- 
ware or a detailed model of the entire vehicle, is mounted so as to simulate the de- 
sired support condition, and a short-duration pulse, in the desired direction, is ap- 
plied at I. As in modal analysis, the structure is assumed linear, and thus the res- 
ponse at R is directly proportional to the area under the impulse. When the response 
is normalized to a unit impulse, an "impulse transfer function" is obtained. 
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The next step in the procedure, shown in Fig. 8, consists of approximating a 
given time history of load application by a series of impulses. The response to this 
long-duration forcing function is then determined by performing a Duhamel integration. 

While the approach is straightforward, the problems of experimentally generating 
an approximation to an impulse, accurately measuring the area under it. and measur- 
ing the resulting low-level responses required a demonstration test to prove the feasi- 
bility of the method. 
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A typical test setup with a LM structural test article is shown in Fig. 9. The 
LM is supported from the apex of its outrigger trusses through a set of air springs to 
simulate a "free-free" condition. An electrodynamic shaker, programmed using an 
exact-waveform synthesizer, was used to apply impulses as large as 10,000 lb with a 
time base of about 1 millisecond to one landing gear attachment point. A force trans- 
ducer mounted to the shaker head was used to measure the impulse time history. In- 
strumentation to measure both strain and acceleration response was located at various 
points on the vehicle. A few gages were located on the outrigger truss members them- 
selves, very close to the impact point, while others were placed about 15 ft away, on 
the ascent stage aft equipment bay . 
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A typical impulse and several i.ormalized response time histories are shown in 
Fig. 10. The actual responses, as well as the pulse itself, were recorded on magne- 
tic tape, digitized, and then normalized to obtain the impulse transfer functions. The 
degree of variability in the technique, and in the test article, was assessed by com- 
paring impulse transfer functions derived from various short-duration pulses. In 
spite of suspected nonlinearities in the vehicle, repeatability was good for both strain 
and accelerometer measurements over a wide range of input forces and for many 
transducer locations. 
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Application of the impulse transfer functions to obtain LM responses to a long- 
duration tiansient loading is illustrated in Fig. 11. It is seen that excellent correla- 
tion is obtained between the response calculated using the impulse transfer function and 
the actual measured response to the long-duration loading. The results shown arc for 
a response location that is tar irom the input force location, and involve a complex 
load path from impact point to response point. Similar agreement between experimen- 
tal and predicted responses has been shown at other locations, for both strain-gage 
and accelerometer measurements. What has been demonstrated then is a technique 
which consists of measuring and cataloging structural dynamic characteristics in the 
form of impulse transfer functions, storing these intermediate results in a memory 
bank, and then recalling them to predict responses to arbitrary transient loads. Ap- 
plying this technique to the shuttle or detailed replica models will permit the predic- 
tion of loads due to a variety of transients. 
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Summary 


Prompt advances in several areas o£ technology will be required to handle the 
transient loading problems that will be encountered by the shuttle. 

The large- amount oi complex structure which has to be analytically represented 
to obtain reliable re sults will require extensions of current computing capability. 
Therefore, methods such as modal synthesis, which will relieve the size of the re- 
quired programs associated with comprehensive idealizations, will have to be utilized 

more extensively. 

Studies of launch holddown devices and ae rodynamic model studies of ground- 
wind alleviation systems should be pursued to reduce shuttle weight penalties for 
eround-loading conditions. 


Analytical and experimental studies to determine an optimum interstage arrange- 
ment and stiffness should be performed to minimize the adverse transient leads due to 
X-Z coupling, and the effects of aerodynamic interactions during separation. 

Since the size and complexity of the vehicle will limit the type of testing which 
can be performed to verify structural integrity, an impulse -transfer function technique 
should be considered for determining the loads and accelerations associated with many 
of the transient conditions which can not be tested directly. 


Sound judgment will be required to decide the depth of analysis, the extent of 
model testing, and the necessity for full-scale testing. 

Tho challenge, even when one connldero only the transient-loading problems. Is 
inde ed for m idal >le . 
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SPACE SHUTTLE LIFTOFF DYNAMICS 
Donald C. Wade 


NASA Manned Spacecraft Center 
Houston, Texas 



SUMMARY 


Because of the potential impact on structural weights and on crew and 
payload environments, space shuttle liftoff dynamics constitute an im- 
portant design consideration. The analytical simulation of space shuttle 
liftoff dynamics poses both an interesting and complex problem. This 
presentation discusses many factors associated with liftoff dynamics and 
the p'.tential areas affected by then. The merits of discrete versus 
statistical analyses of the condition are discussed and recommendations 
are made to assist in finding an optimum blend of criteria, analysis, and 
operations which will result in an efficient and adequate design. 


INTRODUCTION 


Probably no other space vehicle operation Is characterized by a busier 
schedule of electrical, propulsion, and mechanical systems operations 
that excite the dynamics of the structure than liftoff. For example, the 
Saturn V ignition sequence begins 9 seconds prior to release. During this 
9-second period, ignition signals are sent to engines in exacting order, 
valves open and close, engine thrusts build up, chamber pressures and 
scores of other parameters are sensed for an internal go-rio go decision, 
and holddown and access arms are retracted- to free. the. vehicle for liftoff. 
Variations in these parameters, coupled with a loading environment compli- 
cated by wind, gusts, vortex shedding, and engine gimbaling, result in a 
very complex structural design loading condition. A realistic structural 
simulation of shuttle liftoff is important because the resulting loads 
can affect the design of portions of the structure, as well as sensitive 
systems and payloads. An accurate simulation of liftoff is also important 
in establishing operational wind restrictions and for studying the effects 
on the crew and passengers. The unique configuration of a piggyback, 
winged shuttle orbiter configuration with a winged twelve-engine booster 
adds more complexity to the problem. 
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STRUCTURAL LOADING AT LIFTOFF 
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winds. An additional engine motion is induced when the servo-driven engines responding to these 
conmands experience an oscillation often described as "engine dither." 


FACTORS AFFECTING LIFTOFF LOADING 
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SPACE SHUTTLE STAGING DYNAMICS 


R. H. Schuett, M. C. Clark, and M. J. Hurley 


General Dynamics Convair Division 

Stage separation is recognized as a major Space Shuttle problem area. The 
piggyback arrangement precludes utilization of separation techniques developed for 
tandem vehicle stages. Also, since most shuttle configuration 3 are not symmetrical, 
thereby complicating interactions, experience gained from Til an III C solid motor 
separation and Atlas -Centaur fairing jettison is not directly applicable. Unlike present 
day launch vehicle stage separation, the depleted space shuttle booster is as massive 
as the orbiter element so that large intervehlcular interaction is probable. Abort 
separation is likely to yield the most severe separation condition since aerodynamic 
loading is significantly higher during the abort regime. Aerodynamics, including in- 
terference effects, will dominate the separation dynamics for all but the lowest dynam- 
ic pressures. 

The separating booster must clear or withstand the orbiter exhaust plume, as 
well as clear its hard structure under all probabalistic conditions. Following separ- 
ation, the booster element must initiate a capture mode and reorient to a trimmed 
attitude for subsequent reentry. The orbiter continues on its mission, or, for abort, 
prepares for reentry. 

This paper summarizes separation analyses conducted during the NASA Phase A 
Space Shuttle Study, the Air Force (SAMSO) Space Transportation System Study, and 
company sponsored research which evaluated separating multiple bodies placed sym- 
metrically about a center core. An aft hinge system with separation force supplied by 
inertial and aerodynamic reactions was recommended. With the piggyback arrangement, 
other schemes could be required. Concepts using rails, links, thrusters and pistons 
are discussed with advantages and shortcomings of each system identified. Factors to 
be evaluated for system comparison include system weight, performance losses, 
separation clearance time history, abort capability, and costs. 

Supporting research and technology necessary for selection and definition of a 
space shuttle separation system are identified. Covered in brief discussions are the 
three topics, interference aerodynamics, plume impingement heating and pressures, 
and propellant motions during abort. 


PROBLEMS 

UNSYMMETRIC, PIGGYBACK ARRANGEMENT 
CONTROL & RECOVERY REQUIREMENT 
SEPARATION OF BODIES OF EQUAL MASS 
SEPARATION IN HIGH DYNAMIC PRESSURE 
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SEPARATION DYNAMICS 

Separation is a combination of many transient events over a relatively short time 
period and the implication of one event on succeeding events must be considered in any 
analysis. The sequence is begun based on a pre -established staging velocity or time 
or could be triggered by a sensor which indicates some measure of propellant depletion. 
The first transient is then booster engine thrust decay. Associated with thrust decay is 
a relaxation of the structure responding to the changing acceleration vector. The com- 
bination of structural response and propulsion system characteristics have, on existing 
systems, resulted in thrust oscillations during engine decay. Gimbaling of the booster 
engines during this time period may be required for control in event oi unsymmetric 
engine shutdown or to trim the vehicle under additional forces such as orbiter thrust or 
separation system forces if those are activated prior to booster engine shutdown. Con- 
trol forces after engine shutdown are provided by the ACPS. 

The second series of events is associated with activation of the separation mechan- 
ism. Releasing the structural tie between stages by means of retraction pins, explos- 
ive bolts, or whatever, introduces two transients; the first is the shock associated with 
the mechanism and the second is the redistribution of loads due to breaking the tie point. 
These conditions are followed by activating the separation forces such as thrusters, 
pistons, inertial and aerodynamics forces, and the kinematics of linkages, guide ra Is, 
or hinges as may be necessary to control initial movements 

The next major seriea-of events are associated with orbiter engine ignition. Orbit- 
er engine ignition could be commanded at a pre-selected time after separation initiation. 
The rigid body and transient motions associated with engine thrust and activation of the 
orbiter control system are significant in separation system evaluation. Orbiter exhaust 
plume impingement causes a hostile environment to be imposed on the booster and could 
be a prime factor in establishing separation required at **rae of orbiter ignition, or es- 
tablish design requirements for portions of the booo. structure. 

Present during the separation sequence will be the aerodynamic forces. The inter- 
ference effects must be determined and included in simulations of separation dynamics. 
Aerodynamics will be of major importance for abort separation and of minor importance 
during normal separation. 

The separation dynamics through the preceeding events must consider the separ- 
ation of energy and momentum of one body into two bodies and the shift in points of body 
rotations after separation. For abort conditions, the phenomena of propellant motions 
in the booster tanks will add to the complexity of the analyses. 

Separation dynamics can be divided into two categories - rigid body motions and 
elastic response. The discussions in this paper will concentrate on the rigid body mo- 
tions and required timing sequences and forces necessary to provide an acceptable, 
safe separation. Elastic response should not be neglected in overall evaluation and 
selection of a separation concept. 
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SEPARATION DYNAMICS 
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PHASE A & STS STUDY RESULTS 

During the Phase A and STS studies Convair evaluated several separation systems 
for separating multiple bodies located symmetrically about the central or core vehicle. 

An aft hinge separation of the outer bodies from the core element while the core is under 
thrust was selected as being the most promising scheme for symmetric vehicles. Other 
concepts evaluated included a translational -rotational linkage system, an inclined ramp, 
and a system of lateral thrusters. The bases for selection of the aft hinge were; exper- 
ience on Titan III C solid motors and Atlas Centaur fairings; the system is passive; read- 
ily extendable to abort; and maximum energy was Imparted to the payload. Aft hinge sys- 
tem separation was initiated by releasing forward attach points. The boosters then are 
free to rotate about the aft hinge and will do so under the combined aerodynamic and in- 
ertial reaction loads. To determine aerodynamic interference effects, captive trajectory 
wind tunnel testing was conducted. 

The captive trajectory system permits the trajectory of a body separating from an- 
other (stationary) body to be determined experimentally in a wind tunnel. An analog com- 
puter uses force and moment data from a balance within the separating body to compute 
the resultant trajectory. Trajectory motion is simulated by the six-degree -of-freedom 
support shown. 

The simulation includes the aerodynamic characteristics of the separating body dur- 
ing and just after separation, as influenced by the flow field generated by the stationary 
body, together with the mass properties and propulsion characteristics of the body itseli. 

In this way, the effects of body release position and attitude can be studied. 

The model support, an electromechanical positioning system with all axes of moHon 
contained within a single mechanism, is independent of the stationary body. This mechan- 
ism has an envelope of motion lying within a cube about 30 inches on a side. Drive motors, 
located in a case below the tunnel floor, are printed-armature electric motors with ex- 
tremely fast response characteristics. 

After the free-stream flow has been established, the analog computer computes the 
trajectory using model strain-gage balance data in conjunction with body mass, moment 
of inertia, rocket thrust, altitude, and other data. This trajectory, transformed into 
velocity components, is then supplied to the support control drive motors, thus position- 
ing the separating model in a smooth, accurate simulation of the separation trajectory. 

The analog program is time-scaled; thus, what is seen is an apparent slow-motion move- 
ment of the separating body through its separation trajectory. The balance outputs are 
converted to aerodynamic coefficients, and the actual position and angle outputs are pro- 
cessed into full-scale parameters. 
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PHASE A & STS STUDY RESULTS (CONT’D) 

In conjunction with the analysis of space -shuttle separation, a digital computer 
simulation was developed to analyze the separation dynamics of hinged or linked 
lifting-entry vehicle clusters. The program computes the kinematics of separation 
of as many as five auxiliary (booster) elements attached around the periphery of a 
central core (orbiter) element in six rigid-body degrees of freedom for the core 
plus three for each auxiliary element. After separation, all elements are computed 
in six-degrees-of-freedom. 

Aerodynamic forces provided by the captive trajectory testing were used in the 
computer simulation. A typical abort separation display is shown on the opposite 
page. The correlation with the captive trajectory results was acceptable. 

Most present space shuttle configurations are not suitable for aft hinge separa- 
tion because the aft hinge system does require that the aft ends of all stages be 
reasonably flush. Also, the hinge reaction on the core, which.for symmetrical veh- 
icles was canceled by the body on the opposite side, could cause vehicle control 
problems. Therefore, this staging system evaluation must re-examine all concepts 
in view of the present shuttle configurations. 
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CONCKPTS 
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^ links free lo velooity. 

ative to the orbiter and the link rotati n - orbiter, the booster thrust 

By the time the booster plume comes n ^ on ^ tlme orbiter engines are 

:=r~ 

r, ■— - *nt time before Water eneine 

shutdown. 

BMiffllgl- Staging e^Z^ o" 

propellant depletion « to reduce veMole 

11 also'provlding reiative lateral velocity. Plutne 
pressures will impart lateral velocity to the booster. 

RAMP/RAILS - This sequence is very similar to the preceding system with initial 
separation travel being guided by ramps or rails in place of linkag . 

PISTON - Staging is initiated by an accelerometer which indicates zero or very small 

= 2 ^ 

both ends of the system. 

SOCKET - Very similar to the preview wncept. £££"££ 
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PRELIMINARY CONCEPT EVALUATION 

r inks-ROOST'ER THRUST - This system has some potential for normal separation, 

the aerodynamic Torres during the — t * 

cepts in payload capability . 

, nnmTFR THRUST - Primary advantages of this system are no coast time and 

1^”™^ abortoondltlons. Performance Is mmtlmlse - -I* the 
elimination of coast time. The improved probability of satisfying [abort t con 1 
is due to having the orbiter under thrust immediately, reducing the tim 
is in the vicinity of the booster which in turn reduces the impact of aer^yna 
force i. The support links are in tension during normal ascent flight tending to re 
duce weights. SUL impingement will be significant, but preliminary calculations 
have indicated pressures and temperatures are not intolerable. Primary l 
a^onhe system is the timing and controllability required between the various 
engine shutdown and startup characteristics. Also, vehicle control ^ 

ignition and shutdown may be delicate. The separation dynamics a„d ^er/or 
controllability are quite sensitive to booster thrust decay mid l orbiter ^hrust nse 
time histories. The time phasing between booster engine decay and orbiter engi 
ignition will also affect controllability and separation dynamics. 

RAMP/RAILS - No particular advantage can be cited for this concept over the pre- 
^TsT^T. Ramps/rails have more simplicity than linkages but also have more 
weight. Lateral separation would be less, and therefore the plume impingement 
effects would be more important. 
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PRELIMINARY CONCEPT EVALUATION 
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IMPINGEMENT SEN 


PISTON - This concept appears to be the most straight forward, simple approach. 
Sensors and timing requirements are very practical. Angular rates imparted to the 
orbiter can be minimized by placing the piston at the orbtter e.g. Assuming booster 
piston location between the LO and LH tanks, the angular rates imparted to the 
booster are not excessive. Negative aspects are: the relatively long coast time of 5 
to 8 seconds reduces performance, large force will be required to separate the veh- 
icles under abort conditions with additional complications due to booster/orbiter 
differential drag forces at ab^rt which could require guides of some type to react 
moments and axial forces; redundancy for fail safe, or fail operational, fail safe re- 
quirements will increase the system complexity and weight. 

ROCKET - Most of the comments of the piston concept apply to this system; however, 
the system does have more complexity and has the added problem of rocket plume 
impinging on the orbiter. The concept could use solid rockets or engines fed from the 

LO and LH tanks. Replacement or refurbishment would be necessary with solid 

rockets. An additional LO tank and considerable LOg plumbing would be required 
with the second option. A preliminary reliability assessment utilizing solid rockets 
indicated a triple redundant system would be necessary. The abort condition would 
require much larger thrust than normal separation. 
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PRELIMINARY CONCEPT EVALUATION (CON'T) 



TYPICAL SEPARATION TIME-HISTORY 

Lateral and axial separation time histories are shown for the Link-Orbiter Thrust 
and Piston concepts. Time zero is taken to be the time at which restraining members 
are disconnected. For the link concept, this means that the orbiter is up to full thrust 
at time zero, when the booster engines have decayed to about 10% of full thrust and the 
system is released. The lateral separation velocity is due to the rotating links and the 
impingement pressures on the booster. Axial separation is, of course due to T/W of 
the orbiter. 

For the piston concept, to match the lateral separation of the link system will re- 
quire a 200,000 pound force applied for 1 second. A 100,000 pound force applied for 1 
second will provide half the lateral separation distance, A piston stroke of approxi- 
mately 10 feet would be required with the 200,000 pound force and 5 feet for the 100, 000 
pound force. The orbiter was assumed to be up to full thrust at 4.0 seconds after 
separation is initiated, or 3.0 seconds after the piston has reached the end of its stroke. 

The significant separation distance difference between the two concepts is the 
axial displacement. At 5 seconds the relative axial movement has been 600 feet for 
the link system and about 100 feet for the piston system. The greater separation dis- 
tance reduces hazardous conditions during abort in event of an explosion. 

The piston system, with its lateral separation at orbiter ignition reduces plume 
impingement effects. The timing of the sequence is easily established. 

The link system maximizes performance, has abort potential, but suffers from 
the hostile impingement environment and could be sensitive to engine shutdown and 
startup characteristics. 
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TYPICAL SEPARATION TIME HISTORY 
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TIME (SEC. FROM SEPARATION 
SYSTEM INITIATION) 




IMPINGEMENT PRESSURES 

Typ'cal orblte, 

;rzr,r:: t * -rsrzx— trirrrr- 
— *» « — zi:i:z:L My 

nozzles ££ and the Meet, number, the static pressure was determined 

downstream of the oblique shock which exists at the Intersection of the flow stream 
and the element. 

From the Shows'" ' T he ptane'oMntersectton of the nows from the 

**» — 

S‘,'S“ reamr £ Cstim of’J ezines, bu, probably less than . psls. ™. 
pressure would act over a strip possibly 2 feet wide. 

=S:5S=Si= 

pounds. 

at u»o 0 .a, v-t-iii were calculated for one orbiter engine. The 
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>* *— forward movement 
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SUPPORTING RESEARCH & TECHNOLOGY 
Several key areas have been identified as requiring definition before selection 

Sr 

“r^rrezr:r=^=r. 

om in the vicinitv of the vehicle, these must also be simulated m the tests. 

* freedom captt ve trajectory capability, at leas, 8 or 9 degrees of freedom 
iTbe devetoped for Dual system verification. Buffe, response testmg tvtll be 
necessary if separation is to occur in high dynamic pressure. 

2 Plume impingement - Several separation concepts will incur considerable plume 

' rrr e r:: " rrre^irZst'^n^p^atio,, 

E^rrr^r;. 

Analyses techniques, verified by the above tests, should be developed to p e 

^th! Iks which in turn, will cause the booster to undergo erratic behavior. 
These muttons, although probably too slow to affect suparatl»eleare^e^u 
nlaee exce^ive demands on the vehicle control system. Analyses techniques 
should be developed which describe the propellant motions under abort conditi . 


SUPPORTING RESEARCH & TECHNOLOGY REQUIRED 





SUMMARY 


Investigations should continue on the alternative staging methods such as separ- 
ation rockets, pistons, rails, linkages or aerodynamic separation. Each technique 
has advantages and disadvantages. Fcr example, pistons reduce plume heating effects. 
Rockets will require additional attach points and doors. Aerodynamic separation will 
require a higher staging dynamic pressure (i.e. , performance reduction) but perhaps 
simplifying separation system hardware. These techniques must be reevatuated con- 
sidering both nominal separation and abort separat. m. Abort, in addition to large 
aerodynamic forces, will introduce unknowns in engine controllability (for propulsion 
system failures) and many other additional factors to be considered in separation sys- 
tem requirements and design. 

Evaluation criteria must consider separation distance time history from a clear- 
ance as well as potential explor, >n hazard for the abort conditions. Weight penalties 
associated with each concept must be identified, including systems and structures re- 
quired in the concept and such factors as added thermal protection to account for plume 
effects. AV losses incurred during the coast period or by application of separation 
forces are to be considered. Disturbances to the orbiter and booster should not re- 
quire excessive additional vehicle control capability. The ability of the system to 
perform under abort conditions for the total ascent phase or portions thereof muut be 
established consistent with the abort philosophy . 

Research must be conducted to establish interference effects on aerodynamics 
during separation. Exhaust plume impingement analyses must be developed and 
verified by tests for adequate representation of separation. For abort conditions, 
analysis techniques describing motions of the unused propellant in the booster must 
be developed and verified by tests. 
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GROUND-WIND-LOAD CONSIDERATIONS 

FOR SPACE SHUTTLE VEHICLES 

by Wilner H. Reed HI 

NASA-Langley Research Center 
Hampton, Virginia 


i::;roductic'n 


Four years ago a technical meeting on ground-wind load problems in 
relation to launch vehicles was held at Langley Research Center (ref. 1). 

In the concluding session of that meeting, a panel addressed the question: 
"Where do we go from here? - the designer's viewpoint." We probably 
shouldn’t be too critical of that panel for failing to include in their 
crystal-ball-look into the future, consideration of vehicle launch opera- 
tions such as the following: Take an orbiter vehicle that resembles a 

Douglas DC-8 jet transport; strap it piggyback to a booster vehicle that 
resembles a Boeing 747 super jet transport; mount the combination vertically 
on a launch pad; design the system to withstand a 70-knot wind, and then 
be launched vertically in a 35-knot wind. Today, we call the concept 
Space Shuttle. 

In the area of ground-wind loads, as in other technical areas, the 
space shuttle concept presents a full measure of challenging problems 
For example, in addition to the usual ground-wind load problems dealt 
with on previous launch vehicles, such as steady drag loads, gust loads 
and vortex shedding, we also foresee new problems, some with interesting 
names such as "stop-sign flutter" and "galloping" stability. In this 
paper, we will attempt to identify and focus attention on those arc s 
which ' space vehicle designers have heretofore been able to ignore, '<ut 
which are likely to be of key importance on Space Shuttle configurations. 
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SPACE SHUTTLE GROUND- WIND-LOAD CONSIDERATIONS 


This figure presents a brief list of factors to be considered in 
designing a vehicle to withstand the effects of ground wind. First, one 
must specify an appropriate ground-wind environment for various phases . 
of operation. The design wind environment for space shuttle vehicles is 
specified in reference 2 in terms of steady wind profiles and turbulence 
spectra. A requirement stated in reference 2 is that the space shuttle 
shall be capable of withstanding peak ground winds up to 72 knots at the 
60-foot level during stay-periods on the launch pad. The grcund-wind 
environment creates loads on the vehicle which are both steady and dynamic 
The dynamic load components can be associated with fluctuations in the 
wind itself, vortices shed from the structure, or some form of aerodynamic 
instability such as "stop-sign flutter" or "galloping" stability. 

These loads and responses create 3uch problems as the following: 
stresses in the vehicle and the associated tie-down structures; clearance 
requirements to allow for relative deflections between vehicle and ground 
service structure; fatigue due to the accumulation of long wind-exposure 
times on reusable vehicles*; ground handling of winged vehicles during 
towing and erection operations; and finally, guidance alignment difficul- 
ties encountered on a swaying structure. 



*Note: If the average time on the launch pad were one week per launch, 

a shuttle vehicle would accumulate approximately two years of 
wind exposure time during its desigh life span of 100 missions. 
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PACE SHUTTLE GROUND-WIND-LOAD MODEL STUDIES 


This fioure depicts various ground -wind- load model studies v*i ich are 
planned for Langley wind-tunnel facilities in support of the space shuttle 
program. The approach is step-by-step, progressing from simple models in 
low-speed wind tunnels to the complete neroelastic modeling at high 
Reynolds numbers of the vehicle-launch pad complex. 

First simple rigid models with interchangeable lifting surfaces will 
be utilized to isolate and study some spec if ic aerodynamic instaoi. it es. 
Galloping instabilities - the tendency of noncircular blunt bodies to 
oscillate in the cross-stream direction - are being explored m the Langley 
12-Foot Low-Speed Wind Tunnel. More will be said about these studies later 

in this paper. 

In a second phase, with the same model used in galloping studies, 
stop-sign flutter - a stall-flutter type phenomenon which involves single- 
oscillations of s lifting surface -ill be investigated 
in the 7- by 10-Foot High-Speed Tunnel — wi 1 1 emp oy a 
special oscilla ting-sting mount developed for measuring airplane ro 
damping derivatives. 

The third phase of ground-wind-load model studies will be conducted 
in the Transonic Dynamics Tunnel in freon. In these studies, light -weigh 
rigid models of coupled booster-orb iter configurations will be 
to^a flexible mount system. The mount has variable stifi fness and damping 
features in order to simulate the fundamental cantilever bending and 
torsion mod. of th. vehicle. During this phase of the study. In <*ich 
Marshall Space Flight Center and Langley Research Center will be joint 
n^rticioants it is expected that considerable effort will be devoted to 
the evaluation of methods for alleviating excessive wind-induced oscillations 

of the shuttle vehicle. 

It will be recalled (ref. 3) that in similar ground-wind-load studies 
„„ Saturn V a requirement for additional structural damping was establish- 
ed anSwas subsequently implemented at the launch site in the form of 
islwr LJ ths vehicle end Lunch touer. After th. sp.c. 
stottl. vehicle design hes bsen defined, studies of complete 
models of the vehicles and launch structures will be conducted in th. 
Transonic Dynamics Tunnel. Included as an important aspect of these 
Si studies will be evaluations of wind-load effects Airing vehicle 

erection. 
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AERODYNAMIC FORCES AFFECTING GALLOPING STABILITY OF BLUFF BODIES 


Gall opinp instability is so-named because it is descriptive of the 
low-frequency;’ high-ompl itude oscillations that are sometimes observed 
on ice-loaded transmission conductor lines exposed to wind. For such 
instabilities, an essential feature is that the body cross-sec * 

is such that a steady-state lltt force is developed and that the rate of 
chance of this force with wind incidence angle (angle of attack) . 
negative. With reference to figure 3, it can be seen that when this con- 
dition exists, the change in lift force associated with the change 
relative wind direction due to body motion, y/v, acts in the ® ame d 
tion as the motion, and therefore tends to reinforce, rather than ampe , 
the motion. As a necessary condition for galloping, Den Hartog gave the 
following criterion in reference 4: a section is dynamically unstable 

if the negative slope of th4‘$orce is greater than the drag force, 


i.e. , 


&L + 

Ax 


D < 0 


Although the galloping phenomenon is well-known and has been studied 
extensively in relation to transmission lines, towers, and bridges (see, 
for example, references 4 through 8), it has received little attention 
by launch vehicle designers. This is not surprising since ® os t existing 
launch vehicles have basically circular cross-sections which cannot pro- 
duce the necessary steady-state lift forces required for galloping excita 
tion. On of the other hand, looking at the typical 

section of some space shuttle configurations, one might conclude that the 
prototype design had been Den Hartog's galloping demonstration model 
(ref 4 page 301). Thus, it is expected at the galloping phenomenon 
will play a prominent role in space shuttle grc :nd-wind-load studies. 


Footnote : 


In the oral presentation, motion pictures were shown which 
demonstrated galloping and stop-sign flutter on small wind 
tunnel models. Copies of the film are available for loan 
from NASA-Langley Research Center. 
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AERODYNAMIC FORCES AFFECTING GALLOPING STABILITY OF BLUFF BODIES 




SPACE SHUTTLE GROUND-WIND- luOAD MJDELS 


To explore the galloping characteristics of representative spate 
shuttle configurations, static aerodynamic forces and moments were measured 
on the models shown ir figure 4. These st lies, referred to earlier in 
figure 3, were conducted in the Langley 12-Foot Low-Speed Wind Tunnel. The 
model consisted of a basic body onto which various lifting surfaces could 
be attached. Fiberglas-foam sandwich construction resulted in a light- 
weight, rigid model structure. The model was mounted on a 5-componeat 
turntable balance located beneath a ground board, and could be rotated 
about its longitudinal axis so as to cover an 180-degree range of wind 
azimuth angles. In addition, two accelerometers were mounted near mid- 
body to sense horizontal accelerations of the model. 

Some pertinent model properties are as follows: The body cross- 

section had an 8-1/2 inch wide flat bottom; flat sides slope inward at 
7-1/2 degrees coming tangent to a 3- inch radius circle irtiich forms the 
upper side; the bottom corners have 3/8-inch radii. Overall dimensions 
of the body are 62- inch length, 8-1/2 inch width, and 12- inch depth. 

The wing of the straight-wing configurations had a 7.36 aspect ratio, 

0.42 taper ratio, 7 degree dihedral angle, and an NACA 0012-64 airfoil 
section. The wing of the delta-wing configuration had a 50 degree swept 
leading edge, 2.66 aspect ratio, 0.14 taper ratio, zero dihedral angle, 
and a 3-percent-thick symmetrical airfoil section. 
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STATIC AERODYNAMIC COEFFICIENTS FOR A SPACE SHUTTLE CONFIGURATION ON 

LAUNCH PAD 


This figure shows the variation of static lift and drag coefficients 
with wind azimuth angle for the straight wing configuration pictured in 
the previous figure. (The reference area is taken to be the total pro- 
jected area of the wings, body, and horizontal tail.) Regions of possible 
instability, as given by the Den Hartog criterion, are indicated by 
shaded sections on the C L curve. (Note that negative C L has been 

plotted.) From the slope of the lift-curve, it can be seen that galloping 
instability should be mpst pronounced when the flat bottom side of the 
body faces upstream ( a = 0) ; conversely, the system should be most stable 
when the flat side faces downstream (a = 1800). (These stability trends 
were clearly displayed in the motion picture demonstrating galloping 
stability. ) 




GALLOPING STABILITY AETODYNAMIC COEFFICIENTS FOR SOME SPACE S 'LITTLE 
CONFIGURATIONS ON LAUNCH PAD 


Thi3 figure shows lift and drag coefficients plotted against wind 
azimuth angle, as in the previous figure, for three conf igirations : 
fuselage alone, straight wing, and delta wing. The reference area in 
each case is the total pr jected area of the planform in question. 


154 





PREDICTED GALLOPING STABILITY BOUNDARIES 


In this figure the aerodynamic coefficients presented in figure 6 have 
been utilized to predict galloping stability boundaries for each of the 
three configurations. The ordinate in the figure is a nondimensional 
velocity ratio V/fD in which V is the wind velocity, D ia the body 
width, and f is the natural frequency (in Hz) of the vehicle for an 
assumed single-degree-of-f reedom translation mode in the cross-stream 
direction. The abscissa is a nondimensional inertia damping parameter 


2 6m 
s 



in which 6 is the log decrement damping of the structure, m, average 
s 

vehicle mass per unit length, and P, air density. These stability 
boundaries can be derived by equating the work per cycle done on the body 
by aerodynamic forces to the work per cycle dissipated by structural 
damping. The shaded area in the figure gives an approximate Indication 
of the region covered by representative cantilevered- mounted space shuttle 
boosters exposed to the design peak wind speed. (Structural damping was 
assumed to be one percent of critical damping.) The right side of the 
shaded area represents the full-fuel condition, and the left side the 
empty-fuel condition. Some key points to be made about this figure are: 
lifting st .-faces attached to the body have a strong stabilizing influence 
on galloping; at the design peak ground-wind-speed for the space shutt e 
(72 knots at 60-foot height), the figure indicates that a straight wing 
shuttle booster configuration will experience galloping for the empty 
fuel condition; and, for the body without wings, galloping would occur 
for all weight conditions. 





ACCELEROMETER TIME HISTORIES 


MEASURED ON SPACE SHU'ITLK (iROUND- WIND- 


LOAD M3 DELS 


Tlii-- * ipure shows sample accelo ration t ime histories measured on the 
under Identienl wind tunnel r.,ndltt..„n (V . 50 «/«*. 
a »”()). " The fundamental frequency of the model tn e aeh ease was app oxr 
m,i te l y 7.0 Ilz, and the cor respond in,- mode slu.pe was essent ui 1 1 y • that 
rlcid model pivoted near its base. It is gratifying 0 0 e ‘ ‘ ' £ 

results support the predictions shown in figure 7 regarding the 
the wing on galloping stability. 
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SUMMVRY AND CONCLUDING REMARKS 


Thf requirement for reusable spare shuttle vehicles to withstand tlu- 
effects of ground winds poses some formidable problems. This paper _ 
attempts to identify arxl focus attention on some new problem areas whi- l. 
designers of past launch vehicle systems have been able to ignore, but 
which arc very likely to be significant on space shuttle vehicles. Among 
such problems are: structural fatigue due to the high accumulated wind 

oxpoHum time, on a reusable vehicle; ground handling during erect ion of 
large winged vehicles exposed to wind forces; vortex shedding,, stop- sign 
flutter, galloping stability, and possibly other wind-induced aer. dynamic 
instabilities associated with the geometric shape of space shuttle con- 
f igurat ions . 

Wind-tunnel results obtained from a preliminary invest igation of 
galloping stability for models of typical space shuttle configurations 
are presented. Indications from this study are: wind-induced galloping 

oscillations could occur on shuttle vehicles within the design ground 
wind speed range; the light-weight vehicle condition is least stable; _ 
and lifting surfaces attached to the body tend to have a strong stabilizing 
influence on galloping stability. 
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WIND TUNNEL SIMULATION OF GROUND WIND 
WITH POSSIBLE APPLICATION TO SPACE 


SHEAR AND TURBULENCE SPECTRA 
SHUTTLE LAUNCH PROBLEMS 


R. J. Templin, National Research Council, Canada 


I NTRODUCTION 

It seems to be generally agreed that wind tunnel 
testing of aeroelautic models is an essential part of 
the investigation of the dynamic effects of ground-level 
winds on launch vehicle configurations. However, a 
deficiency in practically all such tests so far carried 
cut has been the failure to simulate, in the wind tunnel, 
the correct wind mean velocity shear, turbulence intensity 
and power spectrum. This deficiency becomes more serious 
when vehicle configurations depart from simple cylindrical 
shapes, as is the case with the space shuttle launch 
configuration, because of the increased difficulties in 
the analytical prediction of ground wind loads. 

Three techniques are briefly described, which have 
been developed or put to use in Canada, for the wind 
tunnel simulation of low-altitude, neutrally stable winds, 
mostly for non-aerospace structural applications. Each 
has its ■ advantages and disadvantages, but any of them may 
be applicable in future launch vehicle wind tunnel tests. 
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JUPITER MO DEL IN SIMULATED TERRAIN ROUGHNES S 

Probably the most reliable method of simulating 
the neutr ally-stable gound wind ayer in wind tunnels 
is that pioneered by Jensen in Denmark, and further developed 
in specially-designed wind tunnels at Colorado State 
University (Reference 1) and at the University of 
Western o.t.rlo (Uln*. 2). I« this technique, n.t.t.l 
terrain roughness is simulated along the floor o 
tunnel test section, which may be up to 90 fe^t in 
length. This permits the growth of a boundary layer 
between two and four feet thick near the downstream end. 
Measurements of mean velocity profile, turbulence 
intensity variation with height, and the power spectrum 
of turbulence within this layer are in good agreement 
with full-scale measurements over various types of terrain. 

As far as is known, only one launch vehicle 
configuration has been tested using this wind tunnel 
technique. Figure 1 is a photograph of a 130 scale 
aeroelastic .'odel of the Jupiter vehicle in the University 
of Western Ontario tunnel during 1966. This investigation 
was carried out under a NASA grant, as a part of a research 
program to compare dynamic model measurements with full- 
scale measurements of vehicle response. 

Although this method appears to be the most highly 
developed technique for ground wind simulation in wind 
tunnels, it is not adaptable to conventional aeronautical 
tunnels, since it , -quires a test section length about 
10 t imes the section width or height. 






RESPONSE OF AE ROELASTIC JUPITER MO DEL _IN TURBULENT 
BOUNDARY LAYER 


Figure 2 shows the unsteady response, in 
non-dimensional form of the 1:30 scale aeroelastic 
Jupiter model in the University of Western Ontario 
boundary layer wind tunnel. The curve 
shows the R.M.S. deflection at the top of the model 
in the wind direction and the curve labelled o is 
the transver 8 a response. H is the total heig o 
the model and the mean wind speed U H is that at 
height H. The reduced wind velocity has been made 
non-dimensional by dividing by the product of model 
natural bending frequency f and the diameter D. 

The figure illustrates one of the effects of 
simulating, as nearly as possible, the complete v n 
turbulence spectrum. In tests of similar models in 
smooth wind tunnel flow, the phenomenon of vortex 
shedding from tail cylindrical structures often results 
in a large-amplitude response over a narrow resonant 
speed range, with little response outside this range. 

In turbulent flow, especially with substantial low- 
frequency content, the vortex-shedding resonant response 
may be entirely absent or greatly reduced. The small 
oeaks in the curves of Figure 2 occur at a value of 
Jhe reduced velocity (which is a form of inverse Strouhal 
number) where vortex-shedding resonance would be expected. 
In this case, however, the amplitude is almost lost in 
the overall response which now extends over the whole 
velocity range. 
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c I TY OF MONTREAL MODEL IN N,R.C. 6 x 9 FOOT TUNNEL 

Figure 3 is a photograph of a 1:400 scale model 
representing about one mile of downtovm Montreal in 
the National Research/National Aeronautical Establishment 
6 x 9 ft. wind tunnel in Ottawa. The purpose of this 
test was to obtain spectra of unsteady pressures on 
the face of the tall building near the centre of the 
picture, for comparison with full-scale measurements 
that had been made previously. 

The photograph also shows an array of 4-foot 
plywood "spires" at the upstream end of the test 
section. This technique has been developed as a means 
of producing thick, turbulent shear layers in the 
relatively short test sections that are typical of 
aeronautical wind tunnels (Reference 3). At a downstream 
distance of about 3 or more spire heights, the boundary 
layer characteristics appear to settle down so that 
the mean velocity profile, turbulence intensity and 
spectrum are reasonably similar to corresponding full- 
scale data. The spires are tapered in order to control 
the mean velocity profile, and their coarse geometry is 
required in order to generate turbulence of suf t icientiy 
large scale. There is some similarity between this 
technique and one developed in the U. K. by Counihan 
at the Central Electricity Research Laboratories (Reference 4). 

Techniques of this type appear to have reached 
a stage of development such that they can be used with 
some confidence for aeroelastlc response measurements 
on tall structure^, A possible disadvantage is tha 
the drag of the spire array is large, and in some wind 
tunnels may cause fan stall. 










TALL BUILDING PRESSURE SPECTRA 

Using the spire technique, unsteady pressure 
soectra have been determined at several locations on 
the tall building model shown in the previous Figvr re. 

A representative result is shown for comparison with 
full-scale measurements in Figure 4. 

In this graph the spectra are normalised with 
reenact to the mean squared pressure fluctuation, 
and the reduced frequency has been made J° nal 

bv multiplying by the height at the top of the roof 
tower (800 ft. full scale) and dividing by the mean 
wind speed measured at that height. 

The shaded area representing full-scale 
measurements is an envelope drawn around »*“*** 
number of separate spectra, which have considerable 
scatter among themselves. The agreement between 
and full-scale measurements is reasonably good, 
thereare indications that the model spectrum is shifted 
toward higher frequencies. The comparison of the 
spectra of longitudinal turbulence measured at the 
height H (now shewn here) display. * the same trend, 
suggesting that the spire technique has not fully 
reproduced the largest scale motions in the lower 

a tmo sphere . 

The advantage of this technique is its adaptability 
to any low speed wind tunnel, with a potential for 
producing shSar layers of thickness approaching the 
tunnel height. A disadvantage in some cases is the la g 
drag of the spire array. 








UNIVERSITY OF TORONTO TURBULENT JET LAT TICE 

A method for generating shear flow with 
associated large-scale turbulence in wind tunnels 
has been under development at the University of 
Toronto for about two years (Reference 5 ). In this 
method a lattice of jet outlets oriented parallel 
to the tunnel axis produces the turbulent flow, and 
by separate regulation of the jet pressures can be 
made to control the mean velocity gradient. A 
photograph of the small-scale experimental jet array 
is shown in Figure 5 . In experiments carried out 
so far, the jet array also constitutes the sole 
method of propulsion of the wind tunnel, which is 
of the open circuit type. 

Although this technique has not ye* been 
applied in large-scale wind tunnels, it offers 
several possible advantages. It does not introduce 
large drag into the tunnel circuit, and in fact can 
be used for propulsive p»rposes. It would appear to 
be more flexible with respect to the control of 
mean velocity gradient. A possible disadvantage may 
be its cost and complexity when adapted to large 
scale, but if coupled to a short, large-scale open- 
circuit tunnel this complexity may be offset by the 
fact that the jet system can also act as the tunnel 
drive sy s tern . 







LONGITU DINAL VE LOCITY SP ECTRUM WITH, 1EU^-JNDUCK.D 
SHEA R FLOW 

When the array of jeta ahown In the previous 
Figure are blown with uniform momentum so ns to 
produce a uniform mean velocity distribution further 
downstream, the turbulence characteristics are 
somewhat similar to those downstream of a square 
lattice of drag elements, except that, for a given 
lattice mesh size, the longitudinal Integral scale 
of turbulence is about twice as large for the jets 
as for a drag lattice. The R. M. S. turbulence 
intensity decreases with increasing downs tr earn 
distance, and is about 10 percent of the mean velocity 
at about 15 mesh widths. 

The character of the turbulent flow changes 
in certain respects when the jet^ are controlled 
so as to produce an approximately uniform mean shear 
across the tunnel. The R. M. S. turbulence, as a 
percentage of mean velocity, now varies through the 
shear layer in a manner somewhat similar to that in 
a boundary layer, and on the tunnel centre-line it 
remains roughly constant with dowstream distance. 
Furthermore, the integral scale becomes comparable 
with the width of the shear layer itself, and is 
much larger than that produced In uniform flows. 
Figure 6 shov s one measurement of the longitudinal 
power spectrum of turbulence in shear flow. Its 
overall shape is in at least qualitative agreement 
with that in the natural wind layer. 
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CONCLUSIO NS 

This paper has briefly described three techniques 
for the wind tunnel simulation of the natural wind layer 
at low altitudes, with possible application to space 
shut tie* launch vihicle problems. All happen to have been 
recently developed or adopted in Canada, but this choice 
is due only to the author's familiarity with them, and 
related developments have taken place in other countries 
as well. They appear to offer means of rectifying a 
deficiency that has charac terieed almost all previous 
wind tunnel testa of aeroelastic launch vehicle models. 

In particular: 

1. Thick turbulent boundary layers that simulate the 
neutrally stable ground wind layer can be generated over 
a long fetch of surface roughness in specially-designed 
wind tunnels with long test sections. Tunnels of this 
type have been in operation for some years at' Coiorado 
State University and at the University of Western Ontario. 
The latter has been used for exploratory aeroelast c 
model tests of a Jupiter launch vehicle. 

2. At NRC/NAE in Ottawa, an array of shear layer "spires 
has been developed, which generates thick shoar layer# 

in conventional tunnels, with velocity profiles, 
turbulence intensity, scale, and spectral shape ain>il a * 
to the neutral wind layer. A similar technique has been 
developed in England. 

3. The University of Toronto Institute for Aerospace 
Studies has under development a technique for generating 
turbulent shear layers by means of an upstream array 

of controllable air jets. In an open circuit tunnel, 
the let array appears to be capable of providing the 
tunnel propulsion system, and of generating shear 
with a depth comparable to the tunnel height, and with 
appropriate turbulence characteristics. 
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LIFTING AND CONTROL SURFACE FLUTTER 


by Robert C. Goetz 

National Aeronautics and Space Administration 
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Hampton, Virginia mr ^ • 
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SUMMARY 

Employing the space shuttle vehicle as a ceuter of discussion, several 
potential problem areas of lifting and control surface flutter are 
surveyed. Topics very briefly discussed include flutter of planar 
lifting surfaces, stall flutter, interference flutter, control surface 
instabilities, and lifting surface - body flutter. Proposed or needed 
research in those areas where the state of the art is not presently 
sufficient to insure a safe and reliable 3pace shuttle mission are 
delineated. 


INTRODUCTION 


The degree of interaction between aerodynamic, elastic, and inertial 
forces - which are the basic ingredients causing lifting and control 
surface flutter - is very much a function of a vehicle's design ard 
the severity of its operating mission. In the concepts for a 
recoverable and reusable space shuttle system, one sees emerging a 
new generation of vehicles: the offspring of a marriage between the 

design technologies of airplanes and launch vehicles. The vehicle 
mission includes vertical launch, separation, orbit, winged reentry, 
and conventional landing. This broad mission, together with the 
criticality of the vehicle's payload to weight ratio, highlights many 
areas of flutter technology that need refinement, and as with any 
new class of vehicle, new potential problem areas. The space shuttle 
proeram will thus push the state of the art of a great many disciplines, 
including the field of aeroelasticity as pointed out by Garrick (ref. 1) 
in a recent international talk. 

In the present talk, several potential problem areas of lifting and 
control surface flutter will be discussed and indications given as 
to their relative importance to the space shuttle program. Proposed 
or needed research in those areas where the state of the art is not 
presently sufficient to insure a safe and reliable space shuttle 
mission are delineated. 
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POTENTIAL FLUTTER PROBLEM AREAS 


* ♦ «.! flutter problem areas and their relationship 

Some of the potential flutter pr ^ indicated in figure 1. 

to the space shuttle opera 6 dynamic pressure (max q), 

2Lf - SS- ... u— - - — 

of the figure. 
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to the current state of t e area will have to be refined to 

For the space shuttle mission, this a , at high angles of 

include blunt, thick, lilting but indicated by the flight 

attack over the complete speesp region and max q, reentry, 

events of ascent through the transonic no* a 
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INTERFERENCE FLUTTER 





PRECEDING PAGE BLANK NOT HIMED. 


The fourth area entitled control surfaces includes such known phenomena 
as transonic buzz (ref. 13) and phenomena referred to more recently 
in the literature (refs. lU and 15) as hypersonic buzz. This latter 
case is an area where the technology learned in more recent hypersonic 
lifting-body vehicle programs is still being assessed. Also included 
in this area i 3 the field of aercservoelasticity. This field concerns 
the interaction of aerodynamic, elastic, inertial, and automatic 
control forces, which can considerably influence the control surface 
limit loads, fatigue load spectra, flutter speed, and the ride/handling/ 
control qualities of the vehicle. 

The fifth and final area is entitled lifting surface - body flutter. 
While the individual stages taken separately do not appear to be outside 
the realm of the current state of the art, the combined configuration 
is a more formidable problem. Obvious difficulties will be encountered 
in wind-tunnel aeroelastic testing programs, which in the past have 
always been conducted for any new vehicle design to prove flutter 
clearance and aeroelastic integrity. From an analytical viewpoint 
the problem is even more difficult. It is sufficient to mention here 
that the analyt cal technology at the present time is not adequate . 

A more detailed discussion of the status of unsteady aerodynamics 
associated with this area will be given by Mykytow in a subsequent 
talk. 

For the remainder of my talk, I wish to briefly present some recent 
results obtained at Langley Research Center that are applicable to the 
space shuttle vehicle in three of these areas: (l) stall flutter, 

(2) interference flutter, and (3) control surfaces. 
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STALL FLUTTER STUDY OF SPACE SHUTTLE WING CONCEPTS 



The first results are from an exploratory experimental study of the 
stall flutter characteristics of two space shuttle wing concepts (ref. 0]. 
The two wall mounted semispan wing models are shown in figure a. ine 
wing configuration of the left of the figure was a 50 clipped delta. 

It had a symmetric airfoil with a constant 3-percent maximum thickness- 
to-chord ratio. The configuration shown on the right was a straight, 
aspect ratio 7 wing which had a 12-percent thick, 6U series airfoil 
shape. The investigation was conducted in the Langley transonic dynami 
tunnel at Mach numbers from about 0.2 to 1.2 and Reynolds number per 
foot from about 0.15 to 2.5 x 10 6 . general testing procedure was, 

with the model at zero angle of attack, flow was established at a given 
Mach number and dynamic pressure. The model was then rotated at a 
nominal rate of 0.6 degree per second up to an angle of attack 
90 degrees. This angle of attack sweep rate is approximately 10 times 
slower than full-scale rotational rate during transition. Strain 
(taces at the model root indicating the relative magnitude of the 
strains in the bending and torsion degrees of freedom were continuously 
recorded, thus allowing correlation between the model response, 
model altitude, and wind-tunnel conditions. 




STALL FLUTTER BOUNDARY FOR AR ■ 7 STRAIGHT WING 


A summary of some of the stall flutter results obtained from this 
investigation is presented in figure 3. The stall flutter boundary 
sE if for the aspect ratio 7 straight ving. (Stall flutter was not 
encountered on the clipped delta wing over the range of the parameters 
of the investigation •, this result being consistent with previous 
studies.) The stall flutter boundary is presented as a ot 

Mach number and angle of attack. A flutter velocity < 
is also shown on the right of the figure. This < e°ef fieient is a 
function of the flutter velocity, V, wing semichord at the 3/4-semi 
span, b, and the circular frequency associated with the wing torsion 

mode, . 

The boundary in figure 3 indicates that the stall flutter region for 
this wing model was confined to angles of attack less than about 20 
and to Mach numbers less than 0.85. For the lower Mach numbers .the 
boundary limit for the stall flutter region is seen to occur near a 
valued the flutter velocity coefficient equal to 1. Thus . no “J 
stall flutter regions were found beyond those which have been studied 
and discussed extensively in the literature. However, it should be 
noted that this stall flutter boundary is not general but ratter 
the boundary found for this particular wing model 

tigation- It has been found in previous studies (ref. 5) that stall 
flutter boundaries may be altered by configuration aerodynamics v 
well as by structural parameter changes. For example, stall boundari 
are fmctionsof Reynolds number as well as Mach number and flutter 
velocity coefficient. 

Shifting gears now, I would like to turn your attention to another 
flutter area: Interference Flutter. 
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ANGLE OF ATTACK, DEGREES 



BOUNDARY FOR AR - 7 STRAIGHT 





BIPLANE FLUTTER PROGRAM 


An exploratory program was conducted In the Langley 26-inch transonic 
uX Jl to letemln. tb. n««r oh.r«t«rl.tlc. of » ^pect 
ratio 8, unBwept wing model in proximity to a go omecrically 
model as shown in figure 1*. Both models had^-peroontthick,^ aeries 
airfoil shapes. Flutter data have been obtained with the models leading 
edges inline and the distance between them varied from about 1/** to 
2-1 A mean aerodynamic chord lengths over the Mach number ra **8® rom 
0.1* to 1.3. Additionally, the coplanar separation distance be tween 
the models was also varied fore and aft one mean aerodynamic chord. 






FLUTTER CHARACTERISTICS OF WINGS IN CLOSE PROXIMITY 


Some results of tils flutter program of wings in close proximity are 
shown in the next figure (5). The effect of biplanar separation 
distance between two lifting surfaces is presented as a function of 
equivalent airspeed over the Mach number range from about 0.4 to 1.30. 
There are two predominant effects to be noted as the separation 
distance between the two wings is reduced. First, the Mach number at 
vhich the minimum flutter velocity is encountered shifts to progressively 
lower values. For example, along the boundary for the single wing, 
which is depicted by the solid curve, the Mach number corresponding 
to the minimum flutter speed (or bottom of the flutter bucket) occurs at 
a value of about 0.75. and for two wings with a separation distance 
h/cf * 0.27 the minimum flutter velocity occurs at a Mach number of 
about 0.1»5. Tbe second effect resulting from reducing the separation 
distance between the wings is the reduction in equivalent airspeed 
or dynamic pressure at which flutter occurs over most of the Mach 
number range. For example, for the lower range of Mach numbers a 
reduction of about 50 percent is seen to result. Perhaps, more 
importantly, are the supersonic results since estimates for the space 
shuttle mission have the vehicle encountering maximum dynamic pressures 
of about 500-600 psf at a Mach number between 1.2 and 1-3- Again, 
the more flexible lifting surface is seen to flutter at lower values of 
equivalent airspeed or dynamic pressure as the wings are spaced closer 
together. The exact amount of the reduction is not known since no 
flutter data was obtained for the single wing in this Mach number range. 
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CONTROL SURFACE INSTABILITIES ON LIFTING-REENTRY VEHICLES 


Another research program was conducted at Langley Re ® ea fch Renter 
(ref. 16) with the final intent being to determine whether a control 
surface instability might exist, for lifting-body reentry vehicle 
configurations at hypersonic speeds. One of the ccnfigurations^tu xe 
was a 1/ 30-scale model of the HL-10 vehicle and is shown in • 

The photograph at the left shows the model at angle cf attack, with 
trailing-edge control surfaces deflected for trim, mounted in the 
M = 15 Langley hypersonic aeroelasticity tunnel The program wa ® 
conducted to determine the aerodynamic in-phase or stiffness coefficient 
and the out-of-phase or damping coefficient for the flaps by the free 
oscillation technique. The variables of the program included vehicle 
angles of attack up to 25^, flap deflection angles up to 60°, and flap 
reduced frequencies from 0.01 to 0.0L. The schlieren photograph on the 
right hand side of figure 6 illustrates the flov field associated with 
the model at 25° - of attack at a Mach number of 15.4 and 

Reynolds number per foot of 7.5 x 10 6 . The flaps are deflected 35°. 
(Trim 8 with vehicle at a = 25° is about 8 = 25°.) Note the 
coalescence of the bow shock, the boundary-layer separation wedge 
induced shock, and the flap induced shock intersecting below but 
in the vicinity of the flap trailing edge. This intersection point 
was observed to move upstream vith increasing flap angles, and to 
unsteady for flap angles greater than about 30°, and become more 
violent for larger flap angles. This would not only affect tte 
stability of the flaps themselves, but also their effectiveness 
control devices. 



CONTROL- SURFACE INSTABILITIES ON LIFTING-REENTRY VEHICLES 



aerodynamic damping results 

Some control surface aerodynamic damping results from reference 16 
are presented in figure 7 . The results are presented as a function ot 
an aerodynamic damping coefficient with flap deflection angle. The 
experimental data depicted by the solid circular symbols indicate 
a destabilizing trend of decreasing values of aerodynamic clamping with 
increasing flap deflection angles. Negative aerodynamic damping values 
were measured for flap angles greater than about 37°. In contrast, 
modified Newtonian theory calculations shown as the solid line predicted 
a stabilizing trend and positive values of aerodynamic damping over 
the range of aeroelastic parameters tested. Other similar two- 
dimensional results obtained, but not presented here, showed „ood 
agreement between theory and experiment, suggesting that t.iese 
destabilizing results are associated with interactions between the 
three-dimensional flow and thick trailing-edge flaps. However, the 
exact cause of the unstable phenomena was not determined. 











ANTICIPATED FLUTTER RESEARCH 

In summary, figure 8 outlines some of the anticipated areas where 
continued effort is needed in order that sufficient flutter technology 
will be available in a timely manner to ensure a safe and reliable 
space shuttle vehicle* 

In the area of planar lifting surface flutter (l), component model 
wind-tunnel studies should continue to he conducted on candidate 
space sLuttle configuration designs. Flutter data should be obtained 
over the entire flight regime to provide reliable trends including 
the effects of blunt wings operating at high “^les of attack. The 
transonic and maximum dynamic pressure regions of the flight should 
receive priority. 

The stall flutter (2) problem may need both an analytical and ex- 
perimental effort if the wing torsional stiffness is marginal. 

Wind-tunnel simulation of the transition maneuver from high to low 
angles of attack should be conducted on aeroelastically scaled wing 
models simulating Mach number and, in so far as possible ,Reynols^_ 
number is essential. A parallel theoretical effort should be conducted 
to gain a better understanding of the effects of oscillation amplitude 
and the rate of change of angle of attack on stall flutter. The 
final objective being to develop a theory capable of analyzing stall 
flutter to the degree needed for prediction of wing amplitude response 
and the associated induced wing stresses during transition through the 
stall region. 

For the interference flutter area, an evaluation is needed ofsome 
of the newer, more sophisticated analytical methods such as the method 
of Laschka (ref. 12) as to their applicability to the biplanar unsteady 
problem. For the higher speed regimes a quasi-steady theoretical 
approach seems both Justifiable and feasible. 

For control surfaces (U) a better understanding of the buzz phenomena 
is indicated including three-dimensional effects. Some relate 
problem areas for control systems are found in the field of aero- 
servoelasticity . 

The complete vehicle configuration flutter problem (5) warrants 
considerable attention. As with any new vehicle design, an aeroelastic 
proof-of-concept testing program is essential. For the 2-stage con- 
figuration due to the nature of its size and complexity, new, and Perhaps 
novel, wind-tunnel testing techniques will have to be developed. Test 
programs, as well as testing techniques, should be developed as soon 
as possible. Design and construction of complete aeroelastic models 
should be undertaken by the contractors as soon as design concepts 
become fixed. Testing programs, utilizing HASA and other facilities. 





ANTICIPATED FLUTTER RESEARCH 


1 PLANAR LIFTING SURFACES 

t Component wind-tunnel studies 

2 STALL FLUTTER 

• Wind-tunnel simulation of transition maneuver 

• Develop analysis for prediction cf boundaries, 
amplitudes, and induced stresses 

3 INTERFERENCE FLUTTER 

• Comparison between existing analytical methods 
and new experimental results 

4 CONTROL SURFACES 

• Better understanding of buzz phenomena 

• Aeroservoelastici ty 

5 LIFTING SURFACE - BODY 

• Development of testing techniques 

e Complete aeroelastic model wind-tunnel studies 
e Refine and advance analytical methods 
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therefore, could he completed and the results in-hend prior to hardware 
production. This would allow the results, as they might impact the 
design, to be incorporated into the final configuration rather than 
after-the-fact as detrimental weight fixes or penalties. 

A parallel analytical effort should share the high priority of this 
research effort as the current state of the art for lifting surface- 
body flutter is the least adequate. 
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STATE OF THE AFT FOP PANEL FLUTTER 

AS APPLIED TO SPACE SHUTTLE HEAT SHIELDS 

Sidney C. Dixon and Charles P. Shore 

HASA-Langley Research Center 
Hampton, Virginia 


SUMMARY 


Heat shields considered for the space shuttle generally have 
flexible supports and may have orthotropic stiffness properties. The 
state-of-the-art for the supersonic flutter of this type of panel is 
briefly surmarised. Results for the most general support conditions 
are not available, but considerable work has been done for panels 
with elastic edge restraint on two opposite edges. Results from these 
investigations which indicate the effects of elastic edge restraint, 
panel stiffness orientation, unequal support stiffness on apposite 
edges and elastic coupling to an elastic substructure are discussed. 
Four analytical methods presently considered for analysis of panels 
with arbitrary support on all four edges and experimental programs 
required for verification of these techniques are indicate*.. 
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INTRODUCTION 


Heat-shield designs often involve a lightweight outer skin; 
the major design problem is one of supporting thisoutersK 
it will survive in the airstream and, at the same time, will accommodate 
the thermal expansion due to the large temperature differences between 
shield and primary structure. The design techniques used to solve the 
problem ge» -ly lead to relatively stiff panels connected to a primary 
structure by flexible supports as indicated in the figure. The outer 
skins are often stiffened, for example corrugated sheets. wd t^s may 
have orthotropic stiffness properties. In some designs the supports are 
so closely spaced that they provide essentially continuous elastic 
support, whereas in other designs the spacing is so large that the 
attachments must be treated as discrete line or point supports} often 
the side edges are unsupported. Virtually no flutter^ result sexist 
for panels with such edge support. However, some of .no flutter charac- 
teristics of corrugation stiffened panels with continuous elastic 
support on two opposite edges and simple support on the other two edges 
have been determined, and both theory and experiment reveal a highly de- 
stabilizing effect of the elastic edge support (refs, l to 3). Thus even 
though heat shield panels may be relatively stiff, the influence of « e 
elastic support way lead to panel flutter problems. Panel flutter 
instabilities are usually of limited amplitude, however very rapid 
fatigue-type failures have occurred in wind tunnel tests (refs. 1 ana 
k). Hence, it appears necessary to design heat shield panels that will 
not flutter. 

In the present talk the influence of elastic edge support and 
elastic coupling to an elastic substructure on the supersonic flutter 
of panels will be discussed for panels with simple support on two 
opposite edges, as most information pertains to this type of panel. 
Proposed or needed investigate ans for panels with arbitrary edge support 
on all four edges will be indicated. All theoretical results to be 
presented were obtained from the two-dimensional steady aerodynamic 
app rox imation which yields accurate results for a wide range of mch 
numb er and panel geometry for panels with nondeflecting edges. The 
ranges for which this approximation yields reliable results for panels 
with edges, however, has not as yet been ascertained. 
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TYPICAL HEAT SHIELD CONFIGURATIONS 






SIGNIFICANT FLUTTER PARAMETERS 


Sane of the significant paraaetew that should be | 
flutter of heat-shield type panels indicated in .IKEeil 
nuabers in parentheses indicate pertinent ref ^ces ^h describe 
soilc of the investigations that have considered the effects of the 
“«£ paraaeters; the list of references is in no ^ ****** 
recent comprehensive review of the state-of-the-art for pan 
is given in reference 5. 

The effects of most of the parameters listed have been studied in 
o»st detail using approximate aerodynamic theory. Although such progress 

S S ,S «r««nt with «P«i r t U 

many limitations in the analyses still exist. ’ ^ have 

innlane load have usually been restricted to uniform loads „*.—*„*. 

Sss 
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EFFECT OF ELASTIC DEFLECT IONA1 EDGE RESTRAINT 
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The effects of continuous elastic deflectional edge restraint on 
two oooosite edges on the flutter of corrugation stiffened panels ^re 
SEw?iSthI figure, which was taken from reference l6. The abscissa 
narameter^ 1C iTa nondimensional inflectional spring parameter. In 

the a ordinate D flutter parameter b is the panel width, D, the F 

. . ■ *• 


blic - — a J n 

stiffness in the direction of air flow, and 

. . n orf fiTl rtfl 


is the A panel 

twisting stiffness. For'the panel orientation 1 * jhown, where the 

SBsrjstr r 

The panel stiffness in the cross stream direction is D 2 and for 
panel orientation shown Dg »D^ • 

The solid curve was calculated from the normal-mode Galerkin 
solution of reference 3 for the leading and trailing edge simply __ 
supported and C = 7; curves for smaller values of C would be a 

the curve shown* 

The symbols represent experimental data for values of C ranging 
from 5.5 to 7. The panels had a wide variety of side edge-support 
Conditions from very weak supports, with the ends of the corrugations 
SSSStfS vS "Song Z?oA where the ends of the corrugations 
wereattached to a rigid substructure. The Paneiswerenot designed 
specifically for verification of theory and in fact most of the tests 
wSe Se^fore the theory was available. The deflection^ stiffnesses 
at the panel edges were not measured and were difficult to define 
analytically. 

The data are seen to verify the trends of the theoretical curve 
and illustrate the fact that corrugation-stiffened panels tested to 
nutter* .t dynamic pr*.ur* far lean than th. - 
by the conditions of all edges simply supported. Note that th- larg 
experimental value of the flutter parameter is nearly an order of 
magnitude less than the theoretical value for si; -xLy supported edges. 
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EFFECTS OF ELASTIC DEFLECTIONAL EDGE RESTRAINT 



EFFECT OF STIFFNESS ORIENTATION 
ON FLUTTER BOUNDARY 


The effect of stiffness orientation on the flutter boundary for 
a souIJe corrugation-stiffened panel is shown in the figure which was 
taken from reference l6. The abscissa parameter is the ratio oi t 
JSS StS”tiffn. SS Dj2 to tb. uretr of the »oodin^otiff- 

n The flutter ^ parameter also contains Dp . R-sui 
presented 2 for 1C = oo, which corresponds to simply supported edges, 

s-ai-ss ^ to 
s2£S*sj 2 srss ^ 

is the sine function and the governing partial differential equation tnus 
reduces to an ordinary differential equation with J * 

The exact solution is independent of the stiffness parameter . 

The salient fact revealed by the results shown is that a pel with 
corrugations parallel to the air flow is much lesssusceptible to 
flutter than when the corrugations are normal to the air flow, in i , 
flutter has not occurred in the few tests that have been made on ortho- 

tropic panels oriented with the direction of SdT^ 

narallel to the direction of airflow. It should be noted that this _ ^ 
orientation corresponds to the orientation of most proposed orthotropic 
heat shield designs. 

Another significant fact revealed by the results is that the effect 
of elastic deflectional edge support is much less ^ no ® c !j w J® 
corrugations are parallel to the air flow, as compared to the effect 
when the corrugations are normal to the direction of air flow. 




EFFECTS OF UNEQUAL ELASTIC 
DEFLECTIONAL EDGE RESTRAINT 

The results presented in the previous two figures were obtained 
from the analysis of reference 3 which has certain limitations. 

This analysis has recently been extended at the Langley Research 
Center to include the effects of normal inplane loads, damping, elastic 
torsional and rotational restraints, and unequal elastic restraint on 
opposite edges: the restriction.. that.. two. opposite edges be simply 
supported still exists. Preliminary results indicating some effects 
of ureaual deflectional springs are shown in the fien^e in te f m f of 
the flutter parameter and the ratio of the spring stiffness at the 
leading edge to the spring stiffness at the trailing edge; the side 
edges are simply supported. Results are presented for a square plate 
with stiffness ratios D^/D^ = 0.05 and the maximum 

flexural stiffness D_ is in the direction of air flow. Results 
are presented for val&es of the spring stiffness parameter of 

1 5 and 10 . For K-~ = 1 the effect of uneq'.al springs is slight, 
and the results indicate a minimum flutter dynamic pressure for 
of about 0 . 8 . For of 5 and 10 the flutter dynamic pressure 

increases roonotonically with increasing values of . The value of 

the flutter parameter for = 10 and = 5 is whereas 
for K_ =5 and = 10 the value is 130. Thus, these preliminary 

results 1 suggest a panel is less susceptible to flutter if the larger 
spring is at the leading edge. The results also indicate that a 
significant difference in spring stiffnesses is required to affect a 
significant influence on the flutter results. 








effect of elastic coupling to 

ELASTIC SUBSTRUCTURE 


Some effects of the elastic coupling of two identical, isotropic 
1 el plates simply supported on all edges with supersonic flow over 
only one surface are shown in the figure. The results were taken from 
reference 13 and were calculated from a two-mode Galerkin solution. 

The results are presented in terms of a nondimensional spring stiffness 
parameter and the flutter parameter used previously. For k - ° the 
results correspond to only a single plate. As can be Been, increasing 
the spring stiffness parameter to 10 is highly destabilising but a 
further^increase to a value of 20 is very stabilising. Even more 
dramatic variations with k result when the effects of normal inplane 
loads are considered. It is apparent that the elasticcouplingt 
elastic substructure can significantly affect the flutter characteristics 
of panels and should be considered in analyses of hoat shields. The 
analysis is presently being extended at the Langley Research Center to 
obtain an exact solution that includes the effects of orthotropic pane! 
stiffnesses and elastic edge support at the leading and trailing edges 
for both panels. Both theory and experiment reveal that panels on 
continuous elastic foundations can experience subsonic flutter and 
divergence (refs. 13, 27, and 28); hence, extension of the analysis to 
account for these phenomena is needed. 
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WIND TUNNEL TESTS CP HEAT-SHIELD TYPE PANELS 


As a part of general studies of heat— shield type panels several 
Models were tested in the Langley 9-by 6-foot thermal structures 
tunr.el at a Mach number of 3.0, dynamic pressures from 1500 pef to 
3000 psf, and stagnation temperatures from ^50*F to 66o“F (refs. ?9 
and 30). The panel shown in the figure is representative of one of the 
models tested. The outer skin was a thin (.01 in.) corrugated sheet 
which was attached to an elastic substructure by flexible supports. 

The other models also had corrugated skins, but uhe means of attachment 
to the substructure differed from the method shown. The models were 
tested with the corrugations both parallel and normal to the direction 
of air flow. No flutter was detected in any of the tests. It should be 
noted that the test conditions did not reproduce the environment for 
which the heat shields were designed; the temperatures were far too 
low, but the excessive dynamic pressure provided a rather severe test 
of the structural integrity of the models. For one model the corruga- 
tion depth was such that performance was marginal in res.oting thermal 
buckling during radiant heating tests, yet the corrugation was stiff 
enough to perform satisfactorily in the wind-tunnel ttsts. 

No vibration data were obtained for the models and no vibration or 
flutter calculations have been made. Vibration results for similar 
panels are presently being calculated at Langley Research Center by 
finite element methods (NALTRAN). Flutter analyses with existing 
theories would not be realistic as all edges of the panels were supported 
elastically which violates the restriction in existing analyses of 
two opposite edges singly supported. 


REPRESENTATIVE HEAT SHIELD USED IN WIND TUNNEL TESTS 





FLUTTER RESEARCH FOR PANELS WITH 
ARBITRARY SUPPORT CONDITIONS 


2 one of the anticipated flutter research effort YL* 

arbitrary suoport conditions on all four edges are indicated xn the 
fi-rure. ^he* numerical techniques and nethenatical modeiine required 
to ^obtain accurate flutter predictions for complex heat s^eld P* 01 
vrill probably not be amenable to parametric studies of the many 
variables entering the problem. Hopefully, less rigorous analyses 
will be developed that are suitable for generating large amounts of 
information sufficiently accurate to establisn reliable trends. 

Application of at least four analytical approaches is Presently 
anticipated. Flutter results for relatively single panels have been 
obtained from analyses based on finite element cechniques (refs. 31 
and ~< 2 ) and finite difference techniques (ref. 13)* -‘■t 1S presently 

SrtiSpa^d that both of these techniques will be applied to the flutter 
of panels with arbitrary support conditions by personnel at the Langley 
Research Center. The HASTRAN system will be used for the finite element 
calculations; development of a finite difference analysis I tes just begun. 
The vibration characteristics of heat shield panels will be studiedin 
considerable detail at Langley, and both theoretical and 
modes and frequencies will be obtained. These results will be ay^lable 
for normal mode type analyses. Lagrangian multipliers were applied by 
Dowell (ref. 3*0 to the flutter of multibay panels, and he is presently 
developing this technique fer panels with arbitrary support conditions. 

The emphasis of the indicated analytical effort is on obtaining 
the ability** to account for (l) the orthotropic stiffnesses of complex 
stiffened panels and (2) the effects of arbitrary support conditions 
on all fou? edges. However, it should be emphasized that tne e^ets 
of the various parameters listed previously as well asthe effects ®f 
deflecting panel edges on the accuracy of existing aerodynamic theories 

should also be assessed# 

Two types of experimental programs are anticipated; tests on 
relatively simple models to check theoretical trends, and tests of models 
representative 5 of actual heat shields. The latter ‘ J2* "J* 

establish the accuracy of the theory but would also I^^ef lutterdata 
for proposed hardware. The test programs are still in the preliminary 

planning stages# 
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FLUTTER RESEARCH FOR PANELS 
WITH ARBITRARY SUPPORT CONDITIONS 

ANALYTICAL EFFORT 

1. FINITE ELEMENTS 

2. FINITE DIFFERENCES 

3. NORMAL MODE ANALYSIS 

(a) CALCULATED MODES 

(b) MEASURED MODES 

4. LA GRANGE N MULTIPLIERS 

EXPERIMENTAL EFFORT 

1. SIMPLE MODELS 

2. MODELS REPRESENTATIVE OF ACTUAL HEAT SHIELDS 
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VIBRATION RESULTS FOR PANEL 
SUPPORTED AT FOUR POINTS 


Two of the numerical techniques proposed for analysis of flutter 
of panels with arbitrary supports have recently been applied to th 
vibration of an isotropic plate supported at four points (refs. 35 
and 36 ). The figure shows theoretical and experimental results for 
the fundamental frequency of a square plate symmetrically^ supported at 
four points. The ordinate parameter is the nondimensional panel 
frequency where a is the panel length, m the panel mass per unit 
area, and D the panel bending stiffness. The abscissa parameter is 
the ratio of the support position to the panel half length. 

The solid curve represents theoretical results from a two ®° <ie 
Lagrangian multiplier solution (ref. 35), the dashed curve represents 
results from a finite difference solution (ref. 36 ), and thesymbols 
Spresint experimental data (refs. 36 and 37). The agreement between 
experiment aldat a and theoretical results from both types of solutions 
i?U however, the results obtained by lagrangian multipliers are 
in better agreement for support positions near the plate center. 

The results indicate that the magnitude of the frequency is strongly 
affected by support position. The maximum frequency occurs when the 
supports are about halfway between the center of the plate and ^epla.e 
edge, and is about a factor of two times the frequency for the supports 
at the plate center or at the corners. Since flutter results are 
strongly related to panel vibration characteriittics, it appears reasonable 
to JXSe that support position will also have a significant influence 
on the flutter of point supported panels. 


EFFECTS OF SUPPORT POSITION ON FUNDAMENTAL FREQUENCY 

OF A SQUARE PLATE 



CONCLUDING REMARKS 


Results for the supersonic flutter of orthotropic panels with 
elastic edge restraint have been presented. Both theory and experiment 
revealed that elastic deflectional edge restraint is 

Theoretical results indicated that orthotropic panels with the air flow 
in the direction of the maximum stiffness are much less prone to flutter 
than oanels with the air flow normal to the direction of maximum 
SSfSST Vfect of .lMtlo coupling to nn .lantie aubatructur. 
«s shorn, to be either stabilising or destabilising depending on the 
stiffness of the elastic coupling. 


Existing theories are restricted to panels with continuous elastic 
restraint on two opposite edges. In order toproperlyconsiderthe 
flutter of heat-shield type panels, analytical capability is needed for 
nanels with arbitrary support on all edges. Four possible analytical 
^=2s iSSatsT finit. .lent,, fljlt. 

mode techniques, and Lagrangian multipliers. Two types of experimental 
programs were indicated: tests on relatively simple models to check 

theoretical trends, and tests of models representative of actual heat 
shields* 
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ABSTRACT 

In the last ten yeara, progress in the linearized theory of 
unsteady aerodynamics has followed two courses, (a) the improvement 
of the numerical orocedures for isolated lifting surfaces, and later 
(b) the extension of those methods to include the effects of mutual 
interference between two or more surfaces and bodies* This paper 
gives several examples of recent progress to illustrate that the 
prediction of steady or unsteady airloads on some complicated 
configurations, such as components of the Space Shuttle, can now be 
performed routinely and accurately. However, technological gaps 
exist and required research from the Space fliuttle viewpoint is 
delineated* 
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SECTION I 


INTRODUCTION 

Aeroelastic and structural dynamic considerations form a part of the 
moat significant design conditions for large, very flexible and very light 
structures . Examples of problems or technical areas are flutter, divergence 
and control surface reversal and gust response. Others involving feedback 
or active control are load alleviation and mode stabilization. The accurate 
prediction and the reliable and safe control of the above phenomena require 
the capability to predict both the aerodynamic and structural dynamic 
characteristics with considerable confidence. This paper is a rapid review 
and evaluation of the state of the art in unsteady aerodynamics pertaining 
to various potential Space Shuttle designs and concludes with some comments 
and recommendations in this regard. 

Proposed Space Shuttle configurations encompass many aspects which had 
to be ignored in earlier calculations of unsteady aerodynamic loads. Large 
fuselages, folded wing tips, interference between adjacent or tandem surfaces, 
and control surfaces were configurations which were not amenable to rigorous 
analysis ten years ago. Linearized methods which account for many of those 
factors are now available and in day-to-day usage. They bring the unsteady 
aerodynamic load predictions for a complex configuration such as some 
components of the Space Shuttle to a capability equivalent to the prediction 
for isoleted planar wings of a decade ago. This paper will trace the 
ment of linearized unsteady aerodynamic methods from the early considerations 
of airfoils and planar lifting surfaces to the present day capabilities f>r 
aeroelastic and dynamic analyses of configurations with control surfaces and 
multiple surface or body interactions. The following sections will 
Isolated lifting surfaces, control surfaces, intersecting surfaces, and inter- 
acting (but not intersecting) surfaces. 
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SECTION II 


STATE OF THE ART 

Isolated Lifting Surfaces 

In the last ten years, tremendous progress has been made in the 
ability to use linearised theory to estimate unsteady aerodynamic airloads 
for aeroelastic and dynamic response problems* The progress in subsonic 
compressible flow methods rests on several important contributions. The 
first was made in 1940 by H.G. Kussner (Reference 1) who formulated the basic 
lifting surface theory which gives the now well known integral equation 
relating oscillating downwash to the surface integral of an influence 
function or kernel times unsteady pressure. Other important contributions 
were made in 1955 and 1959 by Watkins, Woolston, Cunningham, and Runyan 
(References 2 and 3), who developed forms of the kernel and mathematical 
procedures which then permitted the exploitation and practical application 
of the kernel function method. Since then, many other important contributions 
such as an early attempt to incorporate control surfaces have been made, and 
other improvements continue (References 4-8). The doublet-1 lttice method, 
of Albano and Rodden (Reference 9) is another method which bypasses some of 
the difficulties of the kernel function method and has shown remarkable 
accuracy and ease of computation which makes it particularly attractive. 

Progress in the prediction of unsteady supersonic airloads rests in 
the U.S. on contributions by Garrick and Rubinow (Reference 10) and, for 
three dimensional applications, the contributions of Pines, Dungundji and 
Neuringer (Reference 11) and Li (Reference 12) who suggested the Mach box 
approach about 1954. There then followed rapid developments of workable and 
practical Mach box methods for predicting the supersonic unsteady aerodynamics 
on isolated surfaces (References 13 - 16). 

Control Surfaces 


Methods to predict supersonic unsteady control surface airloads were 
investigated by Donato and Huhn (Reference 17) using Mach box techniques. 
Subsonic unsteady aerodynamic methods have been delayed because of complexities 
resulting from control surface hinge line and side edge singularities. How- 
ever, Berman, Shyprykevich , and Smedf j eld (Reference 18) have explored the 
use of the kernel function method for the subsonic region using approximate 
expressions for the aileron pressure series. Ashley and Rowe (Reference 19) 
have presented a more rigorous formulation for the full span aileron on a 
subsonic lifting surface using the very notable contribution of Landahl 
(Reference 20). Further concentration and progress should lead to tractable 
methods for predicting the subsonic airload on partial span control surfaces 
oscillating in a subnonic flow. The Albano-Rodden doublet-lattice method can 
be applied to wings with control surfaces with a minimum of additional effort 
and is currently being used in that manner at several companies and at the 
AFFDL. 




Intersfecting Surfaces 


Flutter prevention is a critical design conditicn for T-tall air- 
craft. This Is emphasized by a few cases of actual in-flight catastrophic 
failures. Thus, the considerable effort devoted to prediction of T-tall 
unsteady aerodynamics for subsonic flow which started in the early 19 
is not surprising. Stark of Sweden (Reference 21) in 1964 produced a 
reversed flow method using doublets of integrated acceleration potentials 
to obtain the integral equations. Davies of the United Kingdom 
22) extended his subsonic kernel function method to T-tails anu -eported 
extensive apollcatlons of his method also in 1964. Zwaan of the Netherlands 
(Reference 23) developed in 1967, a kernel method following suggestions given 
by Laschka (Reference 4) and Davies (Reference 22). Many versions of the 
subsonic T-tail kernel function method now exist in the U.S. although 
procedure still requires considerable work. A general purpose subsonic kernel 
function method which is to be a development of the non-planar kernel 
function first given by Ashley (Reference 24) now is beingdevelopedby North 
American Rockwell Corporation for the AFFDX. to handle T-tails, V-tails, 
cruciform tails, folded tips, etc., and is in the final rtages of develcpmen . 
Many Space Shuttle configurations would be covered by this computer program. 

The doublet-lattice method of Alba no and Rodden (Reference 9) can 
also be applied to Intersecting surfaces. Kalman, Rodden, and Giesing 
(Reference 25) have extended the method and have applied it to wing-pylon 
combinations, annular wings, wing and fuselage combinations, and wings in 
ground proximity as well as T-tails. An application to Stark s T ^ ail 
(Figure 1) and a comparison with other methods are given in Table 1. Agree 
went is excellent. 

Andrew (Reference 26) is now de\ eloping a supersonic Mach box 
method for multiple intersecting surfaces. This is an extension of earlier 
work by Moore and Andrew (Reference 27) and Donato and Huhn (Reference 17) 
which could handle folded wing tips using the Mach box method. 


TABU 1 GfNfRAt IZEO FORCES FOR THE STARK T TAIL 
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T nrprArtinp but not Intersecti ng Surfaces 

Examples In this category include some Space Shuttle configurations 
,ith toe booster and orbiter vehicles tied to each other fuselage-to- 
fuselage, producing tandem airfoils having interactions between ® J 

. n 8 flutter model investigations have indicated that the oscllia 

tails. Recent fl ® significant influences on the horizontal 

ting downwash shed by a wing has very s g interactions produce 

d a e l t lentil SSTSStS speeds. Flutter speed trends cannot 

^ "““cUy u^». thi. .lng-on-c.il 1. include. 

The AFFDL has applied the Rodden-Albano doublet-lattice method and 
the Albano kernel function method to a wing-fuselage-horizontaltailflutt 

model tested in the WPAFB five-foot subsonic wind tunnel (Reference 32). 

The Albano kernel Ju, ^g-Srizontal 

*■ “"fruits. 

and tail planes are permissible. 

jr: t vr 6 rr./.° r s: 

SrSs £s 

r.Tl- fregu-ncy r.tlo .“ST ‘c^cie. 

;r sa-iars SsS i.nsr r. aarsa £. 

angle of 45 degrees. Direct compaxs lane of the wing 

sari.* ■r^ssr-srs; 

caused thereby for wlng-fuselage-horizontal tail flutter. 
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FLUTTER 



FIGURE 2. 1HE AFFOl SUBSONIC WING-TAIl aUTTER 
MOOa 





Albano has also applied his kernel function method to swept back, 
constant chord flutter models tested at high subsonic and transonic speeds by 
Cornell Aeronautical Laboratory (Reference 36). Typical expet imeiual and 
computed results are given on Figure 3 for a given vertical and longitudinal 
tail position. Predicted results agree very well with the corresponding 
experimental data. This was found to be true for a wide range of ver 
and longitudinal tail positions. 

Tha AFFDL has applied the Rodden-Albano doublet-lattice and Albano 
kernel function method to several other configurations. Figure 4 shows 
lift coefficient results predicted by both methods for an oscillating wing- 
m 3 .hid. i. «. »“"» “ “ 
evaluative standard. Agreement between both methods is excellent. 

The Albano kernel function method has also been applied to the 
same AGARJ) configuration where the horizontal tail is displaced vertically 
a distance equal to 0.6 times the wing semi-span. Results are shown in 
Figure 5 together with coplanar data. The tail unsteay aero y 
coefficients due to wing deflections are lowered significantly. Fi «“ r ® 3 
also shows the effects when tail dihedral is also added to the vertically 
displaced tail. Tall coefficients due to wing defections “ e 
reduced with positive dihedral but are increased for negative dihedral as 
expected. Trends of these predicted coefficients appear to be correct. 

Figure 6 shows the trend of the real parts of lift coefficients 

r i*.s h xi ss £ £. 

lattice and Albano kernel function method are shown. Agreement amo g 
three sets of data is excellent. 






FIGURE 3. CAL TRANSONIC WING-TAIL aUTTER MODQ 
Aik • 0.4. AX/ b • 1.0 



FIGURE 4. REAL SECTION LIFT COEFFICIENTS DUE TO 

WING TWIST, M*0.8„ Sw_= 1.5. COPIANAR 
CASE V 



WING TA’L 



FIGURE 5. SECTION tlFT COEFFICIENTS OUE TO WING 
TWIST, M = 0.8, a 1.5, NON* PLANAR 
CASES. V 




FIGURE 6. INTERFERENCE EFFECTS ON THE IASCHKA 
SCHMID CONFIGURATION, So; - 1.0 
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Since one limit of tandem airfoils is a biplane configuration, the 
Albano kernel function method was applied to a biplane consisting of two 
identical AGARD wings with several vertical separations. Figures 7 and 
give section lift coefficients due to plunging and pitch of one surface with 
the other fixed. These lift coefficients increase as the moving wing 
approaches the fixed tfing. 





moving surface 


a .0,0.5, 0.3 


— i 1 '4 
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FIGURE 7. REAL UR COEFFICIENT DUE TO 
PLUNGE Of ONE SURFACE OF A 
BIPLANE, M • 0. 8. sw * 10 
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FIGURE 8. REM. LIFT COEFFICIENT 0UE TO PITCH 

Of ONE SURFACE Of A BIPIANE. 0.8, 
S OJ • 1.0 
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INADEQUATE WORK IM PROGRESS MARGINAL 


THIN PLANAR WINGS 


INTERFERENCE 




CONTROL SURFACES 


WING-HORIZONTAL TAIL 


THICKNESS, HUNTNESS, ANGLE-OF- 
ATTACK 




FUSELAGE “WING PYLON-STORE 


TRAN- SUPERSONIC HYPER- 
SONIC SONIC 


FIGURE 9. SUMMARY OF THE STATE OF THE ART 
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SECTION III 


SUM MARY, C ONCLUSIONS. ANP~ RECO MMENDATIONS 


The work of the last ten years in the linearized theory of unsteady aero- 

partial span control surfaces and, xn some .esp * efforts should also 

Research emphasis could fill these gaps quickly. Current efforts s " . 
place the supersonic Mach box method on a firm basis for wi^g-tailand in 
secting surface interference problems. However, the computer Programs are 
complex, and careful development and checking are paramount. Figure 
overall summary of our estimate of the current state of the art. 

The transonic area is generally the most critical area fro® icult 

and structural dynamic viewpoints but theoretical methods are most diffic 

Si r 

aerodynamic methods are also needed to avoid buzz phenomena. 

Subiects requiring attention and emphasis include control surfaces, 

4 . interference effects from large or thick bodies or fuselages, 

nL-sinuaoidal or transient type unsteady aerodynamic coefficientsforroot 
“cus forms of solution, and particularly, unsteady aerodynamic pressure 
XOCU8 iorm» Vfllldate theory or to form reliable estimates in design 
“easurements to ^ t ±oled above, research in the last decade has 

SSJbuUd . remarkable, practical and reliant tachnoloK, for predlctloo and 
avoidance of Space Shuttle dynamic and aeroelastic problems. 
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BUFFET AND AERODYNAMIC NOISE 


By Charles F. Coe 

Ames Research Center t N70-36 603 


INTRODUCTION 

Buffeting i.as been defined by some dynamicists to be the response 
of a structure to unsteady aerodynamic flow. It has also been defined 
as the input ioading to the structure by the unsteady aerodynamic flow. 
Without resolving which definition is more nearly correct, this paper 
will refer to buffeting as the random bending oscillations of overall 
vehicles and major structural components. Aerodynamic noise (which in- 
cludes boundary-layer noise) is caused by the flow unsteadiness, and is 
related to the surface pressure fluctuations causing buffeting and local 
panel response. 

The facing figure (#1) illustrates the main sources of flow un- 
steadiness causing aerodynamic noise and buffeting. The main sources are 
oscillating shock waves that occur at transonic and supersonic speeds, 
turbulent wakes, separated flow, attached turbulent boundary layers, 
vortices and possible flow separation from sharp cornered bodies, inter- 
ference flows, and transitional flow. The sketch depicts the flow over 
the now familiar MSC space shuttle concept in the launch configuration. 
However, the same sources of unsteady flow would occur on any configura- 
tion for both the launch and reentry phases of flight. 




BUFFET 


At present, there is no analytical method available for prediction 
of buffet Intensities. Consequently, such predictions must be based on 
experimental information generally obtained from wind-tunnel tears of 
scale models. 

Previous ej-perience from flight and wind-tunnel tests of aircraft 
and launch-vehicle configurations has shown that the severity of buffet- 
ing is largely dependent upon the extent of flow separation on the ve- 
hicle. Some of these previous data have shown the effects of airfoil 
profile on wing buffet, and that wings having thickness ratios greater 
than 10 percent can be expected to buffet at zero lift. All wings re- 
gardless of planform or profile will buffet when stalled. It can be 
expected therefore that currently envisioned space shuttle concepts will 
suffer from buffet, and that buffet loads will have to be carefully con- 
sidered for each candidate configuration. 

The first obvious question on shuttle vehicle buffeting concerns 
buffet of prospective wing configurations, particularly the MSC proposed 
straight wing. Results of some exploratory Langley tests reported in 
NASA LWP-872 are shown in the facing figure (#2) that confirm the expec- 
tation that a 12-percent-thick straight wing will buffet at zero angle 
of attack. These buffet boundaries for both a 3-percent-thick clipped 
delta wing and a 12-percent-thick straight wing indicate that the thin 
delta wing would be free of buffet at low angles of attack during launch 
but that the straight wing will buffet between M=0.8 and M»0.9. As 
expected, the boundaries indicate that both wings would buffet at high 
angles of attack. 


240 


WING BUFFET BOUNDARIES 



Examination of buffet boundaries is an appropriate first step, but 
boundaries yield no information on the potential severity of the problem. 
In order to gain some insight as to whether buffeting would be seriously 
detrimental to the straight wing, some very preliminary estimates of 
wing bending -moment fluctuations have been made, and are presented in the 
f acing figure (#3). The results on the left show the ratio of maximum 
predicted full scale bending-moment fluctuations to the MSC design static- 
load limit for angles of attack up to 12°. These results, which were 
scaled from measurements on a lt-percent-thick wing, indicate that for 
normal launch conditions (a < A°) the bending-moment fluctuations would 
be less than 2 percent of the design static limit. An extreme launch 
condition angle o£ attack of 12- would result in t .1 1° percent » ‘ 

limit fluctuation. Recent Langley teats reported in NASA LWP-872 yield 
a similar appraisal of the ratio of maximum peak amplitude of bending- 
moment fluctuations to static moment at a - 60 . These results, shown 
on the right, indicate a maximum ratio of about 15 percent. 

It should be kept in mind that the above-mentioned buffet intensity 
estimates are only i.irst approximations. Additional tests of an elas- 
tically scaled model of the MSC straight wing are now in progress that 
will soon yield considerably more accurate results on low angle-of-attack 
wing buffeting intensities* 

The estimation of the wing buffeting is critical to the wing design, 
but because of very large fuselages and empennages, wing buffet may not 
be the key factor responsible for the buffeting of the overall vehicle. 
Exploratory tests of a few launch-configuration concepts are therefore 

needed to appraise the effects of major geometric variables. Such tests 
are scheduled, and are noted in the summary of this paper. 


PREDICTED WING BUFFET INTENSITIES 


O 
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AERODYNAMIC NOISE 


The information needed to fully describe the unsteady surface pres- 
sure field assocfited with aerodynamic noise is as follows: first, the 

intensity or magnitude of unsteady pressure; second, the location and 
area of coverage on a vehicle surface; third, the frequency spectra; and 
fourth, a pattern of spatial correlation. The intensity of fluctuating 
pressures is the easiest of the above listed information to acquire, and 
is consequently the most voluminous among the available experimental 
data that can be used for preliminary estimates of space shuttle unsteady 
pressures. The facing figure (#4) shows typical broadband intensities 
of local pressure fluctuations for different types of flow and the ranges 
of Mach numbers where data are available. The attached boundary layer 
which has the lowest intensity pressure fluctuations is in general the 
most thoroughly described statistically for M < 3.5. Attached boundary 
layer data for M > 3.5 are very 3 parse and lacking in description. Transi 
tional flow which can result in high intensity pressure fluctuations has 
been investigated for intensity and spectra at only the two Mach numbers 
shown. Pressure fluctuations in regions of shock waves and separated 
flows have been studied quite extensively, and detailed spatial correla- 
tion and scaling information is nearing publication at Ames. 

Aerodynamic noise is unfortunately very configuration sensitive, so 
that eventually a model of the final space shuttle configuration should 
be tested to establish a complete description of the unsteady pressure 
environment. In the meantime some exploratory space shuttle tests are 
needed to obtain preliminary data in regions of interference and possible 
vortex flows. Additional general research is also needed to acquire 
more complete statistical description of surface pressure fluctuations 
beneath hypersonic transitional and turbulent boundary layers. 
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INTENSITIES OF LOCAL PRESSURE FLUCTUATIONS 



MACH NUMBER 




SUMMARY 


In summary, buffet and aerodynamic noise are so configuration sen- 
sitive that the moat critical research will be needed after final con- 
figuration selection. Some important interim research is in progress 
or planned, however, at NASA centers that will provide much needed pre- 
liminary design information. The key items of activity are tabulated 
in the fating figure (tf5). 

For BUFFET, exploratory tests are now in progress at Ames to de- 
termine the overall orbiter vehicle buffet loads of the MSC configura- 
tion during launch, and to estimate the unsteady loads at the booster 
attachment points, Elastic models are also under construction and will 
be tested to determine if the buffeting of a complete two-body launch 
vehicle will be significantly altered by changing wings from straight to 
delta planforms. Tests of an elastically scaled MSC straight wing are 
in progress that will yield more precise wing buffet information than is 
currently available. The planned MSFC investigation is to evaluate the 
possibility of critical aeroelastlc effects using a quasi-steady aerody- 
namic technique. The last item under BUFFET is for ad hoc tests of the 
final configuration. Ames is Indicated for support, and the question 
mark refers to uncertainties regarding contractor participation. 

For AERODYNAMIC NOISE, aside from the need for an eventual final 
configuration study that will fully describe the surface pressure fluc- 
tuations, exploratory tests are planned at Ames that will identify and 
measure regions of high-intensity flow unsteadiness on candidate con- 
figurations. Some general research programs at MSFC and Ames that have 
application to shuttle vehicles are continuing with the objective of 
improving the statistical description of surface pressure fluctuations. 
Unsteady pressures due to high altitude plume induced flow separation 
is being studied by MSFC. Basic studies of hypersonic transitional and 
turbulent boundary layers are getting underway at Langley and Ames. 
Lastly an investigation is in progress at Ames to determine the effects 
of panel motion on aerodynamic noise. This is a potentially critical 
item for application to the calculation of response of thermal protec- 
tion system panels to aerodynamic noise. 


SUMMARY OF CURRENT AND PLANNED 
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PREVENTION OF COUPLED STRUCTURE-PROPULSION 
INSTABILITY (POGO) ON THE SPACE SHUTTLE 

S. Rubin 

The Aerospace Corporation 
Los Angeles, California 


INTRODUCTION 

Many liquid-propellant rocket vehicles have experienced 
longitudinal vibration, nicknamed pogo after the jumping stick, due. 
to an instability arising from interaction of the vehicle structure with 
the propulsion system. These vibrations can produce an intolerable 
environment for astronauts and equipment, can overload vehicle 
structure, and can lead to loss of propulsion performance. A NASA 
monograph concerned with the prevention of pogo is well on its way 
toward publication and should serve as a useful guide to design for the 
Space Shuttle. 

The paper will briefly describe the nature of the pogo 
phenomenon, present a summary of past instabilities, relate the 
possible deleterious effects of pogo, and discuss the contents of the 
design criteria monograph. Finally, research tasks are recommended 
for support of pogo development activities on the Space Shuttle. 


The pogo instability phenomenon results from a dynamic coupling of 
vibratory motion of the structure of a space vehicle with oscillations in the 
propulsion svstem. The simplified vehicle schematic at the left of chart 
portravs a propellant tank, a feedline, and an engine. Imagine the engine 
tc be comprised of a pump, a pump discharge line, and a thrust chamber. 
When the vehicle vibrates longitudinally. the pump and the propellant in 
the -lcxible tank undergo oscillatory motions. These two motions produce 
o-cillatorv flow in the 'feedline and in the pump discharge line. The How 
oscillations lead to oscillations in engine thrust and in pressure at the 
pump inlet, which then act as regenerative forcing functions on the vehicle 
structure. Another type of propulsion- system behavior, involving an active 
tank-ullage pressurization system, may also lead to an instabihtv termed 
ullage-coupled pogo as opposed to the above engine-coupled p c 

A block diagram of this positive feedback process is also shown. 

The coupled structural-propulsion system is linearized to consider pertur- 
bations of the svstem variables to determine the degree of linear stabili.y. 

Of great significance is the relationship of natural frequencies of the 
structure with those of the propulsion system The upper- right diagram 
shows the variation of natural frequencies of structural modes A and B 
and of a single propulsion system mode denoted by f~. The frequency f 
changes primarily because the steady pressure at the inlets to pumps change 
with time and the degree of pump cavitation thereby also changes. Coinci 
dence of structural and propulsion natural frequencies strengthens the 
coupling which may lead to instability. Other important system character- 
istics are the gain and damping of the modes of the structural and propulsion 
svstem s. Pogo has been corrected in flight by introducing an accumulator 
into feedlines to maintain separation of significant modes of the structural 

and propulsion systems . . , 

The results of stability analysis can be expressed m terms of the 
variation of damping of the coupled system with time, as depicted at the 
center right. System damping tends to decrease at frequency coincidences. 

shown, the svstem mode involving the structural rooae B has a damping 
ratio Cb which dips, but remains positive. However damping for the system 
mode involving structural mode A becomes negative and the system there- 

fore becomes unstable in a linear sense. 

\ corresponding depiction of vehicle longitudinal acceleration appears 
at the lower right. When Cb dips. narrow-band random vibration in struc- 
tural mode B tends to intensify. Later in time during the period of 
instability involving mode A, there appears a slowly varying limit cycle 
controlled by nonlinear behavior. Modeling the coupled system tor accurate 
prediction cf the limit-cycle behavior is beyond the state of the art. 








CHARY 1 

NATURE OF POGO (ENGINE COUPLED) 



tTHRU ST OSCILLATION 



A compilation of information on engine-coupled pogo instabilities is 
shown on chart 2. Acceleration amplitudes at spacecrafts appear in the 
left plot and at engines in the right plot. The plots are^intended to be an 
overall depiction of pogo history and should not be used for detailed 
comparisons among vehicles because of differences in measurement 

locations. , . , , 

Bands 1-3 show ranges for certain consistently unstable * hor 

family vehicles. Points and bands 4-8 are for Titan-iamlly vehicle of 
the past. Among the Titan it vehicles. N- 1 1 is unique in that is represents 
the first attempt, to correct pogo with a feedline arcumulator devue. 
Unexpectedly vibration severity increased markedly, later it was oc 
nized that the system had been incompletely understood, cr* 

N-ll mission occurred because large chamber pressure gm. 1 ' 
during pogo caused shutdown of the engines. 

The Gemini-Titan vehicles were the first to have a pogo man- rating 
requirement and a limit of 1/4 g was imposed at the base of the Gemini 
capsule. Accumulators were incorporated into fuel and oxidizer feed 
of the Titan Stage I and the limit was not exceeded except on Gemini 
Titan- 5 for which the two oxidizer accumulators were improperly charged 
prior to liftoff. All Titan III vehicles other than Titan IIIC- < have been 
stable. Titan IUB and advanced Titans have accumulators to prevent 

instability. a man - rating requirement, an appropriate working- 

efficiency boundary taken from a proposed standard is shown in the left plot 
This boundary applies to short-term sinusoidal vibration of man m a ^back 
to-chest direction. The degree of impairment of man's working efficiency 
of course depends on the nature of the task and other aspects of the environ- 
ment. Note that the second Saturn V launch with Apollo 6 considerably 
exceeded this boundary; the first Saturn V vehicle did not pogo. Because 
of concern for the astronauts and for structural loading within the space- 
craft, accumulators were introduced into the a- 1C stage of later Saturn 
vehicles and the instability has not reappeared during S- 1C burn. 

Suspicious indications of vibration during S-2 burn of early Saturn V 
vehicles and the Apollo 6 experience during S- 1C burn led to addition 
of accelerometers to detect S-2 pogo on the next Saturn V launch (Apollo 8). 

A 10 g amplitude at the center engine toward the end of S-2 burn was 
rc corded in a local structural mode involving relatively large motion ot the 
center engine. This instability was avoided on later vehicles by shutting 
down the center engine one minute early. However, the se verity of lesser 
vibrations earlier in S-2 burn increased markedly in later flights. About 
3. 8 g was recorded at the center engine on Apollo 12 leading to considers* l n 
of system correction by means of a center-engine lox accumulator. Apollo 
13 Jas not corrected, and even though it was virtually identical to Apollo 12, 
the vibration rose to 34 g! As with the Titan II N-il vehicle, chamber 
pressure oscillations caused premature engine shutdown. T ^ vl . br ^ 1 
mode is so local that less than 0. 1 g was observed at the ApoUo 13 space 
craft. Here is a striking example of the great sensitivity of pogo to small 
system changes - six launches of Saturn V/ Apollo vehicles without large 
early S-2 vibration, then large vibration on Apollo 12, and enormous 
vibration on Apollo 13. Obviously the lack of a pogo problem on the fi - 
several launches is in itself an insufficient basis for a lack conce n. 

The recent flight of the first Diam.int B vehicle is exceptional in tha 
the frequency was high (second longitudinal structvral mode) and it is the 
only known occurrence of pogo during bu r n of aprea sure - fed « 
next paper deals with a new mechanism of feedback which may be impo 

on the Diamant B. 
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CHART 2 

POGO INSTABILITIES (ENGINE COUPLED) 
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Clearly the astronaut environment is a chief concern for the Space 
Shuttle. Man's sensitivity to vibration, and the levels of vibration which 
have resulted from pogo instability, strongly indicate the need to have 
stability in all modes involving significant motion in the spacecraft. 

Although the number of cases of positively identified structural and 
equipment failures is small, such failures are also a major concern 
Whereas man rnay often be the limiting factor in the ®y 8te ^’ £ ° r 
instabilities near engines (socalled minipogos) the 

for structural failure and for the possibility of premature shutdown of 
engines* 


CHART 3 

POGO EFFECTS 

» INTOLERABLE ENVIRONMENT FOR ASTRONAUTS 

AND EQUIPMENT 

• OVERLOAD VEHICIE STRUCTURE 

• ENGINE SHUTDOWN 
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A NASA Space Vehicle Design Criteria Monograph on pogo is 
currently being prepared by the writer under the auspices of the Langley 
Research Center. The monograph, currently in Reader s Copy form, is 
entitled "Prevention of Coupled Structure-Propulsion Instability (Pogo)"; 
the Task Manager at Langley is G. W, Jones, Jr. 

An introductory section states the problem, discusser the need for 
attention to it, and presents the purpose, scope, and approach of the mono- 
graph. Basically the aim is to prevent pogo by modifying the system if 
necessary to achieve adequate linear stability. 

The State-of-the -Art, Criteria, and Recommended Practices sections 
deal with the individual aspects shown on Chart 4. The Appendix deals with 
the assessment of the significance of the structural modes. Some 62 ref- 
erences to the literature are given. 


CHART 4 

DESIGN CRITERIA MONOGRAPH 

PREVENTION OF COUPLED STRUCTURE-PROPULSION INSTABILITY (POGO) 


. INTRODUCTION 
» STATE-OF-ART 
. CRITERIA 

. RECOMMENDED PRACTICES 

* APPENDIX 

• REFERENCES 


1 MATHEMATICAL MODELING 
PREFLIGHT TESTING 
STABILITY ANALYSIS 
CORRECTIVE DEVICES 
FLIGHT EVALUATION 
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The need for pogo stability has already been stressed. That stability 
must be demonstrated by a suitable combination of analysis and test, ^.mce 
the achievement of a stable coupled system is the objective, analyses 
should be performed with a linearized mathematical model. Experimental 
data must be obtained to establish values for those model characteristics 
which significantly effect the resonant frequencies, gams, and dampings 
associated with modes of the structural and propulsion systems. Great care 
must be exercised so that tests are conducted in a manner to insure that 
small oscillation behavior is observed. This is particularly critical for 
determinations of structural-system damping and of turfcopump dynamics. 

Structural modeling for consideration of pogo on the Space Shuttle 
will be more complex than in the past owing to major three-dimensional 
coupling of motions resulting from asymmetric geometry of the vehicle. 

Based on past experience, particular attention must be paid to modeling 
of propellant tanks. Although a number of advanced tank-modeling approaches 
have been proposed in the last five years, these have not been sufficiently 
evaluated to determine the most suitable approach. The need is for a 
method to deal with both longitudinal and lateral oscillations and the ef fee 
of outflow oscillations into the feedlines to the engines. Also, detailed 
modeling will be required for engine -compartment structures to treat 
motions of such items as pumps, line bends, and engine support beams. 

Modal testing will be essential to the experimental verification of the 

structural" system model. , . . 

The overall structural system should be described for the mouel in 

terms of normal modes of vibration. The Appendix of the monograph pro- 
vides a method for assessing the significance of the individual modes. 
Insignificant modes should be excluded to avoid unnecessary clutter in the 
analysis. The significant modes must be defined with considerable 
accuracy. The correctness of the model should be re-evaluated after modal 
analysis. The monograph recommends that modeling improvements be 
assessed and incorporated when they change a natural frequency by more 
than 5 percent or a structural gain by more than 10 percent. 

Modeling of the propulsion system can begin on the basis of known 
properties of the propellants, flow- line geometries, and steady- state 
pressures and flows of the system. Only pump inlet cavitation can not be 

estimated from standard information. . , 

Frequency-response testing of the turbopumps will be required. 
Uncertainties in results of past testing of turbopumps are believed to be 
largely the result of inadequate instrumentation and test procedures. 

Needed are pressure and How instrumentation to detect low-amplitude 
oscillatory behavior. Flow instrumentation has not been available and 
pressure instrumentation ranges have been based on steady and large 
transient requirements. Also needed are test procedures which provide 
sufficient data to define completely the turbopump dynamics withou. requir- 
ing any assumptions of turbopump behavior. That is, the turbopump shoul 
be tested as if it were a "black box" to be described in terms of measured 
properties alone. Tests of the feedlines and of the overall propulsion 
system will also be required. 

Pogo stability analyses which include the consideration of parameter 
uncertainties, should be performed early in the vehicle development process. 
The monograph refers to an approximate methods for an initial investigation 
of the possibility of incurring engine-coupled instability. With this method 
the significance of uncertainties in the system parameters can be initially 
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CHART 5 


CRITERIA 


• COUPLED SYSTEM SHALL BE STABLE 

• ANALYZE WITH LINEARIZED MATH MODEL 

• EXPERIMENTAL EVALUATION OF MAJOR CHARACTERISTICS 

« PRIOR TO FLIGHT, VEHICLE CONSIDERED STABLE IF 

• 26 DB "DAMPING GAIN MARGIN" j NOMINAL 

• 230° "STRUCTURAL PHASE MARGIN" I CONDITIONS 

• PROBABILITY OF INSTABILITY SMALL 

• IF NECESSARY. ELIMINATE INSTABILITY BY SYSTEM 
MODIFICATION 

« ACCURACY OF MATH MODEL SUBSTANTIATED BY VEHICLE 
RIGHTS 
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assessed. As the vehicle progresses through its development phases and 
becomes better defined, increasingly detailed models should be constructed 
and stability analyses performed. Each cycle of stability analysis will 
provide information for a next round of refinements in the modal, including 
thr possible consideration of corrective modifications for achieving stability. 
Moreover, the process will contribute to advanced planning for testing „nd 
ti.e earl,' ( onsideration of alternative corrective modifications for the 
•=•* and w-*int s of implementation, reliability, and influence on the function 
it elements ot the system. Analyses should continue with updated 

p.. r ot » *■ t •* rs and increasing refinement until flight data ano postflight analyses 
; „ve demonstrated acceptable stability of the space vehicle. 

Special gain and phase margins are defined in the monograph and the 
r. dn and 30 degree margins are minimum requirements for the nominal 
coupled system at each time of tlight. These margins are defined so that 
any analytical technique that identifies an instability can be employed. 

This is achieved by prescribing that stability must be maintained for clearly 
specified modifications of the system. Complex mode shapes for the 
"nonciassical” normal modes of the coupled structure -propulsion system 
should be obtained to provide reference information for review of flight data. 

It is recommended that the vehicle be considered stable when analyses 
performed show either (1) that the system is stable under extreme tolerance 
(worst-case) conditions or (2) that the probability of occurrence of an 
instability is less than an allowable value assigned on the basis of overall 
vehicle flight-worthiness considerations. A conservative approach is out- 
lined for determining whether the probability of instability during flight is 
sufficiently small. A Monte Carlo analysis of feedline resonant frequencies 
may be Involved. 

Close-coupled accumulators for engine feedlines should be considered 
in any study of potential corrective modifications to achieve vehicle stability. 
Appropriate for cryogenic propellants are the metal-bellows type developed 
for Titan III vehicles and the gas-injection type developed for the Saturn V. 
Consideration should be given in the initial design phases to allow for possible 
future installation of such devices. A method is given for tentative sizing of 
an accumulator during early development of the vehicle. Requirements are 
identified for a development program to integrate accumulators into a 
propulsion system. 

The initial several flights of the vehicle should be specially instru- 
mented to obtain data which can be used to verify the mathematical model 
used for pogo stability analysis. A practical method to determine experi- 
mentally the degree of stability in flight has not been established. Due to the 
lack of such a method, the absence of an instability on an initial flight is far 
from a guarantee that the vehicle is adequately stable. Therefore only a well- 
verified mathematical model can provide confidence in the inh' -ent stability 
of the vehicle. Special instrumentation is identified and data reduction tech- 
niques are suggested. 

If an instability occurs in flight, the relative amplitude and phase of 
oscillations throughout the structural and propulsion systems indicate the 
"shape" of the unstable mode of the coupled system. This shape may be 
distorted by nonlinear effects at high amplitudes of oscillation. When the 
system is stable and sufficient coherency exists among the random oscilla- 
tions of the system variables, spectral analyses should be performed to 
determine amplitude and phase relationships among variables. Coupled- 
system mode shapes determined from the flight data should be compared 
with analytical predictions for evaluation of the mathematical model. Natural 
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frequencies of structural modes can be readily detected through spectral 
density analysis of acceleration data. Resonant frequencies of feedlines 
have been more difficult to detect, but narrow-range pressure instru- 
mentation should improve chances of success, A technique is available 
for evaluating individual portions of the structural dynamic model in 
cases of purely longitudinal vibration. 



Pogo instability is potentially a major problem for every phase of 
rocket-powered flight of the Space ShutUe-_A-monogr.Hph.is being prepared 
to make certain that the best advantage is taken of past experience. 

There are two intrinsic limitations associated with current prac ices. 
First the coupled structural-propulsion system can not be ground tested 
to directly evaluate pogo stability. Unfortunately, it appears that certain 
conflicting r^uiremem. rule such listing So nn.Wsi. b. 

relied upon to obtain confidence in vehicle stability Wtore Uie i • 

Second, current practice does not permit the determination during 
flight of the coupled- system damping. It appears that the determination _ 
ot system damping in flight requires the application of an a.m c 
banco at various times of flight so that resulting oscillations ,r. •«.. 

modes can be discerned. Unless system damping can be deter; ted 
experimentally for initial flights, the mathematical model 
only basis for assessing the degree of system stability. 

Development of a reliable model demands considerable analysis and 
test activities throughout vehicle development. These activities must be 
carefully programmed into the planning to insure an orderly and timely 
development. 


CHART 6 
SUMMARY 

* 

• POGO IS MAJOR CONCERN FOR ALL ROCKET-POWERED PHASES 

• MONOGRAPH PROVIDES CRITERIA AND RECOMMENDED PRACTICES 

• RELIANCE ON MATH MODEL TO ASSURE ADEQUATE STABILITY 
EVEN AFTER SEVERAL FLIGHTS 

• STRONG INTERPLAY OF ANALYSIS AND TEST 

• POGO ACTIVITIES MUST BE CAREFULLY PROGRAMMED INTO 
DEVELOPMENT PLANS 

§ LIMITATIONS OF CURRENT PRACTICES 

• REQUIRE FLIGHT TO REALISTICALLY TEST COUPLED SYSTEM 

t INABILITY TO DETERMINE COUPLED -SYSTEM DAMPING IN 
FLIGHT 
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Four specific research objectives appear to be attainable and could 
contribute markedly to the technology for treating pogo instability. The 
first two listed on chart 7 would significantly improve our ability to 
define the propulsion system by test. The need for dynamic flow trans- 
ducers is great. In addition practical means for testing a turbopump as 
a "black box” should he defined so that the pitfalls of deducing its behavior 
from overall propulsion-system behavior are avoided. For example, 
sinusoidal pulsing tests for two hydraulic boundary conditions at the 
turbopump discharge, with instrumentation for inlet and discharge 
dynamic pressure and flow, can yield the parameters of a "four-terminal 
network" which describes the turbopump. 

A thorough evaluation oi propellant' tank modeling methods is 
necessary. A number o* modeling methods have been proposed, but their 
comparative i apahilitv tor pogo needs have not been explored. 

Because o: a potential v large payoff, methods of inflight disturbance 
to determine svstevr. damping should be studied. It would appear that 
di*turbanc es generated both within and external to the propulsion system 
should be considered. 


CHART 7 

TECHNOLOGY NEEDS 


• DEVELOP TRANSDUCERS FOR PROPELLANT FLOW 
OSCILLATIONS 

• DEVELOP TEST PLANS FOR COMPREHENSIVE MEASUREMENT 
OF TURBOPUMP DYNAMICS 

• EVALUATE PROPELLANT-TANK MODELING METHODS 

• INVESTIGATE METHODS FOR INFLIGHT EXCITATION TO 
MEASURE COUPIED-SYSTEM DAMPING 
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1 - INTRODUCT ION 

The vibrations known as "POGO effect" result from a coupling between the 
structure, the tanks and lines, and the motor of a liquid propellant rocket 
stage. One of the elements of the coupling is due to the forces applied on the 
structure by periodical variations of the combustion chamber pressure. This paper 
aims at laying stress on the influence that ae: dynamic phenomena in the after- 
body may have on the value and the very nature of these forces, (--&• 

After pointing out the main parameters of the problem, two particular cases 
will be briefly examined : a single nozzle rocket operating at low level (L. 17 
booster of the DIAMANT B launcher), and a multi-nozzle rocket operating at high 
altitude. 

2 - SETTING THE PROBLEM 


In the methods generally used for forecasting the POGO effet, the efforts 
resulting from the periodic variations of the chamber pressure are computed by 
taking into account only the thrust coefficient of the propulsive nozzle in 
vacuum, but neglecting the effect of pressure variations on the other parts of 
the propulsive stage outside the nozzle proper. This assumption may be quest ionned, 
particularly when the base cross-section A A is large, relative to the throat sec- 
tion A t . If we call <&h, A and 0^ the simultaneous pressure variations at 
the aft cross-section and in the chamber, 0 **a and o")y the thrust variations 
due to the forces on the aft section and on the intunnel wells of the nozzle 
(thrust in vaccum) , it can be written; (Fig. 2) 

, A* it* 

y Afc *5 {Si 

being the thrust coefficient in vaccum. 


With this relation it is possible to appreciate the validity of the 
assumption generally retained, already mentioned, regarding the relative order 
of magnitude of 3T A / STy for various frequencies. 


The A a / AfcCy depends only on the geometry of the after-body and the 

nozzle. In the case of L.17 (D1AMANT B), its value is 8 ; for the SATUHN 11 stage, 
relative to 5 engines, it is 80. 
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The estimation of the ratio ( ^ (•'A /^N- > is very diificult by 

theoretical means ; moreover its introduction in the forecasting the stability 
of the POGO coupling entails that, for each frequency, Us amplitude and phase 

arc known. 

A first semi-empirical approach is given hereafter for a single nozzle 
rocket operating in an air atmosphere, and a qualitative description of the 
phenomena that may take place in a multi-nozzle rocket operating at high 
altitude is indicated. 

3 - SINGLE-NOZZLE ROCKFT AT LOW ALTITUDE 

The phenomena by which it is possible to express qualitatively the pressure 
variations “Sh* in the aft section in correlation with the bft c, variations 

is as follows • any E\ t c variation modi! as, in a quasi-steady fl •/. the sepa- 

ration line between the jet and the surrounding atmosphere- ; it may liso change 
the nature and the Intensity of the vortices issued from the Jet. All these move- 
ments induce vibrations in the surrounding atmosphere , as a noise at the excita- 
tion frequency. 

Bench experiments o'? 'he L.17 motor, aimed at measuring the noise around 
this rocket, show that th-ive is an actual correlation between the pressure varia- 
tio.is around the rocket and these of the chauber. 

The noise generated by the nozzle is transmitted to the aft section ; its 
intensity decreasing with increasing distance, but if the aft part of the rocket 
it as a resonant cavity ; in that case. around the resonance lrequency of the 
cavity, the pressure variations in the tank rear bulkhead may be very important. 

Vibrations observed during the first flight of DIAMANT B seem to have to be 
attributed to such a phenomenon. 

Acoustical tests of the rear failing have been carried out showing that 
the frequency observed in flight is not far from that of the results 

obtained are however difficult to analyze, as the excitation level that may be 
reached on the ground is low compared with that actually met in flight. 

On the other hand, tests at reduced scall have been performed; they give 
a proof that the phenomenon actualy exists.Visualizations of the phenomena are 
prepared for the coming months. (Fit- • ‘ "-.'.i ’) 

The ratio ( ) actually measured, at about the resonance fre- 

quency, reaches 3. This figure should be considered bu* as an order of 
for the full scale rocket, as it is not possible to respect all the similarity 
laws, especially for the damping, at the 1/40 scale model. 

The analysis of the POGO effect on DIAMANT B has been performed by intro- 
ducing measurements made on the ground, regarding in particular the damping of 
the structural mode. It appears that, taking into account the measurements during 
firing at rest, the acoustical measurements on the flight skirt or measurements 
on reduced scale models with air jets, the open-loop gain is such that vibrations 
of POGO type should diverge. (Fig b) 


4 - 


MULTI -NOZZLE ROCKET IN VACUUM 


It may be considered that an effect similar to that of a cavity at atmosphere 
pressure appears with multi-nozzle rocket operating in vacuum. In effect as a conse- 
quence of the jet convergence the pressure in the cavity limited by the je ex aus 
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and the rear bulkhead may be rather large , as it has been proved by tests performed 
on reduced scale models. (Fig. 6) 

This pressure is piloted by the pressure in the combustion chambers, and the 
problem lies in the determination of ( ) as a function of frequency. 

This investigation has been recently initiated on a reduced scale model of a - 4 - 
nozzle rocket whose chamber pressure, supplied with compressed air, was modulated. 

5 - CONCLUSION 

The first ground tests carried out show that the aerodynamical phenomena on 
after-bodies may have a noticeable influence on vibrations of POGO type. The results 
already obtained on the ground on a single-nozzle configuration must be confirmed 
and made more accurate. It seems necessary to have experimental results available 
before ascertaining this effect on multi-nozzle configurations operating in vacuum. 
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INTRODUCTION 

The space shuttle mission requirements such as high reliability, 
minimum refurbishment, minimum weight, and novel environmental combinations 
impose singular considerations and requirements on the structural design. 
This is particularly true in the design of the thermal protection system 
(TPS). Unlike previous missile requirements the TPS must successfully 
withstand repeated applications of combined and separate thermal and 
vibration environments. Acoustic fatigue will therefore be a major con- 
sideration in the design of the TPS. The subsequent paper therefore 
presents some of the expected acoustic fatigue characteristics of TPS 
concepts, and defines areas where supporting research is needed to provide 
the required reliability and confidence in the TPS design. 
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Basic space shuttle TPS design concepts are shown on figure 1. They 
consist of the following: 

A. Concept A consists of a radiative metallic structure with 
insulation attached either to primary structure or secondary structure, 

but in either case, the insulation does not contact the radiative structure. 
Aerodynamic loads are carried into primary load carrying structure by stand- 
off supports. 

B. Concept B duplicates concept A except that the insulation is in 
contact with the radiative structure. The significance of this design will 
be shown later. 

C. Concept C is the bonded ceramic foam tyne insulation which is 
structurally capable of direct aerodynamic load transfer. 

D. Concept D is representative of an active coding system in 
which the tubes shown here represent the heat exchanger used for transfer 
of re-entry heat from the TPS structure to a rejection exchanger or heat 

sink. 
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Since the purpose of this paper is to compare the vibration response 
and fatigue characteristics of each TPS concept, TPS design concept A was 
chosen as a baseline and preliminary vibration response characteristics were 

calculated. These are shown in figure 2 . 

The overall sound pressure level 'db ref. 2 x 10' U microbar), overall 
G rms and duration of excitation for the Indicated mission environs- ~.t an 
shown for concept A. The levels were obtained from an acoustic s : -ctral 
and structural areal density extrapolation of Saturn V, S-IVB data. Since 
this is a spectral extrapolation, an increase or decrease in overall sound 
pressure level does not necessarily cause a proportionate increase or de- 
crease in overall Gw The objective of this figure is to indicate that 
for each mission cycle a severe vibration response will occur for a relatively 
short duration of time and a moderate response will occur for an extenaed 
length of time. This moderate response will occur simultaneously with a 
severe thermal environment. The levels shown for lift-off and MACH 1/MAX Q 
are approximately equal to the levels obtained on lightly loaded Saturn V 
upper stage structure for the same mission-phase. It should be noted however, 
that Saturn V structure was required to survive only one mission cycle 
environment whereas the TPS system must survive repeated application. 
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VIBRATION RESPONSE OF METALLIC RADIATIVE 
HEAT SHIELD CONCEPTS 


MISSION CNV 

OVERALL SPL DB 

OVERALL ORMS 

DURATION OF 
EXCITATIONS 

LIFTOFF 

157.5 

30 

10 SEC 

MACH 1 
MAX Q 

159.5 

40 

40 SEC 

RE-ENTRY 

152*3 

5 

100-1200 













Figure 3 provides a simplified comparative evaluation of the vibration 
response characteristics and fatigue life potential of the various TPS 
design concepts. Ihe graph illustrates the relative response RMS stress 
levels for the same acoustic excitation source versus time. The graph 
purposely has no values assigned to either axis since the purpose of the 
figure is to provide representative trends rather than absolute information 
which is unavailable. For purposes of illustration the acoustic excitation 
is assumed constant. 

Concept A, with its high area to weight ratio and low damping value 
would be expected to have the highest vibration response (G^) , and result- 
ing RMS stress level for well coupled acoustic excitation. This concept 
would therefore have the least acoustic fatigue mission life potential. 

Concept D with its lower area to weight ratio should have a lower 
response with resulting lower RMS stress. Therefore concept D should have 
a higher acoustic fatigue life potential. 

Concept C is shown as a dashed line since no information could currently 
be obtained for the response of panels protected by this method. Since it is 
a rigidized fibrous material it would be expected to provide significant 
joint damping to any supporting metallic structure and therefore reduce 
structural response significantly. Its ability to survive even a moderate 
vibration environment and, particularly, low frequency vehicle dynamics is 
currently undefined. From the supporting structure standpoint, it does 
subjectively appear to incorporate a high acoustic fatigue life potential. 
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TPS SYSTEM RESPONSE TRENDS 


ACOUSTIC FATIGUE 
MISSION LIFE POTENTIAL 


CONCEPT 


FOTENTIAL 






'" ,G PAGE not FILMED. 


Concept B has an interesting and advantageous vibration response 
characteristic. Since the radiative panel is supported by a limp fibrous, 
highly damped insulation, its initial response to acoustic excitation should 
be significantly lower than an identical unsupported panel. Unless an 
adequate insulation can be developed that will withstand the dynamic 
bearing condition Imposed by the moderate panel response, the insulation 
will structurally degrade in the bearing areas. As the insulation degrades, 
the panel becomes less damped, the panel response increases, and the in- 
sulation degradation rate increases. The panel configuration will finally 
approach concept A and will have a vibration response and resulting stress 
equal to concept A. Stated another way, the acoustic fatigue life potential 
is a function of the rate of insulation degradation. This characteristic 
could be advantageously employed to provide extended miss. 'on life fatigue 
potential. 

All concepts could probably be designed to have an RMS stress level 
below the material endurance limit point 1. The resulting TPS system weight 
would, in all probability, be prohibitive. With TPS system weight a 
primary driving function, the response of large area to weight ratio structures 
will be above the endurance limit (point 2). Therefore, a finite fatigue 
life will result. If this fatigue life is greater than an acceptable 
refurbishment schedule (point 3), then any concept will be adequate. 

However, if a minimum refurbishment rate is a major objective (point 4), 
then TPS design concepts providing for induced panel damping at a minimum 
weight penalty must be developed. 
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To provide for the reliability, inherent in space shuttle 
objectives, the required research and development in the area of 
acoustic fatigue is presented in figure 4. 

Due to the novel design concepts required to accomplish TPS 
systnm objectives, existing data is non-existant. Before fatigue 
life assessments can be performed, representative forced response 
characteristics of TPS type structure must be established. 

Basic S/N curves for TPS material must be defined before re- 
liable estimations of fatigue life can be obtained. 

The effects of combined thermal and aerodynamic environments 
on the mission life of proposed TPS designs and materials must be 
determined. 

To provide for experimental verification of TPS system designs 
and shuttle component life tests, the relative forcing efficiency 
of re-entry aerodynamic excitation and standard plan wave, rever- 
beration chamber test modes must be determined, and test exagge- 
ration factors established, to provide reasonable qualification 
test times. 
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CONCLUSIONS 


Acoustic fatigue will be a significant design consideration for the 
space shuttle thermal protective system. Although not an insurmountable 
problem, significant research and development effort is required to obtain 
confidence in the reliability and mission life potential of any selected 

concept. 
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:r t structural integrit.. at tru Space- Strutt le 
;\:m the ope.rat to*\a l aspects at the* GSE in the 
launch, (2) siting and personnel hazard con- 
sidcr.it ion- with ropn L to the surrounding civilian womnumtA due to o.i 
sound being propagated into these areas. (3) inuunt iluctuating pressure 
environment during the sub. trails, super and hypersonic flight regime and 
it f s e fleet upon the structural integrity ei the Space Shut tie vehicle, (*) 
blast environment due to a catastrophic failure of the Space Shuttle, (5) the 
sonic boom associated with orbiter reentry and booster return flight, (6) the 
structural dynamic response characteristics of the Space Shuttle vehicle and 
(7) acoustic' fatigue and it's effect on the repeated use of the Space Shuttle 
vehic lo . 


Each of these problem areas must be examined in detail in order to realisti- 
cally assess their influence on the total mission reliability and success of 
the Space Shuttle system. 

Due to the limited time available, it is not possible to discuss all of these 
problem areas, therefore, 1 would like to briefly discuss two of these areas 
and present the results of the preliminary analysis of the acoustic environmen- 
tal characteristics of the Space Shuttle system. The two areas which will be 
discussed are the rocket engine environment developed during holddown and in a 
general sense the inflight fluctuating pressure environment developed during 
the trans and supersonic Space Shuttle exit (light regime. Further, 1 would 
like to present a specific aerodynamic phenomena that is ot a growing concern 
to those of us at MSFC and to discuss the possible influence of this phenomena 
on the Space Shuttle system. This phenomena is commonly t*et erred to at MSFC 
as the "High Altitude Plume Induced Flow Separation Phenomena”. 
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PRECEDING PACT: BLANK NOT HI. MED. 


Rocket Engine Acoustic Environment 

Since the initial formulation of the jet noise theory by Lighthill [ll many 
sophisticated theories have emerged. The most notable ones are those developed 
by [2] Ribner , L3] Phillips and [4] Williams. Each of these theories emphasized 
a different point of view. Lighthill by nis acoustic analogy approach developed 
a volume of acoustic sources in terms of a stress tensor whereas Ribner utilized 
the 'dimple-source" technique, i.e., monopole, dipole and quadrapole. Phillips 
presented his theory in terms of a converted wave equation the solution of which 
is restricted to high speed flows, i.e., Mach wave radiation. Williams in ex- 
tending the theory of Lighthi 1 l emphasized the effects of high speed convection 
and Mach wave radiation. Recently, Peters [ and Pao [ b I have introduced new 
aspects into the jet noise theory. Peter’s theory emphasizes a non- isent ror ic 
converted wave equation. Pao’s theory is actually an extension of the Phillips 
approach, the solution of which is valid for all transonic wind supersonic Mach 
numbers . 


All of these theoretical approaches requireddetai led knowledge of the mean and 
turbulent characteristics of the exhaust flow field and their practicability is 
restricted by the fact that they only apply to free undeflected exhaust flows. — 
These theories have provided valuable insight into the basic noise generation 
mechanisms of jet or rocket exhaust flows and much has been accomplished through 
their study. In the final analysis however, whenever it is desired to provide 
acoustic environmental estimates of actual situations, such as the Space Shuttle 
during holddown , some form of empirical techniques and/or model tests have to be 
employed. 

The comments which follow will be directed towards the techniques used for esti- 
mating the acoustic environmental characteristics of the Space Shuttle developed 
during holddown, (or static firing). These comments are not necessarily restric- 
ted to only this area but apply equally as well in the other areas of acoustic 
prediction; i.e., ground plane environment in the mid- and far-field areas, 
launch environment, etc. 

Definition of the decailed acoustic environmental characteristics of space vehicle 
systems is usually accomplished ia three distinct phases. The first phase is 
the preliminary environmental definition by using an empirical technique to ob- 
tain an estimate of the environment. Once this has been completed, the second 
phase can be initiated. This phase consists of planning, designing and conducting 
of model test programs utilizing scale models of the propulsion system to be used 
on the full scale vehicle. The third phase then consists of scaling the resulting 
model data to the full scale configuration. The second and third phases are 
normally performed as soot, as firm design information on tne propulsion system 
is known. Full scale testing is additionally employed to verify the estimated 
environments. 

Many empirical techniques have been devised throughout the years. All of these 
techniques relv in one manner or another upon the principle of an "apparent 
source model" C 7 ,8 ] and are based upon the principles of dynamic similarity. For 
the case of the Space Shuttle environment, this apparent source model consists of 
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a non-dimensional function which describes the location ot the apparent 
sources" and another nor -dimension.! I function tor determining the strength 
of the acoustic source. This technique was employed in predicting the bpacc 
Shuttle environment and the results are presented in Figure 1. A nominal 
thrust of 400.000 pounds per engine and 13 engines was used in computing these 
environments. Figure 1 depicts the overall sound pressure level, in B 
(RE* 0 00003 N/m2), as a function of Space Shuttle station. In the aft g_ 
of ihe ShuUle; tL environment is on the order of 1*5 to 170 « s and in the 
crew compartment area the environment is on the order of 152 to I5w dB *. Un. 
spectra distribution of the acoustic energy, in one-third octave hand sound 
pressure level form, at two selected Space Shuttle stations is presented in 
Figure 2. In the aft region of the Shuttle, i.e., Station 160, the energ> is 
primarily high frequency in nature and is distributed over a wide ranj. ot 
frequencies with the maximum energy being located m .he 5C0 H.rt 
octave band. In the crew compartment area, the overall sound pressure level 
ha* decreased to approximately 154 dB with most of the energy decrease occur- 
ring at the high frequencies and the frequency distribution ot the energy s 
centered about the 100 Hertz one-third octave band. These preUminary est - 
mates were obtained utilizing existing data from conventional rock*t ’ 
i e., those with chamber p-essure on the order ot bOO to 1-00 p.-i. The Sp 
Shuttle propulsion system will however, have chamber pressures operating at 
more than twice the pressures of these conventional engine., i.e., 3000 psi, 

and the effect of this increase chamber pressure on the inherent acoustical 
characteristics is not known at this time. This is one ot the factors which 
needs emphasizing during the second phase of the environmental detinitron, 
model test phase. Plans for this phase are in progress and the primary objec 
tives will be to: (1) assess the influence of the increased chamber pressure 
un the acoustical characteristics, (2) delineate the influence ot the non- 
conventional cluster arrangement of the rocket engines on the resu C 

pressure field, (3) study the effects of deflectors; i.e.. could be beneficial 
to TPS, (4) define the spectral characteristics of the radiated pressure tiel , 
(5) define the spatial correlation characteristics of the acoustic pressure 
field. Definition of the spatial correlation characteristics is required 
order to provide a realistic assessment of the dynamic structural response 
electronic design and qualification characteristics of the Space shuttle. Wt 
the reuseability concept being applied to the Shuttle* the necessity o . 

the %»elaLon characterise ies is even more acute from a fatigue standpoint 
Therefore, auequate assessment of dynamic response is mandatory (hence pro pe 
environmental acoustic tests need to be conducted) m order to insure t e 
structural integrity of the Shuttle and avoid costly redesign of Shuttle compo- 
nent due to fatigue. 
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Inflight Fluctuating Pressures 

Definition of the inflight fluctuating pressure environment is I* 1 ™** *"‘*™* y 
determined by empirical means and/or model testing. Delineation of the flue 
tuating pressure environment is also accomplished in several dist J nct 
similar to those employed in defining the rocket noise environment. The first 
ohase is to determine, through empirical means, the total energy an 
soectral content of the energy associated with several different types of aero- 
dynamic flow regions encountered during Space Shuttle flight. The fluctuating 
pressure environment is determined primarily by the external geometric, ^ charac- 
teristics of the Space Shuttle and therefore once the design has been establis e 
the second phase can be initiated. This phase consists of establishing the 
location of the many different types of aerodynamic flow regions as a function 
of the flight conditions to be experienced during Space Shuttle flight, i.e., 

Mach number, angle of attach, dynamic pressure, etc. This phas * is 
accomplished by flow visualization processes such as shadowgraphs, Schlerie 
or some other means. The third phase consists of conducting a scale model wind 
tunnel test of the Space Shuttle and the fourth phase is the scaling of the 
model test results to the full scale configuration. Full scale flight measure- 
ments are then required to verify and modify the design and qualification 
environments. The second, third, and fourth phases of the environmental defini- 
tion are normally initiated after the design has been frozen. 

Many experimental programs have been conducted in order to define the fluctua- 
ting pressure characteristics in distinct flow regimes associated aircraft or 
space vehicle flight. Most notably are those experimental studies of; Bull L9J , 
Bull and Willis [10] , and Serafini [ll] in regions of attached turbulent 

boundary layers; Kistler [12] , in regions of separated flow; Chyu and Hanley 
[13] and Coe [14] on an ogive cylinder in attached and separated flow regions 
and; Robertson [15,16] , in regions of protuberance induced environments. 

These experimental results plus numerous inhouse experimental results at MSFC 
obtained from a 4% scale model wind tunnel test program of the s * tu ' n y 
were used to obtain typical fluctuating pressure environment in four distinct 
flow regions associated with the Space Shuttle flight. These preliminary environ- 
ments were defined for regions of; (1) attached turbulent boundary layers, (2) 
separated flow, (3) oscillating shock and (4) separated flow/oscillating shoe* 
interaction, and the results are presented in Figures 3 and 4 in one-third 
octave band sound pressure level form. These environmental estimates are con- 
sidered valid for maximum dynamic pressures of approximately 500 to 700 psi. 

In regions of attached turbulent boundary layer, Figure 3, the overall fluctua- 
ting pressure level (OAFPL) is approximately 142.0 dB and for the separated flow 
region the OAFPL is approximately 152.0 dB. The separated flow region contains 
considerably more low frequency energy than the turbulent boundary layer region. 
In Figure 4, the OAFPL is approximately 152.0 dB for the separated flow/oscil- 
latinx shock interaction region; the total energy is the same as that for a 
typical separated flow region but the frequency distributions of the energy J* 
shifted an order of magnitude lower in frequency. For the region of oscillati g 
shock, the frequency distribution of the energy is narrow band in nature and 
centered in the low frequency range of the spectrum with an OAFPL of 169.0 dB. 
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The detailed environmental definition will be established after the Space 
Shuttle configuration has been defined. Due to the asymmetric nature of the 
Space Shuttle, the flow visualization and model test phases will have to be 
much more extensive and well defined than usual. During the model test phase 
it is planned to obtain the spectral and spatial correlation characteristics 
of the fluctuating pressure field. Again, complete definition of the spatial 
correlation properties, the cross power spectral density in particular, is 
required in order to realistically assess the structural integrity of the 
Space Shuttle system. 
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Plume Induced Flo* Separation 

During all of the past flights of the Saturn V vehicle, an aerodynamic 
phenomena has occurred about which very little information is known. This 
phenomena is known as the high altitude plume induced flow separation pheno- 
mena and it is illustrated in Figure 5. This figure was obtained from a 
single frame of a movie film which was taken by a camera that was mounted in 
an Air Force "chase plane" that recorded the sequence of events during the 
second flight of the Saturn V vehicle, AS-502. The ti:r.e at which this frame 
was taken was approximately 133 seconds into the flight. It can be seen from 
this figure that at the leeward siue (the left side), the visible portion of 
the flame front is at the interstage region of the vehicle; between the S-IC 
boost stage and the S-1I second stage. The flight and atmospheric conditions 
at this time is shown in Figure 6. The flight Mach number was 5.44 with an 
angle of attack of 2 degrees and a dynamic pressure of 62.6 PSF. Vehicle 
altitude at this time was 150,000 feet with the corresponding ambient pressure 
and density of 2.7 PSF and 0.00011 pounds per cubic foot respectively. 

The start of the movement of the flame front forward on the vehicle occurs 
early during the flight of the Saturn V, typically around 90 seconds into the 
flight. It appears, that as the vehicle travels along its prescribed trajec- 
tory, the rocket exhaust plume interacts with the surrounding air stream and 
thereby creates a large region of separated flow adjacent to the vehicle. T is 
causes the exhaust products of the rocket engines to be recirculated into the 
separated flow region. As the vehicle gains altitude, the expanding exhaust 
plume causes an additional enlargement of the region of separated flow and 
consequently, there is a very rapid and sometimes violent forward movement ol 
the flow separation front on the vehicle. Figure 7 depicts the movement o 
this flow separation front as obtained from analysis of the chase plane movies 
taken of the flight of AS-502. Near the end of the boost phase of the Saturn 
V flight, typically 150 seconds, the separated flow region completely engulfs 
the S-IC stage and has traveled a total distance of approximately 125 feet. 

If a similar movement occurs on the Space Shuttle, approximately 60% of the 
Shuttle would be immersed in the separated flow region. 

In an effort to acquire more information about this phenomena, an inflight meas- 
urement program was implemented on AS-505. This measurement program consisted 
of ten dynamic pressure transducers which were flush mounted in the intertank 
region of the S-IC stage. Six transducers were located in a ring around the 
vehicle and two each were aligned in a longitudinal direction along Position I 
and III of the vehicle. The instantaneous fluctuating pressure time history 
for two of the measurements are shown in Figure 8. The frequency response of 
the telemetry system used for the measurements at Station 757 and 837 were from 
0-84 hertz and 0-14.0 hertz respectively. The fluctuating pressure levels 
that were obtained were unexpectedly high. As the separated flow front moved 
past each transducer the pressure time histories were initially cl 1 ??®* 1 ® nd 
then gradually decayed as the front continued to move forward. The fluctuating 
pressure levels at which clipping occurred were typically on the order of 155 
to 158 dB with the predominate frequency on the order of - to 3 Hertz. There 
wer» manv problems encountered in attempting to obtain measurements in this 
type' of flow region, the major ones being those of the radiant energy from the 
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flow, temperature and vibration. Transient light energy from the flow causes 
extraneous signal output of the transducers that were used, therefore, a thin 
piece of aluminum foil had to be placed between the diaphragm and the protec- 
ting cover cap. The temperature caused a bias to be produced in the data. 

Efforts to compensate for this temperature bias were made prior to launch, but 
the effects of this bias can still he seen in the data. The transducer used 
had a vibration compensation network that served to minimize the efiects of the 
vibration. All efforts have been made to insure the accuracy of the data. 
Admittedly the data does seem u* real 1st ical ly high, however, no justifiable 
reasons have been found to invalidate the results shown in Figure 8. 

The following points should he careluily considered if these levels are indeed 
accurate; (1) what effect does this type of environment have on the aerodynamic 
surface and especially the thermal protec Jon system (TPS) oi the Shuttle?; 

(2) to what extent will this asymmetric flow field load the Shuttle?; (3) how 
does this loading affect the control system?; (4) will the exhaust flow products 
contaminate the orb iter rocket engines and the orb iter surface itself? 

A study program is underway at MSFC in an attempt to gain a better understanding 
of the exact nature of this phenomena. From this study an experimental program 
will be designed and implemented in order to obtain fluctuating pressure data 
and basic flow characteristics under this type oi simulated flow conditions. 




Concluding Remarks 


Efforts arc underway at MSFC to determine the preliminary environmental esti- 
mates in the other acoustic problem areas associated with the Space Shut t le , 
i.c., ground plane environment in near, mid and far- field during holddown and 
launch, blast, etc., that was discussed earlier. The results oi these studies 
will bo reported when they become available. 


With regard to the two area* that were briefly discussed, 
are given: 


the following remarks 


a. Rocket Noise Environment 

(1) Preliminary estimates oi the acoustic environment developed by 
the propulsion system of the Space Shuttle have been determined. These esti- 
mates are state-of-the-art estimates and arc adequate for preliminary design 
requirements. 

(2) A model test program is presently being designed to obtain the 
spectral and spatial correlation characteristics of the radiated pressure field 
associated with the Space Shuttle. A preliminary experimental program is 
planned to determine the effect of the increased chamber pressure as well as 
the non -conventional cluster arrangement of the engines and deflector effects. 


b. Inflight Fluctuating Pressure Environment 

(1) Preliminary estimates of the typical fluctuating pressure environ- 
ments in selected aerodynamic flow regions expected during Space Shuttle t light 
have been determined. Again, these preliminary estimates are adequate ior the 
preliminary design requirements of the Shuttle. 

(2) Flow visualization tests are presently being planned. The results 
of these tests will be used to design a scale model wind tunnel test program 

ot obtain detailed fluctuating pressure information for all Space Shuttle flight 
regimes. 

(!1) An analytical study is presently in progress with regard to the 
plume induced flow separation phenomena. The results of this study will he 
used to design an experimental test program to more clearly and accurately define 
the magnitude of the unsteady pressure and basic flow characteristics associated 
with this phenomena. 
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Sl‘MM A RY 

the dvuam i r i rspon sc s of various proposed desi gn s of 
an ant enna-f ecdhorn roniMtut ion for a communications satellite 
to random excitation in i he launch environment were determined 
using beam and plate f lint r elements idealisations, in both 
normal mode and impedance solutions, and applying power spectral 
density techniques. 

Statistical results in the form of mean square values, 
exceedances of load levels and cumulative damage formed the 
design criteria . 

The same digital programs can be used to determine space 
shuttle responses for various dynamic loading considerations, 
particularly during launch and staging. Sets of compatible 
elements which could be used in the finite element idealisation 
of the complex space shuttle structure are presented. 




I .0 1 vrKOtH’CTION 

Ihis paper deals with methods of analysis and criteria 
for the design of a flexible structure in a random loading 
environment. The particular problem of concern is the design 
of an antenna for a communications satellite, considering t lie 
random character of the launch vibration environment. 

A sketch of the antenna mounted on the spacecraft and 
enclosed by the payload envelope is shown in Figure 1. Sig- 
nificant points in the spacecraft as well as the antenna are 
shown. The random excitation during the launch phase is 
assumed to be applied in the form of accelerations in all 
directions it the base of the spacecraft. That is to say, it 
is applied at the point of attachment to the launch vehicle. 




2.0 ANTENNA DESIGNS 

Various antenna desists were proposed, using different 
materials and different basic configurations. A typical de- 
sign is shown in Figure 2. Here the antenna dish is curved, 
and the feedho n is connected directly to it, as well as 
being cantilevered from its base. Another design was such 
that the feedhorn needed only to be connected to the ase, 
at the despin motor assembly. 

Three different configurations were considered for 
the antenna dish itself. 

a) The first was a formed thin plate attached to a grid- 
work of beams. Materials considered were aluminum, 
beryllium, titanium, and carbon. 

b) The second design used a honeycomb sandwich material. 

The formed sandwich dish was mounted on a central 
post. Proposed materials for the honeycomb and face 
sheets were aluminum, and titanium. The central post 
had to be quite stiff and a carbon composite material 
was chosen. The feedhorn was also of the same material. 

c) The third design was a stretchable alloy wire mesh 

supported by a tubular space frame. This frame could 
be made from aluminum, beryllium, titanium or carbon. 
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ANALYTICAL Ml 1H<>D 


3.0 

3 . 1 St ructural Rep fps»?ntal ion 

In order to perform a dynamic analysis of a 
structure such as that of the ant enna-feenhcrn com- 
bination, the first requirement is a set of stiffness 
or influence coefficients, and a correspond in*? set 
of masses. The method chosen yields the stiffness co- 
efficients by simulating the structure by an asseml lv 
of finite elements, using the direct stiffness method. 
For preliminary design, where a great many models were 
run, lumped masses were used, because of the substan- 
tially fewer number of operations required. However, 
for final design, consistent masses were used, with 
added lumped masses at certain discrete node points. 

The consistent masses allow for a much more accurate 
representation of the dynamic properties of a structure. 

lu-sci ipt ions of work previously carried out by Canadair 
in finite element analysis are contained in References 
“ arid 


finite Elem ent Idealizations 

ar d pi a** element* **ed in the finite element 
a* ior<6 ' * e mm antenna designs are 

5 ) OWI If I n: IT* • '».(** * 1 »-!■»• i * s * 1 * decrees of 
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I nr f he 1 3 ! S’ * r f 
on a i i d * « » n k . 1 1 


a f r 
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the gridwork, and the triangular plate elements were 
used to represent the curved dish. For the second 
configuration, flat triangular shell elements were 
used to represent the formed sandwich dish, and the 
beam elements the central post. 

For the mesh designs, the beam elements were used to 
represent the back-up truss structure, as shown in 
F igure 4 . 

The feedhorn was represented by a series of beams, as 
shown in Figure 5 , for all the designs. 

Since the vibration of the spacecraft itself substan- 
tially affects the vibration of the antenna structure, 
the spacecraft was represented by a series of tapered 
beams, as shown in Figure 6, the representation being 
obviously coarser than that of the main structure of 
interest . 

3 .3 Solution Procedure 
3 . 3.1 Reduction 

The only effective limitation on the size of the original 
stiffness matrix and hence the number of nodes in the 
representation is the bandwidth. However, the eigen- 
value routines available to determine the natural fre- 
quencies and mode shapes of the structure are not 
efficient* for large order matrices, and impose a limit 
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of somewhere near the order 100. Also, if the imped- 
ance method of solution is to \ e used, only small 
order matrices may he solved efficiently, timewise. 

Hence, the original matrix si/e must he reduced by 
retaining only a certain number of degrees of freedom 
and eliminating the rest, by application of the principle 
of virtual v^ork. 

The equations used for the reduction are shown in 
Table 1. The reduction of the stiffness matrix is 
carried out on-line, as given by equation (1). No 
matrix inversion is required. 

If lumped masses are used, the mass of the structure is 
lumped at the retained decrees of freedom. 

If consistent masses are used, a simultaneous or-line 
reduct ion of the mass matrix is performed according to 
equation (2). 

3.3-2 Eigenvalue Solution 

The natural frequencies and mode shapes of the structure 
for the reduced, number of degrees of freedom are then 
found using the eigenvalue routine. 

3.3.3 Transfer Function Solution 

The vibration level of the structure over che entire 
frequency range is found by first determining the res- 
ponse of the retained degrees cf freedom to unit sinus- 
oidal accelerations at the base of * h i* sp^cpr i nf t for 
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the natural frequencies, as well as for intermediate 
frequencies, using either the impedance or the normal 
mode method of solution. 

Details of the two methods are as follows, while more 
complete discussion may be found in reference 1. 

3 . 3 . 3 . I Impedance Method 


The cquat ions 

are as ^hown in Table 2. Equation (J) 

i s the 

matrix 

formulation for the equations of motion 

of the 

system 

, with base acceleration excitation, where 

O] 

* 

mass matrix, either lumped or consistent 

w 

= 

stiffness matrix 

{*} 

- 

accelerations 

« 


displacements 



structural damping 


- 

base acceleration, zero in directions 
where there is no excitation. 


The complex impedance matrix t z 3 is f° rmec * as * n 
equation (4), for a sinusoidal excitation of fre- 
quency Itf . 

The equations of motion thus become equation (5) for 
a unit base acceleration excitation. 

The solution for the displacements and/or accelerations 
is straightforward inversion of the complex £ z J matrix 
and multiplication by the right hand side as shown by 
equations (6) and 17) • 


IMPEDANCE METHOD 

[m]!x! + (l+ig)[K]!xj = -[l 

[z] = -W 2 [m] +(l+ ig)[K] 
[Z]\x\ = -W 2 [M]il! 

?x<= — w 4 [z] _1 

or ;*;= - w 2 [z] — 1 [m] ji * 

TABLE 2 


3 


3- 3. 3. 2 Normal Mode Method 

The equations are as shown in Table 3- 


The basic equations of motion for the system are the 
same as before, represented by equation (3). 


A transformation to generalized coordinates W • 
given bv equation 1$), is applied using the matrix of 
mode shapes [*] obtained from the eigenvalue solution. 


The equations of notion are then premultiplied by 

J ^ , The diagonal generalized mass matrix ^ M^ 
and the generalized stiffness £ KjJ are defined b> 
equations l?) and 110) where the are the natural 

frequencies. 


The left hand side being now diagonal, the equations 
of motion are decoupled, and the generalized coordinates 
q* can be solved independently using equation (12). 


Equation ($) is then used to obtain the actual dis- 
placements 

9 

This is a much more satisfactory method for the solution 
of large order dynamic problems. It allows decoupling 
of the equations of motion, so that no inversion of 
large complex matrices is required, and also allows 
reduction of the number of degrees of freedom. 
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NORMAL MODE METHOD 


[M]|xi + (l + ig) [k]{x!= -[m]!S| 
ix! = [0]{q| 



[Mr] = M T [M]W 

[Kr] = [w 2 n] [Mr] 

[-W 2 [Mr] ♦ (1 + ig) [W 2 n] [Mr]]iqi = 

-W 2 [g] T [M]{1) 

!Ql “-W 2 [M r] + (1 + ig) [W 2 n] [Mr] 

TABLE 3 


(9) 

( 10 ) 


-W 2 M T [M] (1) (11) 
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Discussion of Methods 


The preliminary designs were examined on a remote 


access terminal connected to a high speed IBM 3 A 0 — / 5 


computer. The main problem was one of storage, as 


only 100K were available to the user, although tem- 


porary storage disks could also be used. The impedance 


method was deemed to be more practical for the solution 


of these problems, and reductions were made to the 


original stiffnesses to allow final solutions of approx- 


imately 16-20 decrees of freedom. For final design. 


solutions were run on Canadair’s own IBM 360—6 5 


computer with 36 5K working memory and vast tape and disk 


capability. For larger order problems (i.e. for a larger 
number of retained degrees .of freedom), the normal mode 


method was and should be used, so that there is no need 


to invert large order complex matrices. The impedance 


method may be used as a check on the normal mode solutions 


Sinusoidal launch vibration specifications are then 
satisfied by merely scaling up or down the transfer 


functions found for one g excitation- 


Response to Random Excitation Usins Power 
Spectral Density Methods — 


Once the response (i.e. transfer function) to a unit 
sinusoidal excitation has been computed, the response 


to the random excitation of the launch environment is 
computed usin^ power spectral density (PhD) methods. 


* VA'' oirr 






1 


The 


basic equations ate shown in Table 4- 


Equation (14) relates the random character of the 
response, as a function of frequency, to that of the 
excitation through the modulus of the transfer function 
where 

& - power spectral density of response 

-*our 


\tp\ * transfer function modulus 

* - power spectral density of excitation 

assuming that the excitation behaves as a stationary 
Gaussian random process. 


The root mean square value A of any response, as well 
as the number of zero crossings N Q can then be found 
by numerical integration of equations (15) and (16), 
respective lv . 

The exceedances N of any acceleration or load level y 
is then obtained from equation (17), commonly known 
as Rice’s equation. 


Descriptions of related work on power spectral density 
methods carried out by Canadair are contained in 
References 3 and 4* 
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The response PSD of the lateral acceleration ol t ».*■ 
tip, for one antenna model, is as shown in Figure " . for a 
Thor Delta launch vehicle. Natural frequencies can he seer 
at IS, 39 , 107 cps, etc. 


From this PSD, the root mean square level of the 
acceleration was found to be 7-58 g, and the number of zero 
crossings 325-4 per second. 

The plot of the number of exceedances of each g level 
is shown in Figure 9* 
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RANDOM VIBRATION LEVEL OF SPACECRAFT 
AT SEPARATION PLANE 




THOR DELTA 
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PLOT OF POWER SPECTRAL DENSITY OF ANTENNA TIP LATERAL 
ACCELERATION IN RESPONSE TO RANDOM EXCITATION 
AT BASE OF SPACECRAFT 
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Two different methods were followed in applying these 
results to fatigue, one for preliminary design, and one for 
final design. lie two methods both used the theory of Cumulative 
Damage as a final point in assessing the fatigue life, hut 
arrived at the result in slightly different ways. 

lor the preliminary design, stresses were eomputed for 
a 10 g loading at the tip of the antenna. The exceedance curve 
tor g’s was entered at successive g loadings to give the number 
of t imes these g’s were exceeded per second. These numbers 
were then multiplied hy the time in the environment, to give 
r ],e tot al number of exceedances of that level. This was re- 
p. ii.-d for increasing g levels until the total number of exceed 
. . v tor a level was less than one. The static solution for 
. « e K'c load was then factored to these discrete g loadings 

o , , s p,-., t i um of fatigue stresses obtained, with rectangular 
jr , ; « ri -jugMlai load patterns assumed. The number of stresses 
it • ! - v at 1011s st toss levels was given hy the number of exceed- 
M,,- u rwlfr.1t ion levels. 1 be cumulative fatigue 
da-rag*- vis tie! ,ynrati-ri ft<vn equation 11*1. using a scatter 
I a* f or o I 4 . * . 

« N, 

where n ^ r. number o| < a' 1 * - v e 1 ) 

V i = numl er of < 1 * s ' ■ ■ 1 a 1 1 - n • a * 1 * ' • 1 i 

Conventional wisdom sa\ s the total -••mtrai 1.0 ^ » 

less than one in order that failuie dues no* >. .» 


323 


For final design, a more sophisticated procedure was 
required. Once the displacements ^ are computed at the 
retained degrees of freedom hy either t .e normal mode or im- 
pedance method, the displacements Sg at the eliminated 
decrees of freedom can he found from the oricinal stiffness 
matrix, which has been stored on tape, by usinc Fquat i 'u (1">, 
where K 2 i and K 22 are elements of tlie oricinal part it ned 
st iffness mat rix . 

C « ik k 21 (19) 

® K2 2 

Once the displacements for all the degrees of freedom 
are found, the stresses can be determined on each element by 
cring t ht? element subroutines. These stresses now be 
come the transfer functions per unit base acceleration. When 
multiplied by the input PSD and the numerical integrations 
performed as before, the rms stresses, zero crossings and the 
exceedances of various stress levels per unit time are obtained. 
The stress exceedance curve is entered for a series of success- 
ively increasing stress levels until the exceedance for the 
time in the launch environment is less than one. The cumulative 
damage equation is then formed as before, using these stress 
levels. Thus, an estimate of the fatigue life of each element 
in the structure is obtained. 
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lor elements where there is effectively no interaction 
between stresses, the basic exceedance equation is applicable, 
lor elements where combined stresses are significant, the 
exceedances of a series of factored structural interaction 
boundaries are found, using the method of Houbolt, as given 
in Reference 2 . Ihe number of cycles to failure snould then 
be that value as specified for the combined stress level. 
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These methods have proven must satisfactory in the 
design of t he w t <*nna -f eedhorn combination. Ihe same pro- 
■trams used to analyse the antenna response during launch 
can equally veil he used to analyse dynamic structural vib- 
rations in the space shuttle, paitiiularly during launch 
and st ail in a* However, a much larger’ choice of elements is 
necessary than were required for the antenna, in order' that 
the complex shuttle structure he represented accurately. 

A large number of elements have been thoroughly 
developed and tested under research grants from the Canadian 
Defence Research Board. 

Two different sets of compatible elements that could 
be used in idealizing the shuttle structure have been assembled 
and are presented here. 

The first group, shown in Figure 10. have corner nodes 
only. They allow linearlv varying displacements within the 
element, and hence have constant strain. 

The second group of elements, as shown in Figure 11. 
have intermediate nodes along their edges. This allows the 
displacements to vary quadrat ical ly : hence, the elements 

have a linear strain capability. The hexahedron is also iso- 
parametric. allowing the faces of the element to represent 
curved surfaces. 
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SPACE VEHICLE RESPONSE TO ATMOSPHERIC DISTURBANCES 
Robert S. Ryan 

NASA -Mars ha 1 1 Space Flight Center 
Huntsville, Alabama 

The objective of dynamics analysis is to predict the vehicle response 
(load) statistically (probability of expectancies of load) using a vehicle 
model which is known only in some bounded fashion (nominal plus spreads), 
and an input force (winds) either as an ensemble of individual winds or as 
a statistical quantity. The prediction and interpretation of the response 
of a space vehicle to atmospheric disturbances are very complicated because 
the vehicle response depends on, and interacts with, the various character- 
istics of the wind (magnitude, shear, gust), the vehicle dynamics (rigid 
body, elastic body, propellant oscillations, aerodynamic loading), the control 
system (attitude control, modal suppression, path control), and the statis- 
tical characteristics of each of these subsystems or variables. Technology 
associated with space vehicle response is, therefore, concerned with many 
facets or disciplines. These can be categorized as: (1) dynamic models 

for structure, liquid propellant, aerodynamic forces, control system, 
environment (atmospheric disturbances) , and overall combined system; (2) 
analysis techniques for frequency response, time response, and stability 
analysis, including appropriate techniques for statistical description or 
interpretation of results; (3) criteria for evaluation of results such as 
flutter boundaries, handling qualities, response goals, stability goals, 
design goals (probability of launch in worst wind month, etc.); and (4) 
methods for alleviation or suppression or excessive loads. 

Obviously, one cannot treat all of these areas in one paper. Therefore, 
only major items will be mentioned with emphasis placed on the analysis 
area. The approach taken will be to present state of the art approaches 
and to identify the basic technology development needed for the Space Shuttle 
concept. In addition, special inflight dynamic problems which are expected 
for the Space Shuttle vehicle due to its varied mission and innovated design 
will be identified. 
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SPACE VEH1C 1 F BI-SPOUSE TO ATMOSPHERIC DISTURBANCES 


The objective of dynamic analysis is to predict the vehicle response (load) 
statistically (probability of expectancies of load) usinn a vehicle mo e a 
?s tooii only in some bounded fashion (predicted nominal plus tolerances), and 
In input force (usually winds) either as an ensemble of individual winds r . 

a statistical quantity. 

The decree of accuracy with which we can make these predictions determines 

th . structural 6 integr i t y of the vehicle which can he traded directly into greater 
the structural g y accuraC y is parlayed into the control system 

2 ™ ^ - r " ult 

-h. soace Shuttle vehicle, as now envisioned, is payload sensitive, fatigue 

SSfeSSSSSH^SSK* 

developed if the Space Shuttle is to become a reality. 


SPACE VEHICLE RESPONSE TO 
ATMOSPHERIC DISTURBANCES 


PROBLEM: 


SIMULATION OF DYNAMIC LOAOS AND RESPONSE OF SPACE 
SHUTTLE VEHICLES 


OBJECTIVES: EHEOICT VEHICLE OVHAHHC RESTOHSES AND LOADS WITH AH 
ACCURATE. DEFINABLE. PROBABILITY LEVEL. 


IMPACT: 


VEHICLE STRUCTURAL INTEGRITY. PAYLOAD. AND CONTROL 
SYSTEM REQUIREMENTS 


REQUIREMENTS FOR DETERMINING DYNAMIC LOADS AND RESPONSE 

Accurate and efficient dynamic analysis is dependent upon five basic 
factors: (1) Analytical or test derived models of subsystems which are 

combined into the overall system model. This overall system model must 
include accurate descriptions of the subsystems, liquid propellant, 
structural dynamics, control system, aerodynamic forces, and flight 
mechanics trajectories; <2) An accurate description of the environment 
compatable with the analysis technique. This is as important as the model 
because it is the major excitation force of both static and dynamic lateral 
loads; (3) Analysis procedures which are efficient, lead to understanding, 
and provide accuracy of the results; (4) Performance criteria. These are 
necessary for performance goal settings, both in the design and verifi- 
cation phase; otherwise, ultra-conservative designs result; (5) Active 
load reduction techniques (mainly control system) is necessary as a final 
means of meeting design goals and mission constraints. 

REQUIREMENTS FOR DETERMINING 
DYNAMIC LOADS AND RESPONSES 

I. MODELS 

A. SUBSYSTEMS 

B. SYSTEM 

II. DESCRIPTION OF ENVIRONMENT 

HI. ANALYSIS PROCEDURES 

IV. PERFORMANCE CRITERIA 

V. LOAD REOUCING TECHNIQUES 
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FLIGHT PHASES REQUIRING LOADS CALCULATIONS 





All major areas of the flight profile from the launch pad to the return 
landing require accurate predictions. Each phase can, and usually does, design 
a particular area or structural subsystem of the vehicle. The chart llst ® 
the major areas of concern; however, only the ascent phase will be covered In 
this presentation (asterisk) because of time limitation. Many of the points 
covered will be applicable to the other mission phases. 



FLIGHT PHASES REQUIRING 
LOADS CALCULATION 
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• ON PAD 

• LIFTOFF 
H t m ASCENT 

• SEPARATION 

• DOCKING 

• REENTRY 

• CRUISE 

• LANDING 


334 




— V ** “ — : * ' fcV .. I P ’ : ~~V ’V " ' vt? I 

: ■ .0 . • ■- u t • . ‘ m - * ■ i '“i’:'" ‘ / , ‘ y « ‘ 

--T- h v t v ’ i ’ r-'-^ — ? ~-r,v'r ' • 1 .u, w w i ws f w wi wni we^ai^. ■ i rt I t H RWF N^RiiBMNM^iRfthiNaiNiB 


-iVi: * *-»**»■ . 1 - ij, 

■a tttfciatik iiaa 


HIM 



MODEL HKOUIKLMLNrS 


The system model should include, hut not necessarily be limited to, the 
listing shown on Chart 4 . Accurate models of these various subsystems allow 
a choice to be made as to the degree of sophistication needed for a particular 
analysis. In the final design phase and assessment phase, all the model details 
must be included. 


MODEL REQUIREMENTS 


60 RIGID BODY TRAJECTORY 
STATIC AEROELASTIC EFFECTS 


NON-IDEAL CONTROL SYSTEM (III ll l<\ I.AdSi 


ENVIRONMENT <ltl\l> l>l S(l(inin\i 


BENDING MOMENT OR L0A0S 


30 STRUCTURAL MOOES 
PROPELLANT SLOSHING MODES 
NON-LINEAR, LOCAL AERODYNAMIC FORCES 


PROPULSION SYSTEM 





SYSTEM MODEL .NC 

The state-of-the-art of dynamics system modeling has reached a high 
level of sophistication in the last few years. System equations for a 
six-dearee-of-freedom trajectory, using three-degree-of-freedom elastic 
bodv descriptions, have been formulated. Nonlinear quasi-steady aerodynamic 
distributions were incorporated, as well as nonideal control systems (filters 
and lags). Programed control system gains, time-varying co-efficients . 
anc some means of account for data tolerances, are available. 

Lven wjt:. these advancements, adequate models are still the major 
-TO! le- s recent publication from NASA Electronics Research Center 

Tr * n ,« JT . control Research and Technology," surveys the most 
I, r rchlems facing control engineers. Modeling was listed by a 

:: .. . ... . , e rJK ., rt * as being one of the major areas of research need- 

.... . j 0 r.r. b. lewis, Pennsylvania State Vniversltv, says, 

‘-rre'is'arsc.-tlev nc substitute for a thorough knowledge of the system. 

t ’is a teaic-s anc time-consuming "recess requiring much ingenuity to obtain 
‘ * u - s ,- s . e - acce.s or. wMct: t-.e control d.sier, can be based. Good genera 

.; s . rocecures are neeoef that even complex systems can he satisfactorily 
“ *. .efuovitt, -ase astern Reserve. «av S , W need much more 

c . Active atans cf moaeiin ar.c a-stractir . ., f r om t-e w° _ a e±JL.± 

attributes mat are relevant :c me oecision-maUng anc centre, nrobler. 


c -"i 


; b a tecics..* an 


& e * u - s v *> 
est : roce 
escrit'en . 
! iective 
rinrutfs 


-U-. -ne aooeii-.t teenne lor* v-ic- t- neeoec for ver icle response 
tc d i st urn ances is covered m other oaoers oealmg witr. the bolt 
rheir importance dictates e restatemert cf me ^ Mlc effects at 

Hire accurate elastic body cnaract-r ist ice that include .ocal e * f * c ** . 

sensor locations and m** cross-cour ling, -articular e-nhasl* is 
on mints and localised damming. Nonstationar* aerodynamics m a cractica 
fora°f or S response is a dir. n~d. as -li « w«* in defining 
the data spreads associated with the characteristics o t . t,e 
so that they can be statistically accounted for in the response analysis- 
Finally an efficient statistical proceflure is needed for ana.vilng these 
data tolerances in conjunction with the response analysis. 




SYSTEM MODEL 


STATE OF ART: 

• B-D TRAJECTORY 

• 3-0 ELASTIC BODY 

• QUASI-STEADY AERODYNAMICS 

• PR06RAMME0 CONTROL SYSTEM GAINS 

• RSS PARAMETER VARIATIONS 

TECHNOLOGY NEEDED: 

• NON-STATIONARY AERODYNAMICS 

• IMPROVED SUBSYSTEM MODELS 

• METHODS FOR ESTIMATION OF SUBSYSTEM 
DATA ACCURACIES 

• EFFICIENT STATISTICAL PROCEDURES FOR ANALYZING 
DATA TOLERANCE EFFECTS ON RESPONSES 




VEHICLE LOADS CHARACTERISTICS 



Dynamic analysis of present launch vehicles has produced some basic 
results applicable to future space vehicle design and the development of 
technology • The basic wind loads acting on a vehicle can be expressed 
as a lateral bending moment with the longitudinal loads added to these 
effects. This bending moment can, in general, be simplified to be a 
function of the angle of attack, control force deflection, and elastic 
body accelerations. Other effects are generally small for vehicles of 
the Saturn type. Three basic dynamic effects influence design of the 
vehicle structure: ground winds, inflight winds of high dynamic pressure 

recions, and first stage shutdown transient, Typical structural design 
areas for the Saturn V/Apollo are shown on figure 1. Figure 2 shows two 
verv imnortant effects: (1) That from a rigid vehicle standpoint, the 

design of the nose of the vehicle is predominantly influenced by the 
aerodynamics . This is illustrated by the solid curve R(x) which is the 
ratio of the bending moment partial due to a^gle of attack to the bending 
aoment partial due to control force R(x) » w a . (2) '"he percentage of 

M’ 

n a 

lotal r>enc;nr moment resulting from bending dynamics. Here again, the 
v« icle no** is the culprit with the major portion of the bending moment 
i- t*;* region resulting from bending dynamics. Additional effects on the 

'Apollc vehicle loads were gust penetration vhich had no appreciable 
effects anc static auroelasticl-tv , vhich increased the engine deflection by 
aorroximatel* 1 1. urn t>e large aerodvn<»wic surfaces present on the 
: - wT, tie ver.u> anc its compi** structural design, these effects can he 
lomrcunaeG fro* : * xtxjntrnced on tne Saturn vehicles. Technology 

t~.er ®ust L-e exttuuec it u»e aeroeias licit v area of local nonstationary 
aerodvnaadcs . ^narration. static ae roe last! city , and structural 

d vti am* c Mooe 1 ms. 
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VEHICLE LOAD CHARACTERISTICS 

■ENDING MOMENT MgU' * [ R(i to(D ♦ -Vt)] ♦ I M^W «• I » s fJA * % 

1. SATURN V OESISN LOADS 2 AERODYNAMICS, THRUST FORCE, 

AND RENDIN6 DYNAMICS 

STATION IN METERS STATION IN METERS 




3 STATIC AEROELASTIC EFFECTS FOR SATURN V. J INCREASED 2«\ 
«. 6UST PENETRATION EFFECTS FOR SATURN V. NEGLIGIBLE 
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ENVIRONMENT 1 


Four distinct types of wind inputs are available for the appropriate 
response analysis: (1) discrete, (2) power spectra, (3) nonstationary 

(stochastic) , and (4) individual wind soundings. The accuracy of the 
vehicle response is obviously directly proportioned to the accuracy and 
understanding of these inputs. 

Discrete winds are used mainly as synthetic profiles and l~cc~ >e 
gu«*ts. The state-of-the-art profiles are based on many individua- .ound- 
ir.gs and are available for ETR, SAMTEC, White Sands, and Wall mt. .ind 
shear is a conditional shear based on a reference level wind --pee-. 
Technology needed for these discrete profiles is. a development of joint 
statistics of the shears and wind speeds. 

Power spectra exceedance models for longitudinal, lateral, and 
vertical gust components are available on a worldwide basis. The power 
spectra vertical wave lengths of lOOn to 2,000m are available for the ETR 
test range. Additional development is needed to determine the cross 
spectrum from these same data in order to determine more a curate response 

data. 


A nonstationary wind representation is available based on the 
Rawinsonde profiles (1000m increments) for ETR, and contains interlevel 
correlations. The major development in this area is the shaping filters 
and inter-level correlations for the high frequency wind characteristics 
(100 <. X < 2,000m). 

Detail wind profiles (Jimsphere) based on 25m increments are available 
for ETR, SAMTEC , White-Sands-, and Wallops. The present samples are not 
adequate to duplicate the wind speed and wind shear statistics of the 
Rawinsonde ensetf>le. The turbulence portion is adequate. Either a larger 
wind saaq>le is needed or a statistical procedure developed to correct the 
present sample. 
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ANALYSIS PROCEDURES 


Analysis procedures fall Into the same categories as the environments 
(1) discrete inputs, ( 2 ) harmonic analysis, ( 3 ) nonstationary (stochastic), 
and ( 4 ) Monte Carlo. The analysis using discrete inputs is well known and 
covers various degrees of sophistication from simple planar point time analysis 
to full 6-D trajectories with 3-0 elastic modes and nonlinear aerodynamics. 

This analysis procedure is adequate. 

Generalized harmonic analysis is well developed for both loads and 
lifetime analysis of linear point time systems. This approach has not been 
generally extended to include cross spectrum wind models or analysis 
techniques. 

Honstacionary (stochastic) procedures are state-of-the-art for linear 
time— varying systems which include a wind filter for 1,000m wave lengths. 
Methods for determining response to mean wind and adding the results from 
the stochastic process are available. Extending this approach to nonlinear 
systems, witli the development of the high frequency wind model, would produce 
a very valuable tool in loads analysis. This approach also needs a method 
for lifetime prediction. 

The final approach, the most desirable because no basic restricting 
assumptions are necessary, is the modified Monte Carlo which flies the 
vehicle through an ensemble of individual wind profiles, computes the vehicle 
outputs, and analyzes them in a statistical manner (response level versus 
probability). (See slide 9 .) This approach is state-of-the-art for a 

6-D rigid body trajectory, nonlinear time-varying coefficients, nonideal 
control, propellant oscillations, and planar elastic modes. Present state- 
of-the-art computation time for one trajectory using a hybrid computer is 
seven seconds. Technology development involves more efficient models and 
computational procedures that produce desired accuracy while reducing 
machine run time. Additional developments include procedures for in- 
corporating data spreads (tolerances), nonlinear local angle-of-attack 
aerodynamics, and a procedure (mentioned under "winds' ) for correcting 
the wind sample. 


342 


ANALYSIS PROCEDURES 




PERFORMANCE CRITERIA 


Vehicle dynamic response analysis and response reduction techniques 
require accurate performance criteria as a base for interpellation and 
design goals. A state-of-the-art optimal approach (Honeywell s) is based 
£ a wind .hd . wind <l.000„ U.,t *“““»>• “«* 1 

state and wind sensing. with constraints on glnbal awl*. 8«>*1 
and vehicle drift. Bending moment response is usually cnosen fo p 

alzation. 

Technology is required for nonideal state estimation which includes 
sensor choice^nd location criteria, an improved optimal P^e^re to \ 
computation efficiency and greater model detail, and * v ®****® 

Dresent optimal criteria. Extension of the criteria to includ- other 
responses^ conjunction with the bending moment is needed (example, 
acceleration at cockpit), and the development of a time-varying ®nalog 
of the frozen point time criteria is highly desirable as a simplifi 
means of evaluating system response. 


PERFORMANCE CRITERIA 


STATE OF ART: 


• WIND - MEAN WINO PLUS STOCHASTIC (1Q00M ACCURACY) 

• I0EAL STATE AND WINO SENSING 

• CONSTRAINTS ON GIMBAL ANGIE, GIMBAL RATE, AND VEHICLE DRIFT 

• OPTIMIZED TO BENDING MOMENT 

• FROZEN TIME POINT CRITERIA 


TECHNOLOGY NEEDED: 

• NON-IDEAL STATE ESTIMATION 

• 2SM STOCHASTIC WINO MODEL 

• IMPROVED OPTIMAL PROCEDURE FOR COMPUTER EFFICIENCY AND 
GREATER MODEL DETAIL 

• ESTABLISH VALIDITY OF PRESENT CRITERIA ANO MODIFY TO 
CORRECT DISCREPANCIES 

• DEVELOP CRITERIA FOR BENOING MOMENT PLUS OTHER RESPONSES 
SUCH AS CREW COMFORT (ACCELERATION OF CREW STATION) 

• TIME VARYING ANALOGY OF FROZEN POINT CRITERIA 
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nKNniNC MOMENT VERSUS PROBABILITY 

Typical results for the bending moment for an ensemble of 1,000 winds 
is shown for the Saturn V/apollo type launch vehicle. The results are 
given for the total wind ensemble and for the wind ensemble with the 
turbulence filtered out. Accelerometer feedback gain, g2» is used as 
the parameter. The only restriction in this probability statement is not 
in the method of analysis but in the accuracy of the wind ensemble since 
the meth"'* is not restricted by linearity or stationarity assumption:: . 


BENDING MOMENT, STATION 90 (10 4 N m) 
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LOAD RELIEF AND MODAL SUPPRESSION 


Thia topic is covered In detail by other authors. The Inclusion of the 
chart Is merely for additional emphasis to the importance of reduction 
techniques and the emphasis It places on other technology requirements such 
as model accuracy* 

The basic technology areas are passive reduction techniques (mean wind 
biasing), active reduction techniques, performance cr ^ er * a » atate ® s *V^* lon 
and separation, and efficient analysis techniques. All ot these should 
pursued simultaneously* 


LOAD RELIEF. AMD MODAL SUPRESSION 


STATE OF ART: 

• PROGRAMMED GAINS 

• SENSOR CHOICE: ACCELEROMETERS, RATE GYROS. POSITION GYROS 

• MONTHLY MEAN WINO TRAJECTORY BIASING (ALL PLANES) 

• MIXED STATE ESTIMATION 

TECHNOLOGY NEEDED: 

• ADAPTIVE GAIN SCHEMES 

• PREFLIGHT WIND BIASIN6 SCHEMES 

• INFLIGHT WIND SENSING AND WIND BIASING 

• TECHNIQUES FOR DESIGNING PRACTICAL OPTIMAL SUBSYSTEM 
CONTROLLER USING OPTIMAL PERFORMANCE CRITERIA AS GOAL 

• SEPARATE (MODES) STATE ESTIMATION 

• TECHNIQUE FOR MINIMUM INTERFERENCE (COUPLING) THROUGH 
CONTROL SYSTEM 

• SENSOR CHOICE AND LOCATION CRITERIA 

• MORE EFFICIENT ITERATION PROCEDURES 


346 





■— V 


SPECIA L SHUTTL E LO ADS PROBLEMS 


The final chart lists special Shuttle load problems that are potential 
technology areas. Whether state-of-the-art techniques can he used to 
analyze and solve these problems will he determined only through analysis. 
High priority must be given to these areas early to determine their design 
problem potential and, hence, determine if additional technology is needed. 


SPECIAL SHUTTLE LOADS PROBLEMS 

(POTENTIAL TECHH0L06T AREAS) 

1. AERODYNAMICALLY STABLE LAUNCH VEHICLE FORCES TRADES 
BETWEEN PERFORMANCE. DRIFT, AND LOADS. 

2. MASS AND AERODYNAMIC BIAS FORCES CREATE A STATIC OR 
TRIM LOAD DURING ASCENT. 

3. MASS OFFSETS AND DYNAMIC MASS COUPLING CREATE DYNAMIC 
LOADS THROUGH LATERAL LONGITUDINAL COUPLING. 

4. AEROELASTIC PROBLEMS DURING ASCENT MAGNIFIED OVER 
PREVIOUS SPACE VEHICLES DUE TO LARGE LIFTING SURFACES. 

B. LIFTOFF TWANG LOADS ARE MAGNIFIED DUE TO MASS OFFSETS 
ANO LARGE AEROOYNAMIC SURFACES. 

8. ASCENT LOADS ARE A TRADE OFF BETWEEN BODY LOADING 
ANO WING LOADING. 
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CONCLUSIONS 


The critical technology areas for Space Shuttle vehicle response to 
atmosphuric disturbances are presented on the chart. Other less critical 
but highlv desirable ateas were identified in each basic area * ,,an Y ® 
the areas' shown fall into the various disciplines of Space vehicle design 
and should receive more in-depth discussions by specialists in these areas 
however, they are of utmost importance to response analysis. 


CONCLUSIONS 


CRITICAL TECHNOLOGY 

• PERFORMANCE CRITERIA 

• NONSTATIONARY AERODYNAMICS 

• ADAPTIVE LOAD RELIEF SCHEMES 

• CRITERIA FOR CONTROL SENSOR 
CHOICE AND LOCATION 

• MORE ACCURATE SYSTEM MODELS AND 
FASTER COMPUTATIONAL TECHNIQUES 

• MOOAL OR SUBSYSTEM DECOUPLING 
TECHNIQUES (HARDWARE ANO CONTROL) 

• TECHNIQUES FOR ACCURATELY INCORPORATING 
DATA TOLERANCES IN RESPONSES. 
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ASSESSMENT OF THE POTENTIAL FOR LOAD ALLEVIATION 
CONCEPTS FOR SPACE SHUTTLE 1T EH1CLES 

R. H. Lassea and J. H. Wykes 

North American Rockwell 

INTRODUCTION 


The space shuttle vehicle (SSV) concept utilizes reusable booster and 
01 otter vehicles to place a variety of payloads into earth orbits for an 
ordor-of -magnitude lower cost than that of previous launch-vehicle/ space- 
craft combinations. The structural weight of reusable booster and orbiter 
vehicles is a key item in determining the eventual payload and consequent 
economic feasibility of the SSV. Thus concepts that might result in saving 
vehicle structural weight must be given serious consideration as early as 
possible in the design process. 

The mated booster-orbiter configuration functions as a combined launch 
vehicle and aircraft and experiences the structural dynamic response problems 
of each. The many aerodynamic surfaces of proposed SSV configurations will 
result in response to atmospheric turbulence over a broad frequency band. 
Additional excitation of the mated vehicles, as well as each craft separately, 
will be induced by separated flow or buffeting, interference effects, and 
interaction of shock waves with the boundary layer. The possibility of high 
stress, low cycle fatigue must be considered due to the multiple reuse 
requirement of this program. High loading conditions will occur during 
boost flight, in maneuvers, and at entry-cruise transition, and the aero- 
dynamic surfaces may be flutter-prone either during the transonic-high dynamic 
pressure flight regime or in transversing the stall region following high 
angle of uttuck entry. One way of coping with many of these problems is to 
strengthen and/or stiffen the basic load carrying structure. The effect 
of the resulting weight increase, however, could so seriously Jeopardize 
the puyload capability that the SSV mission could not be accomplished 

During the pust decade or so, automatic control systems have been designed 
to alleviate induced loads and/or control structural mode dynamics of flexible 
vehicles. Many studies have been made and several flight programs have 
demonstrated a variety of benefits and advantages ior both aircraft and launch 
vehicles. The Duplication of such systems and concepts to the aforementioned 
shuttle vehicle problem areas is suggested as one potential means of meeting 
the very btrtngont structural mass fraction requirements of the SSV. This 
paper briefly reviews the state-of-the-art in load alleviation and mode 
control, discusses some of the recent techniques, and assesses the potential 
of such systems for the SSV. Finally, some of the problems requiring solution, 
and the supporting programs t.o supply these solutions, are outlined. 
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w * P ALLEVIATION AND STRUCTURAL MODE CONTROL 
som RELENT STl DI ES AND EXPERIMENTAL PROGRAMS 

Becau^r of Rission requirements ana basic physical differences, the approaches 
u tht studv of load relief and structural nodi control systems havo been 
flifferem for Uun<h \ eh u 1 1 s and aircraft, »e haw arbitrarily selected 
several a irv raf i -t-r lertrcj programs which h*u proi ided both analytical and 
flight test validiit iur. _f important 1* ad relief and structural mode control 
concepts, These types ;f svstewis hav * direct application to the SS\ conf igura- 
tions because of the- use * * * aeroov r»air j t surfaces. 

One analytical study initiated in 1«*M for the Air Force Might Pyna*«cf 
Laboratory usee the XP-Ti to studv iicvibli vehicle- effects. It was demon- 
strated that accelerations and structural loads du* to structural action 
can be alleviated by a control systen that increase!* the damping ratio ol the 
structural modes. The effects on sensing, control forct application, coupling 
and adaptability to changing conditions were studied. It t?=a shown that a 
relatively simple, practical, conventional control system can bt designed to 
provide phase and gain stabilization of the structural modes under varying 
flight conditions. This system was called 1LAF for Identical location of 
Accelerometer and Force. (Fig. 1) 

Another Air Force-sponsored program (15*65), best known by the name Gust 
Alleviation and Structural Dynamic Stability Augmentation System <GASDSAS), 
used analytical and fixed-base simulation techniques to explore gust allevia- 
tion and structural mode control systems, ride quality effectiveness and 
interface problems with handling qualities, displays and terrain following. 
Building on the results of the XB-70 study, practical solutions were achieved 
for problems of force generation, system stability, and interference with the 
regular stability augmentation. Theoretical applications were made to the 
Advanced Manned Strategic Aircraft (AMSA), the C-5A, and F-lii vehicles. 
Practical applications of modern optimal control theories were alto developed 
during this study. 

Under a jointly sponsored NASA/ Air Force program, a comb in- d analytical- 
flight test project to verify the ILAF concept on the XB-70 was undertaken 
in 1S69. An ILAF system was designed as an add-on to the airplane's stability 
augmentation system to control only the symmetric structural modes using the 
existing elevons as force generators. Analytical evaluations showed that the 
system would be stable and effective in controlling the first three structural 
modes without degrading handling qualities. Pilot comments and flight test 
records verified the fact that the system improved crew ride quality in 
atmospheric turbulence. 

Another important analytical-flight test program was conducted for the 
Air Force Flight Dynamics Laboratory using the B-52 aircraft. This program 
Is best known by the name Load Alleviation and Mode Stabilization (LAMS) 
and was conducted to demonstrate the capabilities of an advanced flight 
control system to alleviate gust loads and control structural modes on a 
large flexible aircraft using existing aerodynamic control surfaces as force 
producers. A direct lift system was also studied. Flight demonstration of the 
LAMS flight control system was conducted to provide a comparison of analytical 
and experimental data. The results obtained showed that significant reduction 
in fatigue damage rates were obtained* (Fig. 1) 
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TYPICAL SENSOR & FORCE GENERATOR LOCATIONS 





SYSTEM ADVANTAGES 

The chart shows the major benefits to be obtained from using various load 
alleviation and/or mode control schemes. Obviously, not all systems can 
provide all these advantages, nor are they designed to do so. However, it 
is important to note that when fly-by-wire concepts are employed, as on the 
SSV, implementation of multiple functions is much easier. Early participation 
of the control system designer and the structural analyst is especially 
important in the requirement for force generators where use of already 
existing control surfaces generally compromises the system effectiveness. 

Implementation of an automatic control system to reduce gust and buffet 
response can reduce peak accelerations, loads, and bending moments with 
potential saving in weight due to reduced fatigue damage rate. Improvement 
in ride quality is also possible. Systems to control gust responses on 
flexible vehicles using basic stability augmentation, direct lift and 
structural mode control techniques have the largest technology background of 
analysis and test. Little has actually been done in stuoying the control 
of buffet response; in principle structural mode control techniques should 
be applicable. The potential savings in weight to a given vehicle design 
due to *ust alleviation and buffet control can only be assessed after all 
factors impacting on strength and stiffness have been determined. A number 
of gust load and fatigue analyses of large flexible aircraft have shown that, 
in general, these vehicles have not been fatigue limited to a great extent. 

This means that- structural material provided for discrete loads (gust and 
maneuver) and for stiffness (stability, control, and flutter) have provided 
more than enough material to cover most fatigue requirements. 

Obtaining required flutter margins through use of automatic control of 

1. b.co.™ l»cr...li. ( ly ..tr.clv.. This Is p.r.lcul.rly 
SO since fly-by-wire technology is developing the redundancy techniques and 
actuation hardware required to obtain reliability equivalent to basic struc- 
ture. Thus, weight saving may be possible through reduction of a lifting 
surface's basic stiffness. Next to gust response control, this concept has 
developed a considerable analytical and technology background. The ‘TP* ° f 
flutter to be controlled, lightly damped or explosive, will affect the details 
of the approach to be utilized* 

With the advent of fly-by-wire technology, the possibility of saying weight 
and reducing drag by reducing the horizontal and vertical * 8il ® “ 8f 
thus reducing static stability has received considerable attention. The 
basic stability required would be achieved using automatic controls. A 
balance between the reduced tail size and control surface size ® nd d * f 
limits set by takeroff rotation, maneuver, and stability augmentation requ* 

■ents »ust be achieved. 

By proper attention to the location of control forces required to maneuver 
a vehicle the aerodynamic loadings can be shaped to reduce design bending 
moments and thus save weight. Combined direct lift concepts together with 
structural mode stabilization concepts show attractive potential for imple 
menting this idea. 
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KEY DESIGN CONSIDERATIONS 

ss = ~ss ;:i=crr,-E-;: 

,ess is dependent upon the control laws utilized. 

S&^JEiSa&S&Z&t 

ss-s™'— ■ i :;s" — -• 

found gust sensing vanes and pressure sensors effective. 

ailerons, and elevator^ Leu g ^ K ^ be#a ullU /,.d. Often times 

attention. Up controls and sp< t , lf the basic 

these controls have been designed to do do f ^ *“*£ , IS pf>rt of the 

^ ^j R r ,^sr 8 Tri. 

for Structural mode control con b. effi- »nt and efUctlvc. 

i i i i ri i> n u s t utile? ti ifxtrnh tvi. ly 

On launch vehicles and miss..cs, -njine m • 1 . nc ^ ^ ^ ^ 

to obtain thrust vectoring nr <n that reaction 

Jet. nl»" nfnil'l b. ilftK l lv- (.hi ' , — coni ttiuoim il-nmnd i.nMn, 

sr^r rr^r: iz'r,^ ~ - - 

actuating system needed for aerodynamic ourfaees. 

over the oast decade considerable research and development effort haw *»-en 

ESS 

Mint and taking advantage of normal actuation and aerodynamic Ug«. with 
^ assist f^ electronic compensation, a UM UJdl- 

- rr" 

is representative of a more complex but ^r^litegrated^-bt Uty 

uai,,g rs <*— 

control "theories were effective In syntheal zi ng these types of systems. 
^rreer'uUnS^^ru^^^hrthru;! vector control 

. 1 .. ...... thrust accelerations make acceleration alignment a problem. 
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ASSESSMENT OP LOAD ALLEVIATION AND STRUCTURAL MODE CONTROL 

POTENTIAL POR SSV 

Having reviewed the current atate-of-the-art, examined aome atudiea and 
m«h? program ro.olta, ... «... laport.nt h.rt.ar. toaiorra, 

thia auction examinee the application of all thia technology to the SSV. 

By relying on automatic controla to provide a satisfactory equivalent static 
stability it may be possible to reduce the horizontal and vertical tail 

a .“log. lo ..lgl.1. The po.eotl.l lor Oolo, thl. t. rated (air 
to good. The lark ol ourr.ot lo-d.ptb .tody ol th. aoooot of t.11 
vnHur t inn« x>s Bible contribute* to the cautlou* comment. In addition, the 
current lack of ability to produce large capacity 

unit, adds further caution. The actuation unit is one of the m. " 
to Keep in mind when making compariflons between Apollo and SSV fly by wi 
av a terns. It la, o' course, the technology developed in the redundancy an, 
reliability of the electronic components of fly-by-wire systems like Apollo 
that suggests some optimism in evaluating this approach. 

The use of an active flutter control system on the SSV to reduce rigidity 
requirements has a good chance for success. All of the UftinK ^facx-s 
are expected to have their minimum flutter margins during boost while in the 
transonic flight region. During reentry, stall flutter of theao 
surfaces may also be of concern even though the dynamic pressure will be 
much reduced. Leading edge, trailing edge, and tip controls as well as 
reaction jeta are possible force generation sources. The control would provide 
for margin only. While the vehicles would be vulnerable for n«U«f 
flight, during boost the exposure time is relatively short. The potent a 
for controlling stall flutter also exists and the overall aasoHsmcnt is good . 

Maneuver load control for the booster during entry snd turn and for orblter 
during transition is rsted of fair to good potential. The qualification shad- 
ing of fair is related to the fact that sufficient analyses have not been 
accomplished to determine if a weight penalty exists due to these mane, -vers. 

If a design penalty exists, surface locations snd the low frequency surface 
actuation demands make the chance good for potential weight savings. 

The potential for obtaining benefits from controlling re»p< ^*e to gusts and 
buffet i;o reduce bending moments and Improve ride quality and structural 
fatigue are rated fair to good. The strength designs of the lifting 
could be set by the q magnitudes during boost where the peak values of this 
parameter could be caused by gust. The maximum exposure to gusts and maximum 
dynamic pressure occur almost slmul tanoi.usly. The chances of Implementing an 
automatic control system to reduce these design loads are rated good since 
much of the *t*to-of-th<j-*rt toohnoloKy for thi* kind of applical o . 

The HHV’u are flexible snd the crew will bo subjected to some acceleration 
environment. However, the fnct that the crew will have mainly a 1 Lor 1 «gt 
function except during doTdtlng, and the fact that the exposure time will 
relatively short, indicates gust and buffeting will have llttie effect on 
crew efficiency or safety and ride quality Improvement is not critical. 

Htu r* have shown that structural response due to buffeting is mainly first 
wing bending; thus control surfaces near the tip should be able to damp this 
relatively low frequency motion If thia Is determined to he desirable. 






TO GUST/BUFET STRUCTURAL FATIGUE NOT LIKELY 





CONCLUSIONS AND KBCOMMEN DAT IONS 

The use of load alleviation and structure, mode control concepts for SSV 
should be explored In more detail. Sufficient analytical and flight tea 
experience exlBts for both aircraft and launch vehicles to give an oye'all 
good possibility of success of this approach. In addition, the existing 
f l'y-by-wire avionics technology gives confidence of being able to provide 
the required safety and reliability in this portion of the hardware system. 
However, the unique nature of the space shuttle vehicle results in complex 
nroblems In system design. Adequate verification that these problems have 
beeS solved will requirS the best analytical efforts supported by structural 
and” aerodynamic tests of specific models of the SSV. To assure halm, 
adequate Bystem will evolve, a number of research and development actlvl- 

ties must be pursued. 

To isolate the potential problems of applying load relief and structural 
mode control systems to the SSV, the primary structural design 
must be Identified. This requires In-depth analyses of all cond l lo „s 
affecting strength «nd stiffness as a development level of offoit. These 
analyses must Identify the structural response to maneuver loads, gust and 

ZftsZ and to flutter. Currently, satisfactorily accurate analytical cal- 
culations can be made for maneuver loads, gust loads, and fluUct c 
♦eristics. The analytical capability for estimating excitations turn 
burfet and dynamic interference loads Is not adequate and improved me hods 
are needed as well as experimental data from wind tunnel monel test piograms 
to confirm all calcul ationp. 

Any system design synthesis will require detailed knowledge of the vehicles 
structural dynamics (mode shapes, frequencies, generalised masses and 
r.l „«,>„«> .-or th. t, b. ^ 

This information is required to locate sensors and control f. ,M * " f "‘ 
solutions to the problems Identified in the earlier ssslyses. The deaita 
htllty or the need to use reaction Jets must be determined, tfkitg < 
ace . ml he fact that within the atmosphere the let flow Induces s secondary 
prelu-n field on *he aerodynamic surface from which «. emits Tuts causes 
the force magnitude to become a funet ton of flew regions through whlth 

thi? vrh \ v \ v - BUbaotHt , t , mipnraon tv , atul hyprrmm o. 

As the problem definitions are sharpened and an active contwl system 
dealgu evolves, interface conflicts with the vehicle configuration and 
tlZ ground ntles will arise. One example might he the desire o have 
minimum penetration of the orblter lower surfaces for heat protection 
reasona versus the need for positive and negative forces from a given 
control surface, Trade-off studies may be required to determtne whethei 
adding a separate aerodynamic surface as a force generator is belle, than 
using existing surfaces. 

Finally, detail >d tests will be needed to verify that the aerodynamic 
character 1st lea of the seleeted force generators "11 provide the de* red 
inputp. Also s development program should be Initiated to tnsu.e tbs 
sc tuators of sufficient capacity and redundancy will be available. 
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INITIATE DEVELOPMENT OF MULTIPLE-REDUNDANCY ACTUATORS 


